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SECTION  I 


INTRODUCTION 


1.  OBJECTIVES 

The  objectives  of  this  investigation  are  (1)  to  develop  a  teclmique  for 
analyzing  doubly  corrugated  steel  shells  used  as  protective  shelters  for  air¬ 
craft,  and  (2)  to  verify  the  analysis  by  full-scale  testing  of  such  a  structure. 

2.  DESCRIPTION  Of  STRUCTURE 

Although  the  analysis  developed  and  presented  in  this  report  is  applicable 
to  general  shell  structures  exhibiting  nonlinear  response  when  subjected  to  ar¬ 
bitrary  static  loading,  the  particular  structure  studied  here  was  a  semicircular 
steel  arch  shell  constructed  of  doubly  corrugated  arch  panels  bolted  together 
to  form  a  completely  self-supporting  structure.  These  panels  were  fabricated 
of  14-gage  sheet  steel  that  was  first  deformed  into  somewhat  of  a  U-shape  14 
in,  deep.  Then,  to  obtain  the  proper  curvature,  small  corrugations  were 
rolled  into  the  panels  on  the  bottom  and  sides  of  the  U-shape.  Consequently, 
the  assembled  structure  had  deep  corrugations  circumferentially  and  small 
corrugations  longitudinally.  The  structure  analyzed  and  tested  in  this  study 
was  74  ft  long  and  had  an  inside  radius  of  24  ft.  This  structure  was  manu¬ 
factured  by  Wonder  Trussless  Building,  Inc.  of  Chicago,  Illinois. 

5.  PREVIOUS  WORK 

Since  the  structure  being  studied  was  one  to  which  no  closed  form  mathe¬ 
matical  solution  was  applicable,  the  few  previous  attempts  of  analysis  (refs. 
1,2)  have  necessarily  been  restricted  to  treating  the  problem  as  being  a 
linear,  elastic,  isotropic  arch.  To  the  author's  knowledge,  the  change  in 
structural  response  due  to  the  small  corrugations  lias  never  been  studied. 

Similar  structures  have  been  evaluated  in  the  past.  The  U.S.  Naval  Civil 
Engineering  Laboratory  (NCEL)  lias  published  two  technical  reports  on  the 
evaluation  of  doubly  corrugated  metal  buildings  (refs.  5,4).  In  both  of  these 
reports,  the  characteristics  determined  were  shipping  weight  and  cube,  weight 
of  metal  required,  structural  strength,  ruggedness,  erection  and  reerection 
time,  first  costs,  and  maintenance  costs. 
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Although  recently  constructed  aircraft  shelters  have  a  concrete  cover, 
originally  they  were  designed  with  an  earth  cover.  A  rather  comprehensive 
structural  test  program  (ref.  5)  was  conducted  by  the  Civil  Iingineering  Divi¬ 
sion  of  the  Air  force  Weapons  Laboratory  (AFWL) .  The  objectives  of  that  study 
were  (1)  to  establish  design  criteria  and  standards  for  a  fanily  of  protective 
shelters,  aid  (2)  to  evaluate  currently  available  structures  for  their  poten¬ 
tial  use  as  aircraft  protective  shelters.  The  results  indicated  that  a  48-ft- 
diameter  10 -gage  structure  with  a  14-in. -deep  corrugation  was  capable  ol  sus¬ 
taining  a  10- ft  static  load  of  earth  cover. 

4.  APPROACH 

a.  Numerical  Analysis 

A  numerical  technique  capable  of  treating  the  problems  of  nonlinear 
material  properties,  nonlinear  geometry,  and  orthotropic  response  to  arbitrary 
loading  was  developed.  To  accomplish  this,  a  mathematical  model  which  exhi¬ 
bited  the  same  characteristics  (moment- rotation,  force-defoxmation,  and  ulti¬ 
mate  strength)  as  a  finite  portion  of  the  actual  structure  was  developed.  The 
entire  structure  was  then  analyzed  by  appropriate  assembly  of  hundreds  of 
these  models  within  a  computer  program  (appendix  I)  which  was  developed  as 
part  of  this  investigation.  Basically  the  numerical  approach  employed  was 
as  follows: 

(1)  A  differencial  approach  to  the  basic  displacement  method 
was  employed. 

(2)  A  three-parameter  moment- rotation  (and  force- deformation) 
equation  was  used  to  represent  the  material  nonlinearity. 

(5)  A  bean  model  was  used  which,  when  coupled  to  five  other 
similar  bean  models,  represented  the  behavior  of  a  portion 
of  the  shell  both  in  the  elastic  and  plastic  ranges. 

(4)  A  fourth-order  Runge-Kutta  integration  scheme  was  employed 
to  integrate  the  system  of  differential  equations. 

(5)  The  compatibility  matrix  was  adjusted  during  the  incremen¬ 
tal  loading  in  order  to  account  for  geometric  nonlinearity. 

livery  attempt  was  made  to  retain  the  physical  characteri  sties  of  the  real 
structure  while  developing  a  method  of  analysis  which  had  a  tractable  solu¬ 
tion.  'Hie  end  result  was  a  technique  which  predicts  the  behavior  of  the 
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structure  under  any  static  loading  condition.  This  approach  was  similar  to 
that  presented  in  reference  6. 

b.  Experimental  Program 

The  experimental  phase  of  this  investigation  consisted  of  two  parts. 

(1)  Laboiatory  tests  were  conducted  to  determine  the  moment- 
rotation,  force-deformation  relationships,  and  the  fail¬ 
ure  modes  of  finite  portions  of  the  structure.  This  in¬ 
formation  was  then  used  as  input  to  the  computer  program 
previously  mentioned. 

(2)  Full-scale  tests  were  conducted  on  74  linear  feet  of  the 
structure  using  known  loads  in  order  to  verify  the  numer¬ 
ical  analysis. 

As  a  result  of  this  experimental  investigation,  a  technique  which 
can  be  used  to  predict  the  elastic  moment- rotation  response  of  the  individual 
panels  was  developed.  The  significance  of  this  is  best  realized  when  a  re¬ 
quirement  for  the  analysis  of  a  deeper  corrugated  section  arises.  By  testing 
circumferential  strips  of  the  panels  in  uniaxial  tension,  the  variation  of 
the  modulus  of  elasticity  within  the  cross  section  of  the  panels  can  be  deter¬ 
mined.  This  information  can  be  used  to  develop  a  transformed  section  similar 
to  the  approach  used  in  the  analysis  of  reinforced  concrete  members. 


SECTION  II 


FORMULATION  OF  CONCEPTS 

1.  GENERAL  APPROACH  TO  FULL-SCALE  STRUCTURE 

Although  the  analysis  of  the  structure  being  studied  appears,  at  first, 
to  present  a  relatively  simple  problem,  closer  examination  reveals  several 
complications  that  prohibit  the  use  of  existing  shell- analysis  techniques. 

Hie  s tincture  exhibits  pronounced  orthotropy,  nonlinear  response  to  applied 
loads ,  and  can  be  loaded  to  such  an  extent  as  to  produce  geometric  nonlinear¬ 
ity.  Since  the  spectrum  of  loads  of  interest  encompasses  virtually  all  possible 
magnitudes,  the  following  conditions  must  be  studied: 

(a)  elastic  analysis  (loads  of  low  intensity) , 

(b)  nonlinear  analysis  (loads  of  higher  intensity),  and 

(c)  ultimate  strength  anal>sis  (loads  that  cause  collapse). 

A  procedure  has  been  developed  which  treats  all  the  above  conditions 
simultaneously.  So  that  one  may  better  understand  the  approach  used  in  the 
numerical  analysis,  development  of  the  matrix  formulation  is  presented  in 
detail.  The  differential  approadi  used  to  consider  material  nonlinearity  and 
the  incremental  approach  used  to  handle  geometric  nonlinearity  are  presented. 

The  mathematical  models  for  the  nonlinear  beam  element,  the  orthotropic  shell 
element,  and  the  full-scale  structure  are  also  presented. 

2.  MATRIX  FORMULATION 

Indeterminate  structural  systems  which  behave  elastically  can  be  analyzed 
by  either  tire  "force"  or  "displacement"  method.  Nonlinear  systems  can  also 
be  solved  by  either  method  as  presented  by  Richard  and  Goldberg  (ref.  7)  and 
Goldberg  and  Richard  (ref.  8).  The  displacement  method  prevails,  however, 
because  of  the  ease  of  generating  tire  required  equations  on  the  digital  compu¬ 
ter;  consequently,  this  method  is  applied  to  the  approach  taken  in  this  report. 
Applying  tire  displacement  method  to  nonlinear  systems,  using  a  differential 
point  of  view,  results  in  a  set  of  simultaneous  nonlinear  ordinary  differen¬ 
tial  equations  which  is  solved  by  numerical  integration. 
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The  matrix  formulation  of  the  displacement  method,  as  used  in  elastic 
analysis,  is  based  on  equilibrium,  stress-strain,  and  compatibility  relation¬ 
ships  of  a  deformable  body. 

Consideration  of  equilibrium  results  in  tire  following  expression  relating 
tire  vector  of  applied  forces,  (f } ,  to  the  vector  of  internal  forces,  {pi . 

{f}  =  [A] (p)  (1) 

where  [A]  is  known  as  the  equilibrium  matrix. 

Employing  stress-strain  (e.g.,  moment- rotation  or  force- deformation)  con¬ 
cepts,  tire  internal  forces,  {p},  are  expressed  in  terms  of  the  internal  defor¬ 
mations,  {v}. 

(p>  =  [k] {v}  (2) 

where  [k]  is  the  stress-strain  relationship  (usually  called  the  element  stiff¬ 
ness  matrix) . 

Enforcing  continuity  results  in  the  following  expression  relating  exter¬ 
nal  deformations,  {d},  to  tire  internal  deformations,  (v). 

{v}  =  £B] {d}  (3) 

where  |B]  is  tire  compatibility  matrix  presented  in  figure  1. 

Applying  a  set  of  virtual  displacements  to  a  structure  which  is  subjected 
to  tire  action  of  real  loads  and  introducing  the  theorem  of  virtual  work  gives 

(f }  t(d}  =  {p}*"{v}  (4) 

Substitution  of  eqs.  (1)  and  (3)  into  (4)  results  in 

{p}t[A]t{d}  =  {p}t[B]{d}  (5) 

Therefore,  the  compatibility  matrix  is  the  transpose  of  the  equilibrium  matrix. 


[B]  =  [A]C 


(0) 
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Figure  1.  Compatibility  Matrix  Development 


It  can  now  be  shown  that 


If}  =  [B] t[k] £BJ {d} 


{f}  =  [K]{d> 
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and 


{p}  =  [k] [B] {d}  (9) 

where  [K]  is  known  as  the  structure  oriented  stiffness  matrix  ([K]  = 

[B] C[k] [B]). 

Equation  (8)  represents  a  set  of  simultaneous  algebraic  equations,  the 
solution  of  which  yields  the  linearly  elastic  response  of  the  structure  for 
any  statically  applied  loading. 

5.  DIFFERENTIAL  APPROACH 

In  structural  systems  which  exhibit  nonlinear  material  behavior,  the 
displacements  are  not  linear  functions  of  the  applied  loading  and,  consequent¬ 
ly,  a  differential  approach  must  be  taken.  In  order  to  formulate  the  funda¬ 
mental  equations  of  a  nonlinear  structure,  consider  first  the  set  of  equations 
derived  for  linear  structures  and  introducing  proportional  loading, 

X±P  -  ^(dp  (10) 

» 

where  P  is  the  parameter  of  applied  loads,  X^P  represents  the  generalized 
applied  load  at  point  i,  Ki  denotes  the  force -displacement  function  of  the 
structure  at  point  i,  and  d^  is  the  generalized  displacement  at  point  i  in 
the  j  direction.  The  differential  form  of  eq.  -(10)  is 


or 


A. 

i  = 


F. 

i 


(d.  ,d.) 
J  J 


where 


F. 

i 


<VV 


dK.  (d.) 
-4  J 

dP 


(11) 


(12) 


(13) 
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and 


*1  -  3^9  04) 

Integration  of  eq.  (14)  results  in  a  complete  displacement  pattern  for 
the  nonlinear  structure  for  any  given  applied  loading  and  initial  conditions 
(ref.  8).  This  integration  is  accomplished  numerically  using  a  fourth  order 
Runge-Kutta  technique. 

By  introducing  a  shorthand  notation  for  the  first  term  on  the  right  side 
of  eq.  (13), 


(15) 


K*  is  the  generalized  differential  stiffness  matrix  and  is  determined  by 


K*  =  BHcB  (16) 

where  k  is  a  square  matrix  relating  differential  element  displacements  to  the 
differential  element  forces. 

4.  NONLlNliAR  NDDliL 

In  order  to  develop  a  mathematical  beam  model,  an  expression  for  nonlinear 
moment- rotation  and  force -deformation  is  required.  The  Richard  stress-strain 
formula  (ref.  7)  was  chosen  because  of  its  smooth  transition  to  fully  plastic 
behavior  and  its  asymptotic  nature.  'Ihe  following  equation  is  one  form  of 
the  expression. 


M  = 


R6 


1  4* 

R0 

rf 

M  1 

_ 

P1 

1/n 


(17) 


where  R,  0,  M,  M^,  and  n  denote  the  stiffness  coefficient,  rotation,  moment, 
ultimate  moment,  and  the  curve-fitting  parameter,  respectively.  A  nondimen- 
sionalized  form  of  this  equation  is  shown  in  figure  2. 
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moment,  force,  stress 
ultimate  (or  plastic)  value  of  T 

rotation,  deformation,  strain 
teraiinal  value  of  E  corresponding  to  Tq 

curve-fitting  parameter 


Figure  2.  Richard  Equation 


where  [K]"1  is  the  flexibility  matrix  corresponding  to  the  stiffness  matrix, 
[K] ,  and  [K*]’1  is  the  corresponding  differential  relationship  which  is 
constant  during  each  load  increment  (as  well  as  during  each  numerical 
integration  increment). 

By  using  eqs.  (19)  and  (22),  the  following  expression  can  be  written  for 
a  single  beam  element. 
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(23) 


Introducing  the  notations 


M. 

n 

1  - 

i 

M 

P 

(n  +  l)/n 


and 


► 

M 

n 

A.  = 

1  - 

_L 

j 

M 

. 

P 

_ 

(n  +  l)/n 


and  solving  cq.  (.23)  i'or  the  differential  moments,  the  following  equation  is 
obtained. 


r  *\ 

dM. 

i 

Hp~ 

6EI 

L 

2A. 

i 

jv, 

< 

r-  ■ 

dM 

HP1 

< 

< 

1 

l _ 

A. A. 

_  1  J 

i 

i 

!  2A. 

1  3  _ 

d0 . 

dP1 

V  y 

(24) 


To  illustrate  the  behavior  of  tire  model,  consider  a  large  value  for  the 
curve- fitting  parameter,  n  (thus  making  the  curve  bilinear).  If  neither  end 
moment  has  readied  its  ultimate  value  (Mi  <  Mp  and  Nf  <  M  ) ,  eq.  (24)  becomes 


which  indicates  linear  action  for  the  beam.  If,  however,  the  moment  at  the  j 

end  of  the  bean  element  has  reached  its  ultimate  value  (ML  =  M  ) ,  while  the 

j  P 

moment  at  the  i  end  is  still  elastic  (Mi<  M^),  eq.  (24)  becomes 


which  represents  the  stiffness  coefficients  for  the  beam  pinned  at  the  j  end. 
This  nonlinear  beam  model  is  used  in  the  assembly  of  the  shell  model,  which 
consists  of  six  interconnected  beams. 

In  order  to  incorporate  the  effect  of  axial  force- deformation  in  the 
shell,  tne  element  stiffness  matrix  must  be  expanded.  This  expansion  is 


where  Vj,  f2,  v2,  f3,  v3  denote  the  axial  force  and  deformation,  moment 
aid  rotation  at  the  i  end,  and  moment  and  rotation  at  the  j  end,  respectively. 
The  term  Ar  represents  the  cross-sectional  area  of  the  element  and  f0  denotes 
the  ultimate  axial  force  for  each  element. 
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Figure  3  shows  a  portion  of  the  actual  structure  to  be  analyzed,  an 
equivalent  physical  model,  a  shell  model,  and  a  beam  model.  Since  the  back¬ 
ground  for  tiie  nonlinear  beam  model  has  been  established,  it  remains  now  to 
show  how  these  beam  elements  are  to  be  interconnected  to  form  the  shell  model 
and  to  develop  expressions  for  the  physical  properties  of  these  elements. 

Geometric  nonlinearity  is  treated  with  an  incremental  loading  approach 
that  assumes  finite  deflections.  If  large  deflections  are  to  be  considered 
in  the  future,  an  iterative  approach  must  be  used.  However,  a  rough  approxi¬ 
mation  of  a  collapsed  configuration  may  be  obtained  by  incrementing  the  load 
vector. 

5.  SHBLL  MODBL 

In  order  to  obtain  a  general  method  for  the  analysis  of  nonlinear  shell 
structures,  a  model  must  first  be  developed  which  adequately  predicts  the 
initial  linear  response  and  which  may  be  modified  to  account  for  nonlinear 
action.  Hie  framework  model,  comprised  of  six  interconnected  beam  elements, 
appears  to  fulfill  the  above  requirements.  A  model  first  introduced  by 
Hrennikoff  (ref.  9) ,  later  used  by  Yettram  and  Husain  (ref.  10) ,  and  subse¬ 
quently  extended  to  the  nonlinear  analysis  of  plate  structures  by  Smith  (ref. 
6)  was  used  as  a  guide  to  develop  a  model  for  orthotropic  shell  structures. 

The  model  used  in  this  study  is  shown  in  figure  4.  In  order  to  obtain  the 
properties  of  the  individual  beams,  laboratory  tests  were  conducted  to  deter¬ 
mine  the  initial  elastic  response  and  the  ultimate  (plastic)  capacity  of 
portions  of  the  actual  structure  when  subjected  to  bending  and  axial  force. 

The  details  of  this  technique  and  the  laboratory  results  are  presented  in 
section  III. 

The  full-scale  structure  is  modeled  by  assembling  (on  the  digital  compu¬ 
ter)  many  elements;  consequently,  any  load,  symmetric  or  unsymmetric,  may  then 
be  applied  and  the  results  will  include  the  elastic  and  nonlinear  responses. 
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SECTION  III 


DETERMINATION  OF  MODEL  CHARACTERISTICS 


1.  OBJECTIVE 

Before  the  concepts  previously  developed  can  be  used  to  predict  the 
behavior  of  the  entire  structure,  information  concerning  the  action  of  por¬ 
tions  of  tiie  structure  must  be  obtained.  This  section  is  devoted  to  a 
presentation  of  the  laboratory  tests  conducted  which  are  necessary  for  the 
detennination  of  the  required  parameters  for  moment- rotation  and  force - 
deformation  relationships. 

2.  MOMENT- ROTATION  TESTS 

a.  Apparatus 

In  order  to  obtain  the  flexural  behavior  of  a  portion  of  the  struc¬ 
ture,  a  20- ft  section  of  the  arch  was  tested  as  shown  in  figure  5.  This  test 
setup  consisted  of  center  supports  located  76  in.  apart  which  were  connected 
to  the  panel  with  flanged  cartridge  ball  bearings  and  a  2- inch- diameter  solid- 
steel  shaft.  The  loads  were  applied,  near  the  ends  with  double-acting  hydrau¬ 
lic  rams  which  made  it  possible  to  subject  the  panel  to  both  positive  and 
negative  moments.  Figure  6  shows  the  positive  moment  setup.  It  should  be 
noted  that  the  base  of  the  ram  was  hinged  to  allow  for  rotation  during  loading. 
In  oi'der  to  subject  the  panel  to  negative  moment,  a  1- inch-diameter  eyebolt 
was  extended  through  the  concrete  (used  for  a  load- distributing  device)  and 
the  steel.  The  hydraulic  pressure  line  was  connected  to  the  top  of  the  ram. 

Although  figure  S  shows  the  lower  end  of  both  supports  as  fixed, 
later  modifications  involved  the  installation  of  a  shaft  with  pillow  block 
type  ball  bearings  so  the  effect  of  combined  bending  and  axial  force  could  be 
studied.  This  effect  was  negligible  on  the  initial  response  but  did  tend  to 
result  in  a  change  in  the  ultimate  moment  of  the  section  of  about  15  percent 
(stronger  with  fixed  base). 

b.  Instrumentation 

Since  the  three  parameters  needed  to  use  eq.  (17),  shown  in  figure 
2,  are  initial  slope  of  the  moment- rotation  curve,  the  ultimate  moment,  and 
the  overall  shape  of  the  curve,  instrumentation  was  chosen  that  would 
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continuously  record  the  pressure  in  the  rams  and  the  resulting  rotation  at  the 
supports .  This  information  was  recorded  on  the  existing  data  acquisition 
system  at  CERF.  It  was  recorded  at  1-7/8  in. /sec,  played  back  at  60  in. /sec, 
and  plotted  on  an  x-y  plotter.  Figure  7  shows  the  method  used  to  obtain  the 
amount  of  rotation.  A  12-inch-diameter  aluminum  plate  was  attached  to  the 
support  shaft;  piano  wire  was  placed  around  the  circular  plate  and  attached  to 
a  weight  which  moved  vertically  under  test  conditions.  This  movement  was 
measured  with  linear  potentiometers  on  both  sides  of  the  shaft  in  order  to 
cancel  the  effect  of  vertical  movement  of  the  shaft. 


Figure  7.  Rotation  Measurement  Instrumentation 
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c.  Results 


Initially,  it  was  believed  that  there  would  be  significant  differ¬ 
ences  in  positive  and  negative  bending.  However,  inspection  of  the  results 
(shown  plotted  in  figure  8)  revealed  the  difference  in  elastic  response  to  be 
relatively  small  and,  consequently,  the  effort  required  to  incorporate  the 
difference  into  the  computer  program  was  not  considered  necessary.  The  major 
difference  occurred  in  the  ultimate  moment. 

The  initial  slope  of  the  moment- rotation  curve  is 
M  . 

■q  =  2,000  ft- kips /radian  (28) 


and  the  ultimate  moment  is 


Mq  =  30  ft-kips/panel  (29) 

Although  the  curve-fitting  parameter,  n,  is  approximately  3,  values  of  2  and  4 
would  not  be  unreasonable.  Studies  were  conducted  with  the  computer  program 
to  determine  the  difference  in  final  results  for  a  range  of  2  to  4  on  the 
curve- fitting  parameter.  The  final  results  were  only  slightly  affected  for 
problems  where  structural  yielding  was  not  widespread. 

In  every  test,  failure  of  the  panels  was  of  the  buckling  type. 

Figure  9  shows  typical  types  of  failure. 

3.  FORCE- DEFORMATION  TESTS 

a.  Apparatus 

As  in  the  case  of  moment- rotation  tests,  axial  force -deformation 
tests  were  conducted  to  determine  the  initial  response,  the  ultimate  axial 
force,  and  the  general  shape  of  the  force- deformation  curve  for  a  portion  of 
the  structure.  The  equipment  used  in  this  test  was  a  static  hydraulic  testing 
machine.  The  test  setup  is  shown  in  figure  10.  Figure  11  shows  an  early  trial 
test  on  a  short  portion  of  the  structure  (which  proved  unstable);  it  is  pre¬ 
sented  here  because  the  bottom  support  in  subsequent  tests  was  the  same  as 
tiie  top  support  shown  in  this  figure. 
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Moment,  ft-kips 


Figure  9.  Typical  Buckling  Failures 
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Figure  10.  Force-Deformation  Test  Apparatus 


b.  Instrumentation 

In  order  to  obtain  the  complete  shape  of  the  force -deformation  curve, 
both  the  load  from  the  hydraulic  loader  and  the  displacement  of  the  loading 
head  were  continuously  recorded  on  an  x-y  plotter.  In  addition,  potentio¬ 
meters  were  mounted  horizontally  to  observe  lateral  displacement  (fig.  12). 

This  information  was  recorded  oilpaper  tape  readout. 

c.  Results 

Since  the  panels  were  initially  curved,  the  axial  loads  created  some 
flexure  because  of  the  eccentricity.  This  effect  was  minimized  by  aligning 
the  end  loads  with  the  centroid  of  the  center  (lengthwise)  of  the  panel.  Vir¬ 
tually  all  of  the  sections  tested  failed  by  local  buckling  at  or  very  near  the 
center.  The  type  of  buckling  was  essentially  the  same  as  that  resulting  from 
the  moment -rotation  tests.  The  response  of  the  panels  when  subjected  to  axial 
force  is  presented  in  figure  13.  Several  specimens  were  loaded  and  unloaded 
repeatedly  and  a  typical  curve  for  this  test  is  shown  in  figure  14.  The  re¬ 
peated  loading  had  no  effect  on  either  the  initial  or  the  final  response  of 
the  panels. 

Frcm  figure  13,  the  average  initial  slope  of  the  curves  for  tests 
A-2,  A-5,  A-5,  and  A-6  is 

=  1,950  kips/in.  (30) 

and  the  ultimate  force  is 

f  =  40  kips/panel  =  1.665  kips/in.  (31) 

Also,  the  curves  are  essentially  bilinear  and  can  be  reproduced  with  eq.  (17) 
using  a  large  value  of  n  (say  n  =  20) .  These  values  are  needed  to  develop 
a  model  with  the  same  characteristics. 

In  every  test,  the  failure  was  a  buckling  type  similar  to  those 
shown  in  figure  15. 
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Figure  12.  Lateral  Displacement 

Measurement  Instrumentation 


4.  MODEL  CHARACTERISTICS 

In  order  to  generate  a  system  of  equations  which  represents  the  action  of 
the  entire  structure  when  subjected  to  an  arbitrary  loading  condition,  infor¬ 
mation  must  be  obtained  regarding  the  representative  equations  for  a  portion 
of  the  structure.  Basically,  the  parameters  in  the  equations  shewn  in  matrix 
form  by  eq.  (27)  must  be  detemined. 

First,  the  actual  structure  is  reduced  to  a  cylindrical  shell  which  has 
a  constant  thickness  circumferentially  based  on  the  initial  elastic  response 
of  portions  of  the  structure.  Also,  the  model  will  have  a  constant  thickness 
in  the  longitudinal  direction  which  differs  from  the  circumferential  thick¬ 
ness.  It  should  be  noted  that  the  model  will  have  characteristics  that  re¬ 
present  both  compression  and  flexural  modes. 
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Figure  13.  Force -Deformation  Curves 


Consider  the  flexural  mode  in  the  circumferential  direction.  If  moments 
are  applied  along  parallel  sides  of  the  element  (shown  in  figure  16a) ,  the  re 
suiting  rotation  is  ML/ 2E I  if  purely  cylindrical  bending  is  assumed.  Solving 
in  terms  of  M/6 


I 


Figure  16.  Models 


liquating  this  to  the  experimental  value,  eq.  (28) 


thus 


M  21il  2(30  x  103)  1 

0  "  1)  75 


2,000(12) 


(33) 


I  =  50.4  in.1* /panel 


(34) 


A  very  important  point  must  be  made  here.  In  all  previous  analyses,  the 

moment  of  inertia  has  been  taken  as  I  =  J^y2dA  which  results  in  1  =  92.4  in.1*/ 

panel.  This  approach  assumes  a  homogeneous  cross  section  which  at  first  may 
appear  valid  but  if  the  small  corrugations  are  considered,  it  becomes  obvious 
that  the  slopes  of  the  stress  strain  curves  vaiy  within  the  cross  section. 

It  should  be  noted  that  the  stress  and  strain  are  average  values  and  not  the 
true  stress  and  strain  at  a  point.  The  author  has  developed  a  procedure  to 
predict  the  elastic  response  of  the  section  by  removing  (cutting)  strips  from 

the  section  and  obtaining  a  modular  ratio  which  is  then  used  to  develop  a 

transformed  section.  The  results  of  this  approach  are  essentially  the  same 
as  those  obtained  from  the  moment-rotation  tests. 


Consider  a  rectangular  section 


I  = 


1 

12 


at 


3 

efc 


(35) 


where  t  denotes  the  equivalent  thickness  for  flexure  in  the  circumferential 
direction.  For  our  case  where  a  =  24  in.  (single  panel  test) 

t  fc  =  2.48  in.  (36) 

It  should  be  noted  that  other  approaches  are  possible.  For  example,  an  effec¬ 
tive  El  could  be  deteimined.  However,  since  E  appears  in  both  the  axial  com¬ 
pression  and  the  flexural  modes,  it  was  decided  that  the  approach  used  here  is 
a  more  logical  one. 

Since  the  structure  is  orthotropic,  an  investigation  into  its  flexu¬ 
ral  capacity  in  the  longitudinal  direction  is  required.  The  strength  (both  in 
compression  and  flexural  modes)  in  this  direction  is  very  small  compared  to 
that  in  the  circumferential  direction.  In  fact,  tests  have  shown  that  the 
strength  in  this  direction  could  be  neglected  if  the  grid  size  could  be  chosen 
small  enough.  However,  because  of  the  core  storage  on  the  digital  computer, 
the  grid  size  used  was  such  that  the  flexural  capacity  in  the  longitudinal 
direction  was  required. 

Consider  a  portion  of  the  structure  as 
location  of  the  centroid  is  detexmined  as 

/  ydA  /^Rcos9tRd6 
7  =  T^dT  =  2  2W  = 

The  moment  of  inertia  about  the  x-axis  is 

I  =  /  y2dA  =  2/<f,(Rcos0)2tRd9  =  R3t(<J>  +  |  sin2<J>)  (38) 

X  A  0  2. 

Then  the  moment  of  inertia  about  the  centroidal  axis  is 

JCG  =  \  '  Ay2  =  +  2  Sin2<1’  '  f  sinZV  (39) 

Tins  information  is  used  in  appendix  I  to  determine  the  properties  of  longi¬ 
tudinal  members. 


shown  in  figure  17.  The 
fsin*  (37) 
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Figure  17.  Section  of  Structure 


'Hie  force- deformation  relationships,  as  determined  from  laboratory 
tests  and  shown  in  figure  13,  are  used  to  develop  the  required  model  param¬ 
eters.  To  obtain  an  equivalent  shell  thickness,  consider 


_  a  _  P/A  _  /JP\  (_L\ 
~  e  ~  AI7l  ~  \&L)  \A*  / 


(40) 


/P\/M__(1?950)(76) 
\ AL  /  \  E  /  3C ,000 


0.495  in?/panel 


(41) 
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Then  the  equivalent  thickness  for  axial  stress  in  the  circumferential  direc¬ 
tion  is 


t  =  0.493/24  =  0.0205  in.2/in.  (42) 

'Hie  equivalent  thickness  in  the  longitudinal  direction  was  taken  as  the  true 
material  thickness. 

Considering  the  ultimate  force,  eq.  (31),  and  the  equivalent  area, 
eq.  (41),  the  ultimate  axial  stress,  o  ,  is  obtained. 

%  =  A*  =  0. 493  =  81,3  ksi  (-43') 

'lhe  diagonal  members  shown  in  figure  16b  are  not  actually  required 
for  stability  in  low  ranges  of  loading.  However,  as  the  structure  starts  to 
yield,  the  equations  become  unstable  and  their  solution  meaningless.  These 
members  are  supplied  to  stabilize  the  system  and  permit  a  complete  nonlinear 
solution. 
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SECTION  IV 


FULL-SCALE  TESTS  AND  ANALYSIS 

1.  CONSTRUCTION 

Since  there  was  initially  a  possibility  of  covering  the  structure  with 
either  soil  or  concrete,  a  loading  system  was  designed  with  these  eventualities 
in  mind.  The  loading  floor  was  designed  to  serve  not  only  as  a  de adman  for 
loading  the  structure  but  also  as  a  system  that  could  be  used  to  conduct 
moment- rotation  tests  on  portions  of  the  structure  with  varying  thicknesses 
of  concrete  cover.  This  reinforced  concrete  floor  was  16  ft  wide,  32  ft  long, 
and  8  ft  thick  with  screw  type  inserts  for  tiedowns.  To  minimize  support 
settlement  and  spreading,  wing  walls  and  horizontal  tie  bars  were  included. 
Figure  18a  shows  the  formwork  and  reinforcing  steel  used  in  the  loading  floor. 
The  concrete  was  placed  in  2-ft  lifts  with  the  first  three  lifts  being  2,500- 
psi  concrete  and  the  top  lift  4,000-psi  concrete.  An  overall  view  of  the 
loading  floor,  wing  walls,  and  footings,  prior  to  backfilling,  is  shown  in 
figure  18b.  Figure  18c  shows  in  greater  detail  the  screw  anchor  inserts. 

These  inserts  are  on  2-ft  centers  and  each  one  has  a  "pullout"  strength  of 
approximately  80  kips.  The  completed  loading  floor  after  backfilling  is  shown 
in  figure  18d. 

To  determine  the  best  method  of  erection,  single,  double,  and  triple 
arches  were  assembled  on  the  ground  and  lifted  into  place  with  a  crane.  The 
double- arch  procedure  proved  to  be  the  best  and  was  consequently  used  in  sub¬ 
sequent  construction.  Since  tests  were  plained  for  three  different  shelter 
lengths,  two  arches  were  omitted  during  construction.  Figure  19a  shows  the 
structure  partially  completed  and  figure  19b  shows  the  structure  ready  for 
tiie  first  series  of  tests. 

2.  LOADING  AND  INSTRUMENTATION  SYSTEMS 

The  loading  system  was  constructed  by  attaching  double-acting  rails  to  the 
loading  floor.  Force  links  were  then  connected  from  the  rams  to  the  cables  by 
meais  of  a  pin  at  one  end  and  a  free- turning  swivel  at  the  other  end.  The 
cables  were  then  extended  to  the  structure  and  connected  to  U-nuts  which  were 
attached  to  1-inch-diameter  high-strength  full- threaded  bolts.  The  bolts 
passed  through  both  the  metal  aid  short  sections  of  4-  x  4- in.  timber  which 
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were  partially  encased  in  metal  channels.  The  tops  of  the  bolts  were 
then  double -nutted.  This  arrangement  is  shown  in  figure  20.  Hydraulic  pres¬ 
sure  was  supplied  to  the  rams  by  means  of  the  air- actuated  pump  shown  in 
figure  21. 

Although  the  force  in  the  cables  could  have  been  approximated  by  record¬ 
ing  the  hydraulic  fluid  pressure,  the  use  of  strain  gaged  force  links  provided 
a  much  more  accurate  record  of  the  force  in  the  individual  cables.  Figure  21 
shows  the  two  10- channel  switch  and  balance  units  and  the  portable  strain 
indicator  that  were  used  to  determine  the  cable  forces.  The  force  links  were 
calibrated  so  that  1  pin. /in.  of  strain  represented  10  lb  of  force. 

Horizontal  and  vertical  deflections  were  measured  at  35  locations  on  the 
structure  as  shown  in  figure  22.  This  measuring  method  consisted  of  two  hori¬ 
zontal  strands  of  piano  wire  (vertical  separation  of  approximately  2  ft) 
stretched  perpendicular  to  the  longitudinal  axis  of  the  structure  at  five  lo¬ 
cations.  These  wires  were  attached  to  standards  which  were  independent  of  the 
structure.  Thirty-five  weights  were  then  suspended  from  the  structure  with 
piano  wires  in  such  a  way  that  each  one  intersected  with  a  pair  of  the  hori zon- 
tal  wires.  Short  strips  of  measuring  tape  (1/16-in.  markings)  were  attached 
to  the  vertical  wires  and  to  one  horizontal  wire  in  each  pair.  This  arrange¬ 
ment  (fig.  23)  provided  good  records  when  read  with  a  transit. 

3.  'ATliS  OF  TliSTS 

Since  one  question  to  be  answered  in  tills  research  endeavor  was  whether 
the  structure  exhibited  primarily  arch  or  shell  action,  three  series  of  tests 
were  conducted  on  different  structure  lengths.  The  first  seven  tests  were 
made  on  30  linear  ft  of  the  shelter  with  the  loads  distributed  throughout  the 
entire  length.  The  second  seven  tests  were  conducted  on  52  linear  ft  of  the 
structure  with  loads  applied  on  a  30-ft  portion  from  one  end.  Finally,  seven 
tests  were  made  on  74  linear  ft  of  the  structure  with  loads  applied  on  the 
center  30  ft.  Figure  24  shows  the  arrangement  of  the  loads.  Not  all  hydrau¬ 
lic  rams  were  acting  on  all  tests.  Table  I  gives  the  type  of  tests,  structure 
length,  and  type  and  amount  of  loading  for  the  full-scale  test  series. 
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fa)  Reinforcing  Steel  Layout 


(b)  Loading  Floor  Prior  to  Backfilling 
Figure  18.  Construction  of  Loading  Floor 
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(d)  Loading  Floor  After  Backfilling 


Figure  18 — Concluded 


(b)  Completed  Shelter 
Figure  19.  Shelter  for  First  Test  Series 
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Figure  21.  Load- Producing  and  Load-Measurement  Systems 


Figure  22.  Deflect ion -Measurement  System 
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Figure  24.  Applied  Loads 


Table  I 


FULL- SCALE  TESTS 


Type 

Structure 

Loading 

Load  Numbers  Acting 

Average  Maximum 
Load  per  Ram,  lb 

Total 

Test 

Number 

of 

Test 

Length, 

ft 

(See  fig.  24 
for  load  numbers) 

Side 

Rams 

Center 

Rams 

Load, 

lb 

2 

Symmetric 

30 

1  through  14 

0 

4,956 

69,385 

3 

Symmetric 

30 

All  22 

3,434 

2,940 

68,607 

4 

Unsyinmetric 

30 

1  through  14, 

19  through  22 

3,539 

2,970 

55,690 

S 

Unsymmetric 

30 

1  through  18 

3,598 

2,980 

56,090 

6 

Symmetric 

30 

15  through  22 

5,981 

0 

47,850 

1 

Unsymmetric 

30 

19  through  22 

8,020 

0 

32,080 

8 

Unsymmetric 

30 

15  through  18 

7,962 

0 

31,850 

102 

Symmetric 

52 

1  through  14 

0 

6,044 

84,425 

103 

Symmetric 

52 

All  22 

4,048 

3,450 

80,784 

104 

Unsymmetric 

52 

1  through  14, 

19  through  22 

4,793 

4,010 

75,200 

105 

Unsymmetric 

52 

1  through  18 

4,766 

3,980 

74,645 

106 

Symmetric 

52 

15  through  22 

7,806 

0 

62,450 

107 

Unsymmetric 

52 

19  through  22 

8,130 

0 

32,520 

108 

Unsymmetric 

52 

15  through  18 

8,108 

0 

32,430 

202 

Symmetric 

74 

1  through  14 

0 

6,010 

84,055 

203 

Symmetric 

74 

All  22 

4,158 

5,490 

81,960 

204  • 

Unsymmetric 

74 

1  tli  rough  14, 

19  through  22 

4,783 

4,010 

75,320 

205 

Unsymmetric 

74 

1  through  18 

4,791 

4,020 

75,565 

206 

Symmetric 

74 

15  tli  rough  22 

7,888 

0 

63,100 

207 

Unsymmetric 

74 

19  through  22 

8,175 

0 

32,700 

208 

Unsymmetric 

74 

15  tli  rough  18 

8,088 

0 

52,350 

Tests  were  numbered  as  follows:  the  first  series  (30 -ft  shelter)  is 
numbered  2  through  8;  the  second  series  (52 -ft  shelter)  is  numbered  102 
through  108;  and  the  last  series  (74 -ft  shelter)  is  numbered  202  through  208. 
Consequently,  to  determine  whether  the  structure  behaves  like  an  arch  or  a 
shell,  comparisons  must  be  made  among  tests  ending  in  the  same  number.  To 
facilitate  sudi  a  comparison,  table  II  is  presented  to  show  the  effect  of 
structure  length  on  the  centerline  deflections.  The  last  column  is  a  normal¬ 
ized  load/ deflection  computation  in  which  the  lowest  value  of  P/A  is  taken  as 
the  reference  and  the  other  two  tests  using  the  same  loading  are  then  referred 
to  it.  The  lowest  value  was  taken  so  the  percentage  difference  could  be  auto¬ 
matically  presented.  For  example,  consider  tests  5,  105,  and  205.  With  test 
205  as  the  reference,  it  is  obvious  that  test  105  differs  by  almost  3  percent 
and  test  5  differs  by  .over  6  percent.  Furthermore,  it  can  readily  be  seen 
fran  sudi  comparisons  that  the  structure  exhibits  essentially  ardi  behavior 
and,  consequently,  the  length  effect  is  negligible. 

4.  NUMERICAL  ANALYSIS  AND  COMPARISONS  TO  BXPHRIMliNTAL  RESULTS 

The  numerical  analysis  developed  in  this  investigation  is  based  on  a  dis¬ 
crete  element  idealization  tedinique  which  uses  the  beam  element  developed  in 
section  II.  The  grid  systems  used  for  eadi  of  the  test  cases  were  as  follows: 


Test 

Series 

Test  Case 

Grid  System 

Reference 

1 

5,  8 

20  x  10 

Grid  layouts  are 

1 

2,  3,  6 

10  x  10 

contained  at 

2 

102,  103,  106 

10  x  15 

the  end  of  the 

3 

202,  203,  206 

10  x  20 

report . 

All  test  problems  except  5  and  8  were  symmetrically  loaded  and,  consequently, 
the  grid  represents  only  one-half  the  structure  for  most  cases.  The  grid 
systems  were  chosen  with  as  many  degrees  of  freedom  as  possible  while  utiliz¬ 
ing  only  the  central  core  memory  of  the  computer. 

Since  the  structure  was  constructed  of  doubly  corrugated  steel  panels 

s' 

whidi  have  essentially  no  areas  free  of  small  corrugations  or  holes,  strain 
measurements  were  not  attempted.  Deflection  measurements  were  taken  during 
eadi  test  and  these  results  were  compared  to  the  analytical  predictions  (fig. 
25).  Horizontal  deflections  were  not  recorded  in  test  2  because  the  measuring 
system  had  not  readied  the  final  stage  of  refinement.  However,  in  all  subse¬ 
quent  tests,  horizontal  measurements  were  recorded. 
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Table  II 


CENTERLINE  DEFLECTION  VERSUS  TOTAL  LO.iD  COMPARISONS 


Tes  t 
Number 

Total 

Load  (P) , 
lb 

Centerline 
Deflection  (A) , 
in. 

P/A, 
lb/ in. 

Normalized 

P/A 

2 

* 

69,385 

3.38 

20,528 

1.0159 

102 

a 

70,070 

3.44 

20,369 

1.0080 

202 

a 

69,710 

3.45 

20,206 

1.0000 

3 

a 

61,740 

3.49 

17,690 

1.0000 

103 

* 

61,218 

3.38 

18,396 

1.0399 

205 

* 

62,720 

3.48 

18,023 

1.0188 

4 

t 

55,690 

2.81 

19,819 

1.0104 

104 

t 

56,100 

2.86 

19,615 

1.0000 

204 

t 

56,580 

2.72 

20,801 

! 

1.0295 

5 

t 

56,090 

2.61 

21,490 

1.0635 

105 

t 

56,S80 

2.72 

20,801 

1.0295 

205 

T 

56,575 

2.80 

20,205 

1.0000 

6 

A 

47,850 

2.06 

23,228 

1.0000 

106 

A 

46,610 

1.97 

23,660 

1.0185 

206 

A 

57,048 

2.00 

23,524 

1.0127 

7 

t 

32,080 

1.44 

22,278 

1.0000 

107 

t 

32,520 

1.31 

24,824 

1.1142 

207 

t 

32,700 

1.37 

23,869 

1.07J4 

8 

t 

31,850 

1.50 

21,233 

1.0000 

108 

t 

32,430 

1.37 

23,672 

1.1148 

208 

f 

32,350 

1.38 

23,442 

1.1040 

*  Symmetric  loading, 
i  Unsymmetric  loading. 
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Figure  25- --Continued 
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Figure  25 — Concluded 
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Although  22  full-scale  tests  were  conducted,  only  11  are  presented  here. 
The  first  test  was  considered  a  "shake-down"  and  was  conducted  in  an  attempt 
to  remove  the  effects  of  bolt  misalignment  and  other  irregularities  caused  by 
fabrication  and  construction.  The  results  of  two  of  these  tests  were  disre¬ 
garded  because  of  high  wind  in  one  instance  and  rain  in  the  other.  Because 
of  the  limited  computer  storage  capacity,  analytical  solutions  for  the  un- 
symmetrically  loaded  52-  and  74-ft-long  structures  were  not  obtained.  Con¬ 
sequently,  test  results  from  the  second  and  third  test  series  involving 
unsymmetrically  loaded  structures  will  be  presented  in  the  future. 


Figure  26  is  presented  to  give  the  reader  an  idea  of  the  response  of  the 
crown  when  subjected  to  various  load  configurations.  No  scale  is  shown  since 
these  plots  (from  program  presented  in  appendix  II)  are  for  the  linear  analy¬ 
sis  only  and  essentially  the  same  information  may  be  obtained  from  figure  25. 

It  has  been  noted  that  no  strain  measurements  were  taken  and,  consequently, 
no  experimental  stress  values  can  be  presented.  However,  in  order  to  develop 
an  insight  into  probable  stress  patterns,  an  approximate  stress  analysis  was 
performed  for  a  selected  number  of  the  test  problems.  These  plots,  shown  in 
figure  27,  are  useful  only  for  predicting  regions  of  high  stress  and  are  not 
intended  for  design  purposes.  The  stress  distribution  of  test  problem  3  is 
with  fixed  supports  while  all  others  are  based  on  simple  supports.  The  ques¬ 
tion  of  the  degree  of  fixity  at  the  base  was  studied  and  it  was  concluded  that 
for  the  loading  considered  in  this  investigation,  the  hinged  condition  was 
more  representative  of  the  actual  conditions.  The  main  reasons  for  this  con¬ 
clusion  were  (1)  the  base  plates,  to  which  the  structure  was  attached,  were 
secured  to  the  footings  in  such  a  way  as  to  allow  some  rotation,  and  (2)  the 
panels  were  bolted  to  the  base  plate  in  a  staggered  manner  so  that  the  bolts 
along  the  outside  and  inside  were  2  ft  apart  and  the  holes  were  considerably 
larger  than  the  bolts.  It  is  believed,  however,  that  for  high  load  intensities, 
the  slack  will  be  taken  out  of  the  system  and  the  fixed  base  condition  will  be 
more  representative.  This  presents  a  rather  difficult  problem  if  the  support 
response  becomes  a  function  of  the  load  intensity. 
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(<0  lliird  Test  Series  (74-ft  Structure) 


Figure  26.  Crown  Displacements  (Linear  Analysis) 


Inside  Stresses 


Outside  Stresses 


(a)  Test  Problem  6 


Scale 


(b)  Test  Problem  106 


A 


Figure  27.  Comparative  Stresses  on  Inside  and  Outside  Fibers 


SO 


Inside  Stresses 


Outside  Stresses 


Figure  27 — Concluded 


SECTION  V 


CONCLUSIONS 

The  following  conclusions  have  been  made  as  a  result  of  this  investiga¬ 
tion  of  a  doubly  corrugated  steel  aircraft  shelter. 

(1)  The  initial  response  of  the  structure  can  be  predicted  with  reason¬ 
ably  good  accuracy  by  treating  it  as  a  linear  elastic  problem  with¬ 
out  conducting  moment- rotation  and  force -deformation  tests  on 
individual  panels.  However,  the  cross-sectional  properties  (such 
as  the  moment  of  inertia  and  the  centroid)  cannot  be  taken  from 
existing  manufacturers'  tables  since  this  information  ignores  the 
effect  of  the  small  corrugations.  Rather,  simple  tests  must  be 
made  to  determine  the  apparent  modulus  of  elasticity  of  longitudinal 
strips  at  different  locations  within  the  cross  section.  This  infor¬ 
mation  can  then  be  used  to  develop  a  transformed  section  which 
includes  the  effect  of  the  small  corrugations.  When  a  future 
requirement  for  the  analysis  of  different  doubly  corrugated  sections 
arises,  this  approach  can  result  in  a  significant  savings  since  the 
new  panel  concepts  do  not  have  to  be  fabricated  prior  to  an  initial 

i 

analysis . 

(2)  Ilierc  is  approximately  a  10-  to  15-percent  variation  (depending  upon 
the  loading  conditions)  resulting  from  treating  the  structure  first 
as  an  arch  and  then  as  a  shell. 

(5)  For  the  loadings  used  in  the  full-scale  testing  program,  the  non¬ 
linear  solutions  differed  by  approximately  15  percent  from  the 
linear  solutions.  However,  for  loads  which  could  result  in  inci¬ 
pient  collapse,  the  nonlinear  and  linear  analyses  could  differ  by 
several  hundred  percent. 
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APPENDIX  I 


COMPUTER  PROGRAMS 

1.  INTRODUCTION 

Although  the  computer  program,  SANOS  (Structural  Analysis  of  Nonlinear 
Orthotropic  Shells ),  was  developed  specifically  for  the  analysis  of  the 
structure  described  in  section  I,  it  is  not  restricted  to  such  structures. 

The  program  is  designed  so  that  virtually  any  structural  geometry  (symmetric 
or  unsymmetric)  can  be  handled  with  slight  modifications. 

A  companion  program,  SAN1PLT,  was  written  to  plot  the  results  from  SANOS 
which  expedites  interpretation  of  the  grid  system,  loading,  crown  and  arch 
displacements,  moments,  axial  forces,  and  comparative  stresses  on  inside  and 
outside  fibers. 

Sample  input  and  output  for  the  SANOS  program  are  presented  in  appendix 
III  and  IV,  respectively. 

2.  COMPUTER  PROGRAM  SANOS 

Card  Number (s)  Explanation 

1  Control  card 

2  through  11  Dimension  statements  (These  would  require  changes  to 

increase  the  number  of  elements  used.) 

12  through  18  Common  blocks 

19  through  83  Format  specifications  for  input  aid  output  informa¬ 

tion 

84  through  99  Read  statements  and  comment  cards  for  input  data 

100  through  110  Initialization  of  certain  parameters 

111  through  116  Calculating  dimensions  and  central  angles 

117  through  119  Initialization  of  support  information  vectors 

120  through  129  Reading  and  arranging  support  information 

130  through  144  Initialization  of  parameters 

145  through  146  Reading  load  vector 

147  through  182  Printing  input  information 

183  through  230  Calculating  and  printing  nodal  coordinates 

231  through  280  Calculating  and  printing  element  properties 

281  through  287  Printing  support  and  loading  information 
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288  through  289 
290  through  291 

292  through  430 


431  through  446 

447 

448  through  481 

482  through  511 

512  through  542 

543  through  552 
553  through  562 

563  through  601 

602  through  629 

630  tli  rough  830 


Storing  loading  information  for  later  use 

Calling  subroutine  BAND  to  determine  the  indexing 
matrices  for  banding  the  stiffness  matrix 

This  phase  is  a  DO  loop  that  applies  a  certain  per¬ 
cent  of  the  load,  generates  the  differential  stiff¬ 
ness  matrix,  solves  the  equations,  and  integrates 
the  internal  moments,  axial  forces,  and  displace¬ 
ments.  This  process  is  repeated  until  the  full  load 
is  applied.  During  the  first  pass  through  this 
section,  subroutine  PRNT  is  called  which  calculates 
and  prints  the  linear  analysis. 

Checking  auxiliary  tapes  and  storing  information  to 
be  used  as  input  to  plot  routine  SAN1PLT 

End  of  main  program 

Subroutine  EQSOL- -Solves  the  equations  for  the  dis¬ 
placements 

Subroutine  PACK- -Assembles  the  master  stiffness 
matrix 

Subroutine  BAND- -Determines  indexing  matrices 
required  for  banding  the  stiffness  matrix 

Subroutine  MULT- -Obtains  the  product  of  two  matrices 

Subroutine  NfTRAN- -Multiplies  one  matrix  by  the 
transpose  of  another 

Subroutine  CMPAT- - Generates  the  compatibility  matrix 
for  each  element 

Subroutine  ELEMK- -Generates  the  elemental  stiffness 
matrices 

Subroutine  PENT- -Calculates ,  prints,  and  sets  up  data 
for  plot  routine  SAN1PLT  for  the  linear  solutions 


SANOS 

PROGRAM  SANOS(IM»UT,3JTPUT,  f  APE7,  FILMPD 

DIMENSION  T  « » 3 0 >  , XIXI 8  30) , ARE  A  (830 > ,PM (810)  , JEM (8 30)  ,KEH(833) , 

2  XJ(830)  ,YJ( 130), ZJ (8301 , XK(830> ,YK(813> ,ZK(830>, 

3  IX (It >5)  ,0(1155)  ,  JDin(U55l  ,JBR0(1155>,  NSC  1155)  , 

4  NSP(ti55) ,3( 1153) ,X0 <1155) v  AXI AL  ( 830)  ,PUL T(8U ) , 

5  SOM  (1563),SOM2(1660)  ,  SUM  ( 16'60)  , 

6  S<65  • 11551 , 

7  TKUO.tO)  ,3(10,tQ),P(10>  , SOI 10, 10)  ,ESUO,tO>, 

8  OO(LO)  ,P'5E'4T«3a>  ,C^SC21>  ,SRN(21>  ,NP  (5), 

9  Y (2311 
*,MMTYa(833> 

COMMON  7 2/  X(2U),  2(231),  ?LOTAX(220), 
lDLOT  NJ(220),PLOTMK(221)*STRS0J(220>  ,STRSIJ(220>  ,STRSOK(220) , 
2STRSIK  (223 )  ,MAN3LE 
COMMON/3/N8AYL , NBA YW,  NPR03 
COMMON/4/°LOTN 
COMM3N/5/OA, JPROB, HUM PR 
REAL  .SUM 

232  FORMAT (313,111) 

301  FORMA  T (8110) 

362  FORMAT (1H0) 

420  FORMAT (20X14H3R33LEM- NUMBER  15/) 

421  F0RMAT(16X,41H  NONLINEAR  ANALYSIS  . OF  ORTHOTROPIC  SHELLS  ) 

422  FORMAT  ( 14X , 45H  MATERIAL  NONLINEARITY  *»  RICHARD  FORMULATION  ) 


424  FORMAT (49H  NUMBER  OF  RUNGE-KUTTA  INTERVALS  ******  I 

425  FORMAT (49H  NUMBER  OF  MEMBERS  ******************** 

426  FORMAT ( 49H  NUMBER  OF  LOADS  **********************  NU 

422  FORMAT (49H  NUMBER  OF  SUPPORTS  •**•»*******»••*•**  NSI 

428  FORMAT (49M  BAM3  WIDTH  *•****••***•••*****•*•*•***  NO 

430  FORMAT (49H  CURVE  FITTING  PARAMETER  ************** 

443  FORMAT (53H  ******************SHELL  GEOMETRY*********** 

441  FORMAT (61H0MN  JM  <M  AREA  ULT. FORCE  SE 

IT. MOM,  8X.6HLSN3TM,  8X,5HMTYPE  ,/) 

442  FORMAT (3I4,5F12.3, 18) 

450  FORMAT (18X24H  RJNGE-<UTTA  INTERVAL  *  13) 

451  FORMAT (26X8H  MOMENTS) 

.452  FORMAT (123H  MEM3ER  NO.  MOM.  AT  J  END  MOM.  AT  K  END 
t  AXIAL  FORCE  ULT.  MOM.  MJ/PM  MK/PM 

453  F0RMAT(I7,4E12.5,4F11.3) 

454  FORMAT (27X14H  DISPLACEMENTS) 

455  F0RMAT(12JH  NODE  DELTA  X  DELTA  Y 

1  THETA  X  THETA  Y  THETA  Z 

456  F0RMAT(I6,1E17.5,5E13.5> 

459  FORMAT (19* 1P2E17 ,5) 

460  FORMAT ( 1 0 X29H  THESE  ARE  THE  ELASTIC  FORCES//) 

461  FORMAT (42H  THESE  ARE  THE  ELASTIC  NODAL  DISPLACEMENTS/) 


NRK 

NM 

NLOAO 

NSUPS 

NDANO 

R 


) 

13) 
Id) 
13) 
Id) 
Id) 
F9.2) 
* **  ) 


SECMO 


SHEAR 
ULT.  FORCE) 


DELTA 


494  PQRMA T ( 49H  THICKNESS  OF  SMELL  (BENDING)' 

495  FORMAT (49H  THICKNESS  OF  SHELL  (AXIAL)** 


THICKM  =F9.2) 


497  FORMAT (49H  M03ULJS  0?  ELASTICITY  OF  SHELL 

502  FORMAT (49M  NUMBER  OF  NOOES  ************ 

503  FORMAT (49H  NUMBER  OF  EQUATIONS  ******** 
534  -0RMAH49H  NUMBER  OF  )AYS  LONG  ******** 
505  FORMAT (49H  NUMBER  OF  BAYS  WIDE  ******** 
507  FORMAT (49H  LENGTH  OF  DIAGONAL  ELEMENTS 


THICKA 

E 

NODES 
N0D5 
NBAYL 
NBAYW 
01  AG 


-Fit  2) 
=-9.2) 
=  Id) 
=  Id) 
=  Id) 
=  Id) 
*F9. 2 ) 
•****) 


510  FORMAT (53H3*************THE  NUMBER  OF  EACH  SUP»3RT*************) 

511  FORMAT (1215) 

512  FORMAT (53HO**************LOAD  NUMBER  AND  LOAD**** *•*•******»*•» ) 
514  F0RMAT(I4,F3.3,I16,F8.3,I16,F8,3) 

516  FORMAT (5F12.5) 

517  FORMAT  ( 5 3H0 *** 3E ^(JENT A GE  OF  LOAD  ON  EACH  RUNGE-KJ  TTA  INTERVAL**) 

518  FORMAT (48H  PERCENT  0-  TOTAL  LOAD  TAKEN  IN  THIS  INTERVAL  =  fr.O) 
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^ mw «u*^u WrtWW aasmu^i^hmu  'i J'l^J  ■■pwwsf?^s.yi,^ysiW' u wijy>* fk^yj^>y 1 ^  -v-**-  ~<L-  * *>& 


F9.0I 


f 


a  2 

519 

63 

550 

54 

551 

55 

552 

56 

630 

67 

631 

58 

520 

59 

521 

73 

623 

71 

900 

72 

901 

73 

13304 

74 

10005 

75 

13336 

76 

1000  7 

77 

10008 

78 

10009 

79 

23333 

90 

2000  1 

81 

23332 

92 

20003 

93 

20009 

94 

85 

C 

86 

97 

1 

88 

C 

99 

C 

90 

C 

91 

C 

92 

C 

93 

r 

■/ 

94 

c 

95 

c 

96 

97 

98 

99 

130 

131 

132 

4 

133 

134 

5 

105 

6 

13  6 

10  7 

138 

13  9 

110 

ill 

112 

113 

114 

115 

116 

117 

118 

981 

119 

120 

121 

12  2 

123 

124 

892 

125 

FORMAT (48HOTOTAL  ®£RDSNT  OF  LOAO  INCLUDED  AT  THIS  TIHE  = 

FORMAT < 1  HI » 

FORMAT (15X,40H  ELASTIC  ANALYSIS  **  SMITH  PLASTIC  MODEL  > 

FORMAT (12X,48H  GEOMETRIC  NONLINEARITY  **  INCREMENTAL  TECHNIQUE  ) 


FORMAT (1P5E20. 8) 

FORMAT (1P9E13.3) 

FORMAT ( 49H  HEIGHT  OF  STRUCTURE  *******************  HEIGHT  =-9.2> 
FORMAT (49H  SPAN  OF  STRUCTJRE  ***********************  SPAN  =?9.2> 
FORMAT ( 49M  CENTRAL  ANGLE— 0 EGREES* *************** **•*  CAD  =-9.2) 
FORMAT (8F10.0) 

FORMAT  < II 10 • IF13 

FORMAT ( 49M  THICKNESS  OF  SHELL  (BENDING-WEAK  AXIS)  THICKW  =F9.2) 
FORMAT (49H  HALF  9AN0  WIDTH  PLUS  DIAGONAL  ********  NBSYM  =  la) 

FORMAT  (49H  LENGTH  OF  SHELL  **********************  SENGTH  =.-9.2) 

FORMAT ( 49H  LENGTH  OF  LONGITUOINAL  MEMBERS  *******  SIDE  =?9.2) 

FORMAT (49H  RADIUS  OF  SHELL  **********************  RADII  =-9.2) 

FORMAT (49H  ULTIMATE  MOMENT  (IN-KIPS/IN)  *********  ULTM  =-9.2) 
FORMAT (49H  LARGER  NODE  NUMBER  ON  END  OF  MEMBER  **  NUMBE  =  la) 
FORMAT (49H  SMALLER  NODE  NJM9ER  ON  ENO  OF  MEMBER  *  NUMSE  =  IS) 
FORMAT (49H  NUMBER  OF  NODES  WITH  SUPPORTS  ********  NSUPN  =  la) 

FORMAT (49H  METHOD  OF  NUMBERING  NODES  ************  NCOOE  =  la) 

FORMAT (49H  ULT.  AXIAL  STRESS  (KIPS/SQ. IN/INCH)  **  ULTSTR  =-9.2) 
READ  131, NUMPR 


NUMPR  **  NUMBER  OF  PROBLEMS  TO  BE  SOLVED  IN  THIS  COMPUTER  RJN 
DO  999  JPROB=l, NUMPR 

READ  301,  NPR03, NCOOE  ,NBAVL,NBAYW,NLDAO , NSUPN, NRK, MANGLE 

N°R09  **  PROBLEM  IDENTIFICATION  NUMBER 

NUMBE  **  NOOE  NJM3ER  ON  BIS  END  OF  MEMBERtFOR  BAND  WIDTH  CALC) 
NUMSE  **  NOOE  NUMBER  ON  SMALL  END  OF  MEMBER (FOR  BAND  WIDTH  CAlC) 
N3AYL  **  NUMBER  OF  9AYS  ALONG  SHELL  AXIS 
NBAYW  **  NUMBER  OF  BAYS  along  ARCH 
NLOAO  **  NUMBER  OF  LOAOS  ON  SHELL 
NSU°M  **  NUMBER  OF  NODES  WHICH  HAVE  SUPPORTS 
SENGTH  **  LENGTH  OF  SHELL 
READ  900, SPAN, HEIGHT, SENGTH 
RE AO  903,THICKM,THIC<W,ULTM,THICKA,ULTSTR,E,R 
REAO930,C0UT,CIN,RATMO0 
READ  900, (PCENT(I) ,1*1, NRK) 

MQ=N3AYW*N3AYL 
IF  (N3AYL-NBAYW)4,5,5 
NUMBE  =  NBAYL+3  - 
GO  TO  6 

NUMBE  =  NBAYW  +  3 
NUMSE  =  1 

N9 ANJ= 10*( NUM3E” NUMSE) +9 
N3SYM= (NBANO+1) /? 

HH  =  N3AYL*(4*N3AYW*1)  MBAYW 
NOOES= (NBAYW+1) *(NBAYL+1) 

N005  =  5*NODES 
3IQE=SENGTH/N3AYl 
PI=3. 14159265 

BET  A  =  *!/  2. -2.*  ATANr  (HEIGHT/ (SPAN/2.)) 
RA0II=HEIGHT*39AN*TANF(BETA)/2, 

CA=PI-2.*9ETA 
CAO=CA*180 ./PI 
00  981  I=1,N005 
NSP(I) =0 
NSUP3=0 

DO  893  1=1 ,NSJ°N 

REAO  331, NOOE, (N-*  ( J) , J=1 ,5) 

00  891  K  =  l ,  5 
IF(N»(K))  892,891,892 
NCT=5*N00E-5»NP(K> 

NSP(NOT) =1 
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126 

NSUPS=NSU»Sfl 

127 

NS  (NS JPS) =NGT 

128 

891  CONTINUE 

129 

890  CONTINUE 

130 

IRK=1 

131 

NPRNTsO 

132 

LSUM=3 

133 

RN=(1.0+R)/R 

154 

00  342  J=l,N035 

135 

XQ ( J) =0 . 0 

136 

9(  J) *0.0 

137 

342  CONTINUE 

138 

III=NN*2 

139 

DO  331  1*1,111 

t40 

SUMt I) *0 . 0 

141 

SOM (I) *0.0 

142 

331  33M2(D  =0.0 

143 

90  235  MN=1,NN 

144 

235  AXIA. (MN>=3 .3 

145 

00  236  J=1 , NLOAO 

146 

236  READ  901,1,0(1) 

147 

°RINT  550 

148 

PRINT  423,NPRD9 

149 

PRINT  421 

150 

°RINT  422 

151 

PRINT  551  • 

152 

PRINT  552 

153 

“RINT  424, NRK 

154 

PRINT  425, NM 

155 

PRINT  502, NODES 

156 

PRINT  503.NOD5 

157 

PRINT  5 3 4, NBA ft 

158 

PRINT  505, N8AYW 

159 

PRINT  426, NLOAO 

153 

PRINT  427, NSUPS 

151 

PRINT  20002, NSUPN 

152 

PRINT  428, N8AN9 

153 

PRINT  23033 »NUM3£ 

154 

PRINT  20001, NJM5E 

155 

PRINT  20033, NCOOE 

166 

PRINT  13035, N8SYM  > 

16  7 

PRINT  430, R 

158 

PRINT  495, THI3KA 

.159 

PRINT  494,  THI C<M 

t73 

PRINT  10004,  TMI3KH 

171 

PRINT  10009, ULTN 

172 

PRINT  23039, ULTSTR 

173 

PRINT  623, HEIGHT 

174 

PRINT  621, SPAN 

175 

PRINT  623, CAD 

176 

PRINT  10036, SENGTH 

177 

PRINT  10007, SIOE 

178 

PRINT  10008, RADII 

179 

PRINT  497, E 

180 

PRINT  362 

181 

“RINT  550 

19  2 

PRINT  440 

133 

C 

•  GENERATE  SHELL  3D0R0I 

184 

C 

ORIGIN  ALONG  3ENTERLI 

19  5 

N3AYUN3AYW+1 

186 

AA=0. 0 

137 

PI=3. 14159265 

198 

C 

MANGLE  =  1  ****  ONE-M 

18  9 

c 

MANGLE  =  2  ****  FULL 

3H  (93  DEGREES) 
(180  DEGREES) 
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19C  SO  TO  (863*881) *  MANGLE 

191  881  SEG  =  N3AYW 

192  GO  TO  882 

193  880  SEG  =  N3AVW  *  2 

194  882  CONTINUE 

195  00  9  1=1 »MBAYl 

196  CRS{I>=COSF{AA*3A/SE3+BETA)*RAOII 

19  7  SRN(D  =SINF(AA*CA/SE3*BETA>  *RADII 

198  9  AAsAAH.O 

199  JCT=0 

>3  a  33=0.0 

>31  NBAY2=N3AYLM  «N3A Y3*NBA Yl»N8 A Y2 

>02  C  NCOD£=0 ““NODES  NUMBERED  IN  DIRECTION  OF  ARCH  AXIS 

>33  C  NS00£=1--N00E5  NUMBERED  ALONG  ARCH 

>34  IF(NC00E)16,17,16 

>35  16  CONTINUE 

236  00  11  1  =  1*  N3AY2 

>37  00  18  J=l*  NBAY1 

208  JCT=J0T*1 

>39  X (JCT) =“CRS( J)  . 

>10  Y(JCT) sOB’SIOE 

211  Z(JCT) =SRN(J)-(RAOII*M£IGHT> 

212  18  CONTINUE 

213  10  83=99H.O 

214  GO  TO  112 

215  1 7  CONTINUE 

216  00  110  1=1 i N8AY1 

217  93=0.0 

218  00  111  J  =  1 »NBAY2 

>19  JCT=J3TU 

220  X  < JCT) =-CRS (I) 

221  Y(JCT)=BB*SIO; 

222  Z(JCT)=SRN(I)-(RAOII“H EIGHT) 

223  111  9B=B3fl.0 

224  .110  CONTINUE 

225  112  CONTINUE 

226  “RINT520 

22 7  523  FORMAT (51H  NOOE  X  COORDINATE  Y  COORDINATE  Z  COORDINATE/) 

228  00  11  N=l, NODES 

229  11  PRINT  521,N,X<N)  ,  YIN) ,Z(N) 

>10  521  FORMAT (I5»F11«3|2F16.T) 

231  C  OOTAINING  MEM3ER  PRO’ERTISS  AND  GEOMETRY  RELATIONSHIPS 

212  SaT2=S0RT- (2.0) 

253  0  A44=  .ENGTH  OF  ELEMENTS  ALONG  ARCH 

234  C  OBB  =  LENGTH  Or  ELEMENTS  ALONG  SHELL  AXIS 

>35  AAA=2.9*RA0II*SINF(CA/(2.0*SEG)> 

236  333=  SENGTH/N3AYL 

237  00  13  IJK=1*NM 

238  READ  2 02 , MN, JM* KM , MTYPE 

239  XJ(MN)=X(JM) 

>40  YJ(MN)=Y(JM> 

241  ZJ(MN)=Z(JM> 

242  XK(MN)=X(KM) 

243  YK (MN) =Y (KM) 

244  ZK (MN) =Z (KM) 

245  JEM  (HN )  =  J.N 

246  KEM(NN) =KN 

24/  C  MTYPE  =  1  ********************  INTERIOR  MEMBER  ALONG  ARCH 

248  C  MTYPE  =  2  ********************  EXTERIOR  MEMBER  ALONG  ARCH 

249  C  HTTP-  =  3  ********************  INTERIOR  HORIZONTAL  MEMBER 

250  C  MTYPE  =  4  ********************  EXTERIOR  HORIZONTAL  MEMBER 

251  C  MTYPE  =  5  *********************diag0NAL  MEMBER  . 

252  GO  TO  (250,260»2/0»280,290)MTYPE 

253  250  SECMO=  933*THICKH**3/12. 
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254 

AR£A{NN)=9BB*THI0KA 

>55 

t»JLT<MN>  =  BBB*HICKA*JLTSTR 

256 

°M(MN>=BB9*ULTM 

252 

GO  TO  12 

258 

260 

SECM0=.5*3B8*THICKM**3/12. 

259 

AREA  <  'IN)  =. 5*B33*THIC<A 

250 

PULT(NN) =«  5*B3B*THIC<A*ULTSTR 

261 

PM<MN)=.5*BBB*ULTM 

252 

GO  TO  12 

253 

270 

SECM0=AAA*THICKW*»3/t*. 

254 

AREAINN)  =AAA*THI3<A 

255 

OULT ' IN) =AAA»THI2KA*JLTSTR 

256 

PMIHN)  =AAA*ULTM 

257 

GO  TO  12 

258 

280 

SECM0=.5*AAA*THICKW**3/12. 

259 

AREA ( HN) =.5*AAA*THIC<A 

270 

PULT<NN>=.5*AAA*THIC<A*ULTSTR 

271 

PM<MN>=.5*AAA*ULTM 

272 

GO  TO  12 

273 

290 

SECMO=. 01 

274 

AREA ( NN J  =• 01 

275 

PULT<NN>=10. 

276 

PM<MN»  =10. 

277 

12 

CONTINUE 

278 

C  12 

PRINT  442,MN,JM,<M,AREA(MN>  ,PULTCMN>  ,SECMO,PHINN> 

279 

MMTY*<MN>  ®MTY*E 

280 

13 

XI X ( IV) =SECMO 

281 

PRINT  550 

282 

PRINT  510 

283 

°RINT  511, <NS(I> ,I*1,.NSUP3> 

284 

PRINT  512 

285 

PRINT  516,  (9(11 , I*1»N095) 

286 

►*RINT  517 

287 

PRINT  516, <PCENT(I), 1=1, NRK> 

288 

00  237  1=1 ,N005 

289 

237 

Q(I)=0(I> 

290 

C 

DETERMINE  INOEXING  MATRICES  REQUIRED  FOR  3 AMOINS 

291 

CALL  3ANO  (NODES, NBA  NO, IX, JBLD, JBRDJ 

292 

00  413  MMNMM=1,NRK 

293 

RUNGs&.O 

294 

FACT=0 . 0 

295 

200 

CONTINUE 

296 

I  J=3 

297 

00  333  1  =  1  ,NB5YM 

298 

00  330  J=1,N0D5 

299 

330 

S(I,J)=0.0 

300 

00  549  J=1 , N0D5 

301 

549 

0(J)=Q(J)*PCENT(NMM!H> 

302 

00  14  MN=1 ,NM 

303 

C 

GENERATE  COMPATIBILITY  MATRIX,  B 

304 

CALL  CMPAT(XJ»YJ,ZJ*<<»YK,ZK»T»B»NN) 

305 

C 

GENERATE  ELEMENT  STIFFNESS  MATRIX,  ES 

30  6 

CALL  ELEMX<MN,XIX,T,S,SOM,PM,R,RN,ES  ,PULT,AXIAL( 

30  7 

C 

MULTIPLY  COMPATIBILITY  MATRIX  BY  ES 

30  8 

CA’LL  NULT(3,10,3,  S9,ES,B> 

30  9 

C 

DEVELOP  STRUCTURE  ORIENTED  STIFFNESS  MATRIX,  ZK 

310 

CALL  NTRAN(10,10,3,T<,B,S3) 

311 

15  JMM=JEM  (N't) 

312 

<MM=<EM(MN) 

313 

C 

ASSEN3LE  THE  MASTER  STIFFNESS  MATRIX,  S 

314 

CALL  3ACK  (JMM,KNM,IX,TK,S> 

315 

14  CONTINUE 

316 

C 

SOLVE  FOR  THE  DISPLACEMENTS 

317 

CALL  EQSOLCNOOES.S, J, IX, JBLO, J9R0,NSP) 
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31 «  IF (N’RNT )  678, 680,673 

319  683  CONTINUE 

320  888  PRINT  550 

3>1  PRINT  441 

322  00  339  MN=1,NH 

323  PRINT  442,MN,JEM(NN>,  K£N(MN) , 

324  1  AREA(MN) ;PULTCHN>,XIX(MN),PN<NN),T(MN) 

325  2  » NMTY3 (MN) 

326  889  CONTINUE 

327  PLOTN=0. 

328  R<=1.0/PCENT(t> 

329  CALL  PRNT(NN,J;M,K£N,T,XJ,YJ,ZJ,XK,YK,ZK,XIX,E,RK,N005,0, 

330  1  N°ROB,  AREA, MHTYP,OIN,COUT,RATNOD,Q, CONTROL) 

331  NPRNT  =  NPRNT«-1 

332  678  CONTINUE  '  ’ 

333  PLOTN=l. 

334  IF(NR<)  683,  998,  683 

335  683  CONTINUE 

336  C  INTEGRATING  THE  INTERNAL  MOMENTS 

337  IJ=0 

338  00  42  MN=1,NM 

339  J5=JEN(MN>*5 

340  00 (1) =0  C  J5-4) 

3,1  90(2)=0JJ5-3> 

342  00 (3) =9 ( J5-2 ) 

343  00(4)=0(J5-1) 

3,4  00  <  5) =0 ( K5) 

345  <5=KEN(MN) *5 

346  00(6)  =  0U5-4) 

34  7  00 ( 7) =0 ( K5-3) 

348  )0{8)=9(K5-2) 

349  OOt  9) =0  (K5-1) 

350  00(10)=0(K5) 

351  CALL  CMPAMXJ,  YJ , Z J, XX, YK, ZK,T ,9, MN) 

352  XJ(NN)=XJ  (HN)  03(1)  /  RUNG 

353  Y J( NN) =  VJ (NN) frOO ( 2) 7  RUNG 

354  ZJ(NN)=ZJ(MN)f)0(3)/RUNG 

355  XK(MN)=XK(MN)+00(6)/RUNG 

356  YK( NN) = YK(MN)  *03(7)7  RUNG 

357  ZK(NN)=ZK (MN) *03(8) 7 RUNG 

358  CALL  ELEMK(MN,XIX,T;2,S0M*PM,R,RN,ES  ,PULT, AXIAL, AREA) 

359  CALL  NULT(3,10»3,S9»-<3,8) 

350  C  P  =  (S3)  *  0 

351  CALL  MULT(3,1,13,P,S3,00) 

362  AXIAL  ( MN)‘= AXIAL  (HN)  fP(l)  /  RUNG 

333  00  622  1=2,3 

334  IJ=IJH 

355  5UM(IJ)=SJH(IJ)«-°(I)7RUNG 

336  S0N(IJ)=S0M2<IJ)fP(I)*FACT 

357  622  CONTINUE 

358  42  CONTINUE 

359  C  BEGIN  RUNC-E-KJTTA  INTEGRATION 

370  GO  TO  (700, 701, 702, 733), IRK 

371  700  00  8)3  J=1,N035 

372  800  XO (J) =  X0 (J)*0( J) 76.0 

373  FACT=. 50 

374  RUNG=3 . 0 

375  GO  TO  100 

376  701  00  8)1  J=1 »N035 

377  801  XO (J) =X0 (J) +0(J) 73,0 

378  FACT=.50 

379  RUNG  =  3 . 0 

380  GO  TO  100 

381  702  00  802  J=1,N005 
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382 

832 

X0(J)*XQ(J)*D(J)/3.0 

383 

FACT =1 • 0 

334 

RUNG*6.3 

385 

GO  TO  133 

386 

703 

90  833  J=l,NO<J5 

387 

803 

X0  ( J) *X0 { J) ♦D( J) 76.0 

388 

R'JNG=6.3 

389 

100 

CONTINUE 

390 

IF (IR'<-4)  482,490.493 

391 

493 

CONTINUE 

392 

PRINT  550 

393 

PRINT  420, NPR09 

394 

PRINT  362 

395 

PRINT  450 , MMMNM 

396 

TSUM=*CENT (MMMHN) *  109 • 0 

397 

LSUM*LSUM»-TSUN 

398 

°RINT  519,ISUM 

399 

PRINT  518, TSUN 

430 

PRINT  362 

401 

PRINT  454 

43  2 

PRINT  455 

43  3 

JRO=3 

404 

00  671  JR*5,N005» 5 

43  5 

JRO= JRO+ 1 

436 

671 

PRINT  456, JRO.X3 ( JR-4) ,X3(JR-3),X0 (JR-2) ,X3 (JR-l) ,X3 (JR) 

407 

PRINT  362 

43  8 

PRINT  362 

40  9 

PRINT  451 

413 

PRINT  452 

411 

491 

IBD=3 

412 

NMT2=NM*2 

413 

00  652  I A=1,NNT2» 2 

414 

I-30=I39*1 

415 

ENOM J=SUM (IA) 

416 

EMON<=SUM(IA+l) 

417 

SHEAR®  (EM0MK-EM0NJ)/TCI8D) 

418 

oPoJ®EH0MJ/PM(I99) 

419 

PPP<=EM0NK/PN(I30) 

420 

PAX  =  AXIAL(  130)  /PJLTdlD) 

421 

552 

“PINT  453, 130, EMONJjEMONK, SHEAR, AXIAL (I3D) ,?M (133 ) ,??PJ,?P?<,?Ay 

422 

482 

CONTINUE 

423 

IRK=IRK+ 1 

424 

IF (IR<-5)  203,99,  99 

42  5 

99 

IRK=1 

42  6 

00632  1*1, NMT2 

42  7 

SOM ( I ) =SUN ( I ) 

428 

632 

S0M2(I)*SJM(I> 

429 

399 

CONTINUE 

433 

433 

CONTINUE 

431 

PRINT  900, CONTROL 

432 

PRINT  30000, (X0(I) ,1*1, N095) 

433 

33030 

FORMAT (5(E13.5,5X) ) 

434 

IF  (CONTROL  ,E(J»  0» )  GO  TO  998 

435 

30001 

IF (UNIT, 7) 30001,850,30002,30003 

436 

30002 

PRINT  835 

437 

30333 

PRINT  804 

438 

804 

FORMAT (’PARITY*) 

439 

805 

FORMAT (*EOF*> 

443 

853 

3UFF* ROUT (7,1)  (X3(l),X0(N0O5)) 

441 

21000 

IF  (UNIT,  7)  21000, 8  51, 851, 851 

442 

851 

CONTINUE 

443 

ENO  FILE  7 

444 

998 

CONTINUE 

445 

PRINT  550 
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1*46  999  CONTINUE 

447  END 

448  SUBROUTINE  EQ30L  (  M2, SS  ,P  ,  IX,  JL,  JR»  NSP) 

44  9  DIMENSION  SS (65,1155) , P (1155)  ,1X11155)  ,JL  ( 1155) ,  JR(1 155) ,  NSR(  1155) 

450  M90ES=M2*5 

451  00  423  11=1, NODES 

452  IF  (NS3(II> )  423,425,420 

453  425  ',5  =  IX(II) 

454  =  1 . 0/SS (K3, II ) 

455  JRT=  JR(II) 

456  00  434  JJ  =  II,JRT 

457  434  SS  (K3, JJ)=SS<<8, JJ)*4A 

459  ?(II> =?<II)*AA 

459  00  424  LL=1,JRT 

460  LB  =  IX(LL) 

461  IF(L.-II)  424,424,425 

462  428  IF(NSP(LL) >424,421,424 

463  421  93  =  5SU3,LL> /4A 

464  IF(BB)  423,424,423 

465  423  00  422  JJ=LL,JRT 

466  422  SS(L3,JJ>=SS(LB, JJ)-3S <<9, JJ) *89 

467  P(LU =P(L.)-P(II)*8B 

468  424  CONTINUE 

459  420  CONTINUE 

470  00  452  II=2,N00ES 

471  JJ=N00ES*1-II 

472  IF  <NSP ( J J) )  451,451,452 

473  451  J8=IX ( JJ) 

474  '  <K=II-1 

475  LL=J.(II) 

476  00  450  MM=LL,KK 

477  NN  =  NODES  ♦  l  -  MM 

478  450  P(JJ)=P(JJ)-S5(J3,NN)*P(NN) 

479  452  CONTINUE 

480  RETURN 

491  ENO 

432  SUBROUTINE  ?A3K(  JM,M,  IX,ZK,SS) 

493  DIMENSION  IX (1155 ) , Z< ( 10 , 10) , SS (65, 1155) 

494  IP3=JN*5 

495  I03=KN*5 

486  00  383  IR0W»1,5 

487  JR0W=5-IR0W 

498  IIP3=IP3-JR0W 

499  IIQ3=IQ3-JR0W 

430  <R0H=IX(II?3) 

431  LROW=IX (IIQ3) 

432  00  351  IC0L=1, 5 

493  JCOL  =  5-1  COL 

494  KC0L=IP3-JC0L 

435  LC0L=IQ3- JCOL 

436  IF(I=>3-IQ3>  315,315,314 

497  315  SStKROH,LCOL)=SS< KROd , LCOL) +  ZK (IRON, I COL* 5) 

498  314  IF  < 1*3-1 Q3 )  319,307,307 

499  307  SS (L ROW, XCOL) =SS ( L RON , KCOL) *ZX  (IROWN-5,  IDOL) 

500  319  CONTINUE 

531  381  CONTINUE 

502  00  332  IC3L=IR0W,5 

503  JCOL=5-ICOL 

534  KC0L=IP3-JC0L 

535  LCOL  =IQ3 -JCOL 

50  6  SS ( KROW, KCOL) =SS ( KROW , KCOL)  +Z K (IRON , ICOL ) 

50  7  SS (L ROM, LCOL) =SS ( L ROW*  LCOL) ♦ZK(IR0W+5,IC0L+5> 

5)8  382  CONTINUE 

509  380  CONTINUE 
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510 

511 

512 

513 

514 

515 

516 

51 7 

518 

519 

520 

521 

522 

523 

524 

525 

526 
52/ 

528 

529 

530 

531 

532 

533 

534 

53  5 
536 
53/ 

538 

539 

54  0 

541 

542 

543 

544 

545 

546 
54/ 

548 

549 

550 

551 

55  2 

553 

554 

555 

556 
55/ 
556 

559 

560 

55 1 

552 

553 
564 
55  5 
556 
55/ 

558 

559 
5/0 
5/1 
5/2 
5/3 


RETURN 

END 

S'JBROJTINE  BANQ(X1,N3AND,IX,JL , JR) 

DIMENSION  IX(1155> ,JLC1155> ,JR(1155> 

NOOE5*5*Mt 

NBSYN  *  (NBANO  ♦  l)/2 
KOUNT  =  1 
DO  3/3  11*1, NODES 
IF(II-<OUMT  *  NBSYN)  3/1,3/1,372 
3/1  IX  (II)  *11  -NBSYN  MKOUNT  -1) 

GO  TO  3/0 

3/2  IX<II) *II-NBSYN  *  KOJNT 
KOUNT  =  KOUNT  *  1 
3/0  30NTINUE 

JRT  *(NBAN0  -t>/2 
JS  =N3DES  -  JRT 
IS  * JRT  f  1 
00  340  11*1, NODES 
IF(II-IS)  350,350,352 

350  JLT  *  1 

JRT  =  JRT  ♦  1 
GO  TO  356 

352  IF  (II- JS)  351,353,353 

351  JLT  =  JLT  ♦  1 
JRT  *  JRT  ♦  1 
GO  TO  356 

353  JLT  *  JLT  ♦  1 
JRT  *  NODES 

356  JL  (II)  *  JLT 
JR(II)  *  JRT 
340  CONTINUE 
RETURN 
END 

SUBRD JTINE  MULT(L,M,N,RESLT, AA,8B) 

DIMENSION  RESLTdO  ,13)  ,AA(1Q, 10)  ,33(10,10) 

C  RESLT  =  AA  *  83  AA  IS  L  X  N  BB  IS  N  X  N 

DO  133  1*1, L 
DO  133  J*1,M 
RESLT(I,.J)  *0.0 
00  103  K*t,N 

10  0  RESLT  (I, J)*RE3LT(I,J>* AA  (I,  <> '33  «,  J) 

RETURN 

END 

S'JBROJTINE  MTRAN(L,M,N,RESLT,  AA,B3) 

DIMENSION  RESlT(10,13) , AA(10,10) ,BB(10,10) 

C  RESLT* ( A A) T  ♦  B3  AA  IS  NXL  8B  IS  NXM 

DO  133  1*1, L 
DO  100  J*1,M 
RESLT (I, J) *0,0 
00  133  (<=1  »N 

100  RESLT ( I , J) sRESLT (I,J)tAA(K,I)*BB(K,J) 

RETURN 

END 

SUBRO JTINE  CNF  AT (XJ,YJ,2J,XK,YK,ZK»T , B, MN) 

DIMENSION  B(10,10) ,XJ(830),YJ(R30) ,Z J< 830) , XK (830 ) , 
2  ■ V < ( Q 3 0 )  , ZK(S30) , T (830) 

00  13  1*1,3 
DO  10  J*l,10 
10  3(1, J) *0.0 

XX*XK(NN) -XJ(MN) 

YY*Y<(MN) -YJ(NN) 

ZZ*ZK(MN) -ZJ(NN) 

T(MN)  =SQRTF(XX'XX'YY»YY*ZZ'ZZ) 

XY  =  SQRTF(XX*XX«-YY*YY> 


65 


574 

CA  =  XX/XY 

575 

SA  =  YY/XY 

576 

33  =  XY/T1MN) 

577 

S3  =  ZZ/T  (HN) 

578 

9(1,0=  CA*CB 

379 

3(1, 2>  =  SA*C3 

530 

3(1,3)  =  38 

581 

3(1,S)=-B(1,1) 

532 

3(1., 7)  =-3(1,2) 

533 

0(1, 8) =-0(1,3) 

584 

3(2,1)  =  CA*SB/T(MN> 

33  5 

•3(2,2)  =  SA*S3/T(MN) 

58  6 

3(2,3)  =  -C9/T(MN) 

387 

3(2,4)  =  -SA 

538 

3(2,5)  =  CA 

539 

3(2,3) =-0(2,1) 

590 

3(2,7) =-9(2,2) 

591 

3  (2, 8)  =-B(  2 , 3) 

592 

0(3,1) =-0(2,1) 

593 

9(3, 2) =-9(2, 2) 

394 

3(3,3) =-8(2,31 

595 

3(3,6) =-8(2,6) 

596 

0(3,7) =-9(2,7) 

397 

9 ( 3, 8) =-B ( 2, 8) 

598 

3(3,91  =  -3(2,4) 

599 

9(3,10)  =  -0(2,5) 

30  0 

RETURN 

301 

END 

30  2 

SUBROUTINE  ELEMK < NN, SECMO, T ,E , SUM, P«,R,RN, ES, PJLT , AX IAL, AREA > 

33  3 

DIMENSION  SECMO (3 30) ,  T  (8  30)  ,3MC830>  , SUM (1660)  *ES(  10, 10) 

304 

2  ,  AREA  (3  30),  PULK830)  ,  AXIAL  (830) 

305 

JZ=MN*2-1 

606 

AA=A3SF(AXIAL(NN)/PUlKMN)) 

33  7 

IF(AA-l.O)  50,60,60 

308 

60 

AA=5.3*AA 

53  9 

50 

CONTINUE 

310 

PPIJ=AOSF(SUM(JZ)  /PMCMN)  ) 

311 

PPJI  =  A9SF(SUM(JZ  +  1)/:*'1(MN) ) 

312 

IF  (P*I J-1.0)  10,10,20 

313 

20 

PPIJ=3. 99999 

314 

10 

IF  (P*JI- 1, 0 )  30,30,40 

515 

40 

PPJI=0. 99999 

516 

30 

ES (1,1)= ( AREA (MM)  *E/l  * MN> ) / (1 ,0  +  AA**R) **RM 

5t7 

ES (1, 2  >  =  0. 0 

618 

ES (1, 3 ) =0 . 0 

319 

ES (2, 1) =0.0 

52  0 

ES (3,1) =0. 0 

521 

AI=(l.O-P®IJ**R) **R.N 

322 

AJ=(t.O-PPJI**R) **RN 

523 

QROO  = (5. 0*E*SECMO(MN) /T (MN) ) / (4. O-AI* A J) 

324 

ES(2,2)=PROO*2.-3*AI 

525 

ES  (2,5)  =  -3R0D'1  AI*  A  J 

526 

ES(3,2)=ES(2,3) 

527 

ES(3,3>=FROD*2.0*AJ 

328 

RETURN 

529 

ENO 

530 

SUBRD JTINE  PRNT ( NM, JEN ,KEM, T, X J, YU ,Z J, XK, YK,ZK, XI X,E» RK,N0D5, D, 

331 

1NPR03, AREA, MTYP-,C IN,  COUTjRATMOD, Q, CONTROL) 

332 

DIMENSION  Q ( 1) 

333 

C 

ONE  INPUT  CARO  IS  REQUIRED  IN  SUBROUTINE  PRNT- 

334 

C 

CONTROL 

335 

DIMENSION  JEM(830  ) ,KEN (830 ) , T ( 330 >  , 

336 

2  XIX (830) ,0(10,10) , ES ( 10, 10) ,SS(10,1Q),D3 (10)  , 

337 

3  P (10)  ,3(1155) ,XJ (830 ),YJ(830) ,ZJ(830),X< (830), YK(8 JO), 
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53  8 

539 

540 

54  i 

542 

543 

544 

545 

546 

547 

548 

549 

650 

651 
552 

55  3 

554 

555 

556 
55  7 

558 

559 

560 

561 

552 

553 

554 

555 
666 

557 
568 
559 

570 

571 

572 

573 

574 

575 

576 
5  77 

578 

579 
58  0 

581 

582 

583 
534 

585 

586 

587 
68  8 

589 

590 

591 

592 

593 
694 

595 

596 

69  7 
598 
539 

70  0 
701 


4  74(830) , AREA ( 830) 

*,MTV»E(830) 

COMMON  /2/  X  (231)  i  Z ( 231 ) ,  PL0TAX(220>, 

1PL0T9J(22)>  ,PL0T9<(223> , STRSOJ (223 > »STRSI J (220) , 5 TRSOK  (220 > , 

2ST RSTK ( 220 )  >NA9GLE 
COMM39/1/9BAYL,N3AYW 
COMMON/4/PLOTN 
C0MM39/5/CA, JPR03,NU9PR 
DIMENSION  A ( 13) 

399  FORMAT  (  8X,2H9N,  14X,  6HSTRS0J,  14X,6HSTRS0K,  14X,  5HSTRSIJ, 

114X , 5HSTRSIK  ,//) 

400  FORMAT (I10,4E2Q«5  ) 

420  FORMAT (20X14H»R09LEH  NUM8ER  147) 

452  FORMAT (99H  MEMBER  NO.  MOM.  AT  J  END  MOM.  AT  <  END  SHEAR 

2  AXIAL  > 

453  FORMAT (17, 1P4E17, 5) 

362  FORMAT ( 1H0) 

461  FORMAT (42H  THESE  ARE  THE  ELASTIC  NODAL  DISPLACEMENTS/) 

550  FORMAT (1H1) 

460  FORMAT (10X29H  THESE  ARE  THE  ELASTIC  FORCES//) 

456  FORMAT (16. IE 17.5.5 El 5.5) 

455  FORMAT (123H  900E  DELTA  X  DELTA  T  OELTA  Z 

1  THETA  X  THETA  Y  THETA  Z  ) 

PRINT  550 
PRINT  420,  NPR09 
“RINT  461' 

PRINT  455 

JR0  =  3 

1*0 

00  671  JR*5,N005,5 

JRO* JRO+1 

I*I»1 

XX1*3(JR-4)*R< 

YY1*3(JR-3)*RK 
ZZ1*3( JR-2) *R< 

XX=0(JR-1)*RK 

YY*0(JR)*RK 

571  PRINT  456, JRD.XXl, YYI,ZZ1,XX,YY 
NB AY3*N005/5 
'  »RINT  362 
PRINT  362 
PRINT  460 
PRINT  452 
■  JERK* 3 
00  41  MN=1,NM 
J6=JE9(MN) *5 
OD ( 1) =0 ( J6-4) 

0D(2)=0(J6-3) 

00 (3) *0 ( JS-2) 

00 (4) *0 ( Jo-1) 

00(5) *0 ( J5) 

K6=KEM(MN)*5 
00(6) =0 ( <5-4) 

00(7) *0 ( K5-3) 

00(8) *0 ( K5-2) 

00(  9) *0  (<6-l) 

00 (ID) *0 (K6> 

CALL  3MPAT (XJ,  YJ,ZJ,XK,YK,ZK,T,R,MN) 

ES(1,1)  =  AREA(MN>  *E/T(MN) 

£S (1,2 >*0.0 
£S(1,3)=0.0 
ES(2, t>=0.0 
E3 (3, 1) =0. 0 

ES(2,2)=4.0*E*XIX(MN) /T(M9) 
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732  ES (2*  3)  =  -£S(2*2)/2.3 

733  ES (3*  2) =ES (2,3) 

704  ES(3,3)=ES<2,2> 

735  CALL  NULT(3,10,3,SB, ES,B> 

736  C  P  =  ( S8)  *  0 

737  CALL  NJLT(3»l*13»P»S3»OD) 

738  AXIAl=P(l)*RK 

709  XXX=’<2)*RK 

710  YYY=P(3)*Ri< 

711  IF  MTYPE(MN>-2>  10,13,20 

712  10  JERK=JERK*1 

713  PLOTAX(JERi<)=AXIAL 

714  ?LOTMJ(J£RK)=XXX 

715  PLOTN< ( JERK) =YYY 

716  AA=AXIAL  /AREA(MN) 

717  '  STRS0J(JERK)=AA*XXX*30UmiX<MN» 

716  STRSt J(JERK)=AA-XXX*:iN/(XIX<MN)*RATMOO> 

719  ST.RS3<<JER<O=AA<-YYY*30UT/XIX(MN> 

720  STRSIK  (JERK) =AA-YYY*3IN/ (XIX(MN)*RATNOD) 

721  20  CONTINUE 

7>2  SHEAR=(YYY-XXX)/T(MN) 

723  PRINT  453,MN,XXX, YYY, SHEAR, AXIAL 

724  41  CONTINUE 

725  PRINT  550 

726  PRINT  399 

727  J=0 

728  30  43  MN=1,NM 

729  IP  MTYPE(MN)-2)  50,53,40 

730  50  J=JH 

731  •  PRINT  400,  MN,STRSOJ(J>  ,STRSO<( J) ,STRSIJ( J>  ,ST.RSI<  (J> 

732  40  CONTINUE 

733  C  THIS  ’ART  OF  THE  SUBROUTINE  CONTROLS  THE  ENTRY  OR  NON-ENTRY  INTO  THE 

734  C  PLOTTING  ROUTINES. 

735  C  IF  CONTROL  EQUALS  ZERO  NO  PLOTS  HILL  BE  GENERATED .  CONTROL  IS  READ  IN  AN 

736  C  F10.3  FORMAT. 

737  READ  500, CONTROL 

733  530  FORMAT (F10 .3) 

739  IFICONTROL.EQ. 0.)  GO  TO  42 

740  A  ( 1) =0 A 

741  A(2) =  JERK 

742  A  {  3) sNOAYL 

743  A (4) =NBA YW 

744  A  ( 5)  = NUHPR 

745  A (6) =NPROB 

746  A ( 7) =R< 

747  A { 8)  sNANG'.E 

748  A ( 9) = JPR03 

749  A(10)=?LOTN 

750  N3DE$=N005/5 

751  A (11) =N005 

752  ft ( 12) =NOQES 

F53  H13)=NM 

754  00  700  1=1,13 

755  PRINT  701, I, A(I) 

756  701  FORMAT (5  X, 3HI  =  , 1 5 , 5< , 6H A ( I >  =,F10.4» 

757  700  CONTINUE 

758  PRINT  704 

759  00  732  1=1, JERK 

760  PRINT  703,1, STRSOJ (I) , STR3I Jd ) , STRSOK (I ) , STRSI<< X ) , PLOTM J(I) , 

761  l?LOTN<d),PLOTAXd) 

75  2  70  4  FORMA  T  <7X,  1HI,  13X,5H3Jd),10X,5HIJd>,  10  X,5H0Kd>  ,10X,5HIKdl  , 

753  110X,5NMJ(I) , 10X , 5HMK ( I > , 10X ,5HAX d ) ) 

754  703  FORMAT ( 5 X, 15 , 7 (5 X , F13 . 4) ) 

765  702  CONTINUE 


68 


PRINT  706 
00  705  I=1,N035 
PRINT  707,I,Q(I),I)(I> 

706  FORMAT  (5X, 3H  I  ,10X,4 

707  FORMAT (5X, I5,10X, F10. 
705  CONTINUE 

PRINT  709 
00  70S  1=1, NODES 
PRINT  710,I,X(I>,Z(I> 

709  FORMA T  ( 5X ,  3H  I  ,10X,4 

710  FORMAT (5X,I5«10X«F10» 

708  CONTINUE 

817  IFUNIT,. 7)  817,818,813 

820  PRINT  804 

818  3UFFER0UT(7,1) (ATI), A 

819  CONTINUE 

300  IF  {JNIT, 7) 800,801,8] 
803  PRINT  804 

801  9UFFiR0UT<7,l) (STRSOJ 
304  FORMAT {*?ARITT  ERROR* I 

802  CONTINUE 

805  IFTUNIT, 7)805, 306, 807, 

803  PRINT  304 

806  3UFFER0UT<7,l)«STRSIJl 

807  CONTINUE 

809  IF  (UNIT,  7)809,810, 911, 

812  PRINT  804 

810  9UFFER0UTT 7,1) (3TRSO<( 

811  CONTINUE’ 

813  IFiUNIT, 7)813,314, 813, 
816  PRINT  804 

814  9JFF£TDUT(7,1> (STRSI<( 

815  CONTINUE 

821  IF(UNIT, 7)821, 822, 821, 

824  PRINT  804 

822  3JFFH0UT(7,IMP„0TNJ( 

823  CONTINUE 

825  IFUNIT,  7)825,  826, 827, 
628  PRINT  804 

826  9JFFiR0UT(7,l) (PLOTN<( 

827  CONTINUE 

829  IFCUNIT, 7)829,935. 831, 
332  .PRINT  804 

830  9UFFE  ROUT ( 7, 1) (°LOT4<( 

831  CONTINUE 

833  IF (UNIT, 7)833, 834, 835, 

836  PRINT  804 

534  CUFF- ROUT (7,1) (3(1) ,]( 
835  CONTINUE 

837  IF (UNI T, 7) 837,838,839, 
540  PRINT  904 

838  9UFFER0UT (7,1)  (X(1),X( 

839  CONTINUE 

841  IF  (UNIT,  7)  941,  842, 843, 

844  PRINT  804 

842  3UFF-R0UT(7,1) (Z(1),Z(| 

843  CONTINUE 

845  IF (UNIT,  7)  845,  846, 347, l 
848  PRINT  804 

846  BUFFiR0UT(7,l) (0(1), 3(t 
900  IF(UNIT,7)  900, 847, 847, ( 

847  CONTINUE 
42  CONTINUE 

ENO 


HQ(I),10X,4HO(I)) 
4, 10X, F10 .4) 


MX ( I) , 10X ,4HZ  (I) ) 
4,10X,F10.4> 


(13)) 

2,803 

(1) »5TRSOJ(JERK>) 


(1)  , 5T RSI J (JERK) > 


(1) »STRSOK(JERK>) 


( l ) , STRSI  <  ( JERiO ) 


(1) ,3LOTHJ(JERK>) 


(l),»LOTMK(JERK)) 


(1) , PLOTAX (JERK) ) 


( N005) ) 


(NODES)) 


(NOOES)) 


(N005) ) 
,847 
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3.  COMPUTER  PROGRAM  SAN1PLT 


1  PROGRAM  SAN1PLT (INPUT  , OUTPUT  »T  APE 7 , TAPES  »FILHPL) 

2  CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

3  cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx.xxxxxxxxxxxx 

4  C  THE  INPUTS  TO  THIS  ROUTINE  ARE  THE  TAPE  7  FROM  SAN1  AND  THE 

5  C  FOLLOWING  CAROS 

6  C  (DPLOIC ON -SIGNALS  LAST  SET  OF  PLOTS  WHEN  EQUAL  TO  Q(FIQ.O) 

7  C  (2)  L-LETERMINES  TYPE  OF  LOAOS(IIO) 

8  C  .  L  =1 - X  OIRECTION 

9  C  L=2 - Y  DIRECTION 

10  C  L=3 - 2  DIRECTION 

11  C  L s4 - IHfcTA  X  DIRECTION 

12  C  L=5 - THETA  Y  OIRECTION 

13  C  (3)NUHX — NUMBER  OF  NODES  FOR  CENTER  ARCH 

14  C  NUM2 — NUMBER  OF  NODES  FOR  CROWN  DEFLECTIONS (2110 ) 

15  C  ( 4)  NOO-— NODE  NUMBERS  FOR  CENTER  ARCH  (NUMX  OF  THEM) (110) (LINEAR  SOLUTION) 

16  C  (5)  NOC-NOOE  NUMBERS  FOR  CROWN  DEFLECTIONS  (NUMZ  OF  THEM)  (110) 

17  C  ( 6) NOO-NOOE  NUMBERS  FOR  CENTER  ARCH (NUMX  OF  THEM) ( I 10 ) (NONLINEAR  SOLUIICN) 

18  C  ( 7) SF-SC ALE  FACTOR  FOR  INSIOE  AND  OUTSIDE  STRESSES  (F1Q. 4) 

19  C  A  SET  CF  THE  AfcCVE  ARE  NECESSARY  FOR  EACH  SET  OF  PLOTS 

20  C  THE  LAST  DATA  CARO  SHOULD  BE  BLANK  TO  SIGNAL  THE  END  OF  THE  PLOTS 

21  CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

22  CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXX 

23  COMMON/l/AdOUO)  ,8(1000)  ,C(100Q)  ,0(1040)  ,E(13)  , 

24  1  F(IOOD)  ,6(1000)  , 

25  2  H(220),  0(220),  M220)  ,S  (231)  ,U  (231  >,  V  (1155)  ,W  ( 1155) 

26  COMMON/2/STRSOJ<1000) ,STRSU (1000) , STRSOK (1C 00 ) , STRSIK (1 000)  , 

27  1  PLOTMJ(220),PLOTMK<220),PLOTAX(220),X(231),Z(231), 

28  2  01(1155) ,X0(1155) ,0(1155) 

29  CGMM0N/3/CA,JERK,NBAYL,NBAYW,NUMPR,NPR0B»MANGLE, JPROB 

30  CALL  PLOTS (100  < ,OUH, 5) 

31  CALL  PMOOE  (3) 

o2  CALL  ROTATE  (270.) 

33  CALL  MAP  ( — 1 . , 11 . , — 1 . , 11 . , 0 . , 1 • ,4 • , 1 . ) 

34  C  A(l)-A<10Q0)=STRSOj 

35  C  B(l)-B(1000)*STRSIJ 

36  C  C  (1) -C (100  0) *STRSOK 

37  C  0(1) -0(1000) -STRSIK 

38  C  E  (l)-CA  E(2)=JtRK  E(3)-N8AYL  E(4)=NBAYW  E(5)*NUMPR  E(6)=NPROB 

39  C  E(7)=RK  E(8) sMANGLE  E(9)*JPROB  E(10)»PLOTN 

40  C  E  (11 )  =  N005  E (12)=NOOES  E(13)=NM 

41  C  H  (1) -H (220 )=PLOTMJ 

42  C  0  (1) -0 (220 ) =PLOTMK 

43  C  P  (1) -P(220)=PLCTAX 

44  C  0(1) -Q (1155) =Q 

45  C  S  (1) -S (231 )=X 

4b  C  U  (1) -U (231)-2 

47  C  V  (1) -V (1155)  =  01 

48  C  W  (1) -W il!55) =X  0 

49  4  FORMAT (8F10. 4) 

50  3  REAO  2.PL0TCCN 

51  2  FORMAT (F10.0) 

52  IF (PLOTCON  «£Q • 0. )  GO  TO  16 

53  817  IF(UMi, 7)817, 818, 819, 820 

54  820  PRINT  804 

55  818  BUFFERING, 1)  (E(l)  ,E(13)) 

56  910  IF(UNIT, 7)910, 819, 819, 819 

57  819  CONTINUE 

58  CA=£ (1)  $  JERK=E(2)  $  NBAYL=E ( 3)  %  NBAYH=E (4)  1  NUMPR=E(5) 

59  NPR08  =  E(6)  IRK  =  £  ( 7 )  $  MANGLE=E ( 8 )  $JPROB*E(9)  $PLOTN*E(10) 

60  NOD5=E (11)  S  NOOES=E (12)  $  NM=E(13) 
ol  800  IF  (UNIT, 7)800, 801, 802, 803 
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b2 

803 

PRINT  804 

b3 

801 

BUFF  ERIN <7 ,1) (All) ,A(JERK> ) 

64 

804 

FORMAT (*PARITY  ERROR*) 

b5 

802 

CONTINUE 

66 

805 

IF (UNI It 7) 80 5 ,806 ,807,808 

67 

808 

PRINT  804 

68 

80b 

HUFFERIN (7  »1>  (8(1) ,B (JERK) ) 

69 

807 

CONTINUE 

70 

809 

IF  (UNIT, 7) 809, 81 0,81 1,812 

71 

812 

PRINT  804 

72 

810 

BUFFERIN (7 ,1) (C(l) ,C(JEkK)) 

73 

811 

CONTINUE 

74 

813 

IF (UNIT, 7) 813 ,814,815, 8ib 

75 

816 

PRINT  804 

76 

814 

BUFFER  IN (7,1) (C(l) ,0 (JERK) ) 

77 

815 

CONTINUE 

78 

821 

IF (UNI 1,7)  821,822, 823,  824 

79 

824 

PRINT  804 

80 

822 

BUFFERIN  (7,1)  (HI)  ,H(JERK) ) 

81 

823 

CONTINUE 

82 

825 

IF  (UNIT, 7) 825, 826, 827, 828 

83 

828 

PRINT  804 

84 

826 

BUFFERIN (7,1) (0(1) ,0(JERK) ) 

85 

827 

CONTINUE 

86 

829 

IF  (UNIT, 7) 829,030,831,832 

87 

832 

PRINT  804 

d8 

830 

0UFFERIN  (7 ,1)  (Ml)  ,P(JERK>> 

89 

831 

CONTINUE 

90 

'  833 

IF (UNI {,7)833, 834, 835, 836 

91 

836 

PRINT  804 

92 

834 

BUFFERIN  (7 ,1 ) (U(l> ,Q(N0U5)) 

93 

835 

CONTINUE 

94 

837 

IF (UNIT, 7) 837, 838, 839, 840 

95 

840 

PRINT  804 

96 

838 

BUFFERIN (7,1) (S(l) ,S(NOOES)> 

97 

839 

CONTINUE 

98 

841 

IF(UNIT , 7) 841,642,843,844 

99 

844 

PRINT  804 

100 

842 

BUFFERIN (7,1) (U(l) ,U(NOUES)) 

lul 

843 

CONTINUE 

102 

845 

IF  (UN 1 1,7) 845, 846, 847, 846 

103 

848 

PRINT  e04 

104 

846 

BUFFERIN (7,1) (V(l» ,V(N0D5) )  ‘ 

105 

847 

CONTINUE 

106 

849 

IF  (UNIT, 7)  849,850,851,852 

107 

852 

PRINT  e04 

108 

850 

BUFFERIN  (7,1)  (Mi>  ,N(N0D5) ) 

109 

900 

IF (UNIT, 7) 90 0,851, 851, 851 

110 

851 

CONTINUE 

111 

UC  10  1»1,NQ05 

112 

S  TRSO J (I ) *  A ( I ) 

113 

S  TRSI J (I )  =  8 ( I ) 

114 

STRSOK(I)=C(I) 

115 

STRSIK (I)=0(I) 

116 

10 

CONTINUE 

117 

00  11  1=1, JERK 

118 

PUOTMJ(I)=H(I) 

119 

PU0TMK(I)=0(I) 

120 

PUOT  A  X ( I )  =  P (I) 

121 

11 

CONTINUE 

122 

OC  12  I=l,NOD5 

123 

Q  (I)  =Q (I) 

124 

01  ( I )  =  V  C  I) 

125 

XO  ( I  )  =  M  I) 

126  12  CONTINUE 

127  CO  13  I=l,NOQES 

126  X  tl »  =S  ( I > 

129  Z(I)=U(I> 

130  13  CONTINUE 

131  00  700  1=1,13 

132  PRINT  701, I, Ed) 

133  701  FORMAT (5X.3HI  = , X 5 , 5X , 6HA < I )  =,F10.4) 

134  7GC  CONTINUE 

135  PRINT  704 

136  OC  702  1=1, JERK 

137  PRINT  703,I,STRSOJ(I) , STRSIJ (I ) , STRSUK (I) ,STRSIK(I)  ,PL0TNJ(I) , 

138  lPLOTMK(I),PtOTAX(I) 

139  70  4  FORMAT  <7X, 1HI , 10X ,5HOJ ( I) , 10 X, 5HIJ( I) ,  1QX,  5HCK  (I)  ,  IflX,  5HIMI )  , 

14(1  110X,5HKJ(I)  ,10X,5HMK(I)  ,10X,5HAX(I>) 

141  703  FORMAT (5X » 15 » 7  (5X ,F10« 4) ) 

142  702  CONTINUE 

143  PRINT  706 

144  OC  705  I=1,N005 

145  PRINT  707,1,0(1) ,01(1)  ,X0<I) 

146  706  F ORM AT (5 X , 3H  I  , 10 X,  4riQ ( I) , 10X , 4HD ( I) , 10X , 5HX0 ( I > > 

147  707  FORMAT(5X,I5,lCX,Flu.4,lQX,F10.4,i0X,F10.4> 

146  705  CONTINUE. 

149  PRINT  709 

150  00  708  1=1 iNOOES 

151  PRINT  71 0 , I , X ( 1) , Z (I ) 

152  709  FORMAT  (5X.3H  I  ,li»X,4HX(I)  ,10X,4hZ(I)> 

1»3  711  FCRMAT(5X,I5,1CX,F10.4,1QX,F1Q.4I 

154  708  CONTINUE 

155  CAUU  TUOP(RK) 

156  CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

197  CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
158  C  18-CONTROL  CRCRN  OEFUECTIONS  ANO  CENTER  ARCH  PLOTS 

169  C  I8=0-80TH  PLOTS 

160  C  IO  =  l-CEN TER  ARCH  ONLY 

lbl  C  PLOTN-CONTROLS  LABELS  FOR  CENTER  ARCH  PLOTS 

162  C  PLOTN=0. -LINEAR  SOLUTION 

160  C  PLOTN=l. -NON-LINEAR  SOLUTION 

164  CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXA 

165  CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

166  PLOTN=0. 

167  16=0 

168  CALL  ARCHPLT(RK,X,Z, 01, MANGLE, PLOTN, 16) 

169  PLOTN= 1 • 

170  18=1 

171  CALL  ARCHPLTU. ,X,Z,XJ, MANGLE, PLOTN, 18) 

172  CALL  STRESSP 

173  GO  TO  3 

174  16  CALL  PLOT(0.,0.,40) 

175  ENO 

176  SUBROUTINE  TLOF(RK) 

177  COMMON/2/STRSOJdOOO)  ,  STRSIJ  (100  0)  .STRSOK  (1000)  ,STRSIK  (1000)  , 

178  1  PLOTKJ  (220  )  ,  PLOTMM220  )  »PLOTAX  (220  )  »  X(231)  »Z  (231)  » 

179  2  01  (1155), XUU155)  ,0(1155) 

180  COMMON/3/CA,JEfiK,NUAYL , N8A YW , NUMPR, NPR08, MANGLE, JPRCB 

1«1  CALL  SLAOR(UNEL,ONEW,XL,YH,NBL,NBH) 

182  IF(NBAYH.GT.NBAYL)  GO  TO  1 

lo3  N8AY  L I  =  N0A  YL 

184  N0AYKI=N8AYH 

185  GO  TO  5000 

186  1  NBAYLI=N8AYH  {NBA YHI=N8A YL 

187  GO  TO  5000 

168  CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
189  CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
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190  C  PLOT  NUMBER  2 

191  C  MOMENTS  AT  J-ENO  AND  K-ENO 

192  CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

193  cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

194  5000  XA-0.0 

195  CALL  PLOT  (0. ,0  .,3) 

196  CALL  PLOT ( 0. , 0 • ,2) 

197  00  10  I=i,N8AYLI 

190  CALL  PLO  T(  XA 1 0 • , YH  ,1) 

199  CALL  PLOT(XA+ONEL,YW,i) 

200  CALL  PLOT(XA*ON£L,0. ,1) 

201  CALL  PLOT (XA,Q  <  ,1) 

202  CALL  PLOT <  XA  +  ONEL *  0. ,3) 

203  CALL  PLOT(XA*ON£L,0« ,2) 

204  XA=XAKNEL 

205  10  CONTINUE 

2L6  YA-0 . 

207  LI=NBAYWI/2. 

206  00  11  J=1,LI 

209  CALL  PLOT(0.,YA,3> 

210  CALL  PLOT(Q. , YA+ONEH, 3) 

211  CALL  PLOT(0.,YA-fON£M,2) 

212  CALL  PLOT(XL,YA+ONEW,l) 

213  CALL  PLOT (XL  t Y A+2*0NEM  ,o ) 

214  CALL  FL0T(XL,YA+2*0NEN,2> 

215  CALL  PLOT (0. , YA+2*ON£H,l) 

216  YA=YA+2*ONEW 

217  11  CONTINUE 

216  XAsXL  IYA=Q.  JL=1 

219  C  PLOT  NUMBER  TWO. 

220  C  0EFINING  MOMENTS  AT  J-ENO  AND  K-fcND  OF  ELEMENTS. 

221  00  110  1=1, N8L 

222  00  100  K=1 i NBA YNI 

223  C  PLCTMJ  =  MOMENT  AT  JENO  OF  ELEMENT 

224  C  PLOTMK  =  MOMENT  AT  KEHU  OF  ELEMENT 

225  CALL  NLM3ER(XA-.28«YA+ (.l*ON£H) f.l2,PLOTMJ(L),0.«4hF8«2) 

226  CALL  NLM6ER(XA-.28,YAM.8*ON£W>,.12,PLOTMK(L),0.,4HFB.2> 

2c7  L=L+1 

22 «  100  YA=YA»CNEH 

229  YA  =  0. 

230  XA=XA-CNEL 

231  110  CONTINUE 

232  XFOS= ( XL/2.) -1  • 

233  CALL  PLOT (XPOS ,-.35,3) 

234  CALL  PLCT( XPOS ,-.35,2) 

235  CALL  SYMBOL (XPOS, -.35* .12 , 21HM0MENTS  (  INCH-KIPS  >,u.,21> 

236  CALL  SYM80L(XPCS+.b,-.5U,.lQ ,7HPR08LEM ,0 . , 7) 

237  CALLNUK8ER(XPOS*1.25,-.50,.10, NPROB, 0. , 2HI5) 

236  CALL  PLOT(XL+5.,0.,-3) 

239  GC  TO  €00 

240  500  X A=0  • 

241  YA=0. 

242  L*1 

243  DO  121  1=1 *NBA YLI 

244  OC  ill  K=1  ,N8H 

2<*5  CALL  NLMBER(XA+(,15*ON£L> ,YA-(.2*0N£N> ,.12,PL0TMJ(L) ,90. ,4HF8.2> 

246  CALL  NUMBER ( XA+ ( .85*0NEL)  , YA-( ,2*ONEW)  ,.  12, PLOTMK  (L »,  90.  ,4HF8.2J 

247  L=L*1 

248  111  Y A=Y A+CNEM 

249  YA=0. 

250  XA=XA+CNEL 

251  121  CONTINUE 

252  XPOS=(XL/2.) -1. 

253  CALL  FLCT( XPOS ,-.35,3) 
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254 

255 

256 

257 
250 

259 

260 
261 
262 

263 

264 
2o5 
266 
267 

2b0 

269 

270 

271 

272 

273 

274 

275 

276 

277 

278 

279 
200 
201 
202 

203 

204 

205 

206 
207 
200 
209 

290 

291 

292 

293 

294 

295 

296 

297 
290 
299 
30C 

301 

302 

303 

304 

305 

306 

307 

308 

309 

310 

311 

312 

313 

314 

315 

316 

317 


CALL  PL0T<XP0S,-.35,2> 

CALL  SYMBOL<XPCS,-.35,.12,21HMOMENTS  (  INCH-KIPS  >»0.,21> 

CALL  SYN80L(XPOS+.5,-.50,.10,7HPROBLEK,Q.,7> 

'C  ALL  NUMBER  <  XPOS+ 1.25,  -.50,  .10,  NPROB,0.  ,2HI5) 

CALL  PLOT(XL+5.,0.,-3) 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
C  PLOT  NUMBER  3 

c  AXIAL  STRESSES 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
600  X  A=Q • 0 

CALL  PLOT (  0. «  0 • <  3) 
call  FLCTt  0 • , 0  •  *2) 

00  12  1=1 ,  NBA YLI 
CALL  PLOT(XA+O..YH,1) 

CALL  PLOT(XA*ONEL»YN,l) 

CALL  PLOT (XA+ONEL»0. » 1 ) 

CALL  PLCT(Xrt,Q.,l) 

CAlL  PIOT<XA*ONEL,0. »J) 

CALL  PLOT<XA+ONEL,G. ,2) 

XA=X A+CNEL 

12  CONTINUE 

YA=0 •  , 

lI=NBAYHI/2. 

00  13  J=1,LI 
CALL  PLOT (Q.,YA,3) 

CALL  PLOT  (  0  .  , Y A*ONEH , 3 ) 

CALL  FLOT<0.,YA*ON£W,2> 

CALL  PLOT(XL,YA+ONEH,l) 

CALL  PL0T(XL,YA+2*0NEH,3>  '  ' 

CALL  PLOT  (XL,YA+2*ONEH,2) 

CALL  PLOT(0.,YA«-2*ONEH,1> 

YA=YA+2*ONEN 

13  CONTINUE 

C  CHECK  FOR  FULL  OR  HALF  ARCH. 

IF(MAKGLF,.EU.2>  GO  TO  700 
XA=XL  1 Y  A=  0  •  IL  =  1 
C  PLOT  NUMBER  THREE. 

C  OEFINING  AXIAL  FORCE  FOR  EACH  ELEMENT. 

00  210  1=1, NBL 
00  200  K  =  1 , NBA YNI 
C  PLOT  AX  =  AXIAL  FORCE 

CALL  NUMOER<XA-.20,YA+(.5*ONEH),.12,PLOrAX(L)»O.»4HF5,2) 

L  =L+  1 

200.  Y A  =  Y A ♦  CN  EH 
YA  =  0  . 

XA=X A-CNEL 
210  CONTINUE 

XPOS=  <  XL /2 . ) -1 . 

CALL  PLOT ( XPOS ,-.35,3) 

CALL  PLOT(XPOS,-.35,2) 

CALL  SYHB0L(XPCS,-.35, .12.21HAXIAL  FORCES  (  KIPS  ),0.,21) 

CALL  SYMBOL(XPCS+.5,-.50,.10,7HPROBLEH,0.»7) 

CALLNUMOERTXPOS* 1. 25 ,-.50 , . 10 , NPRO0 , 0. ,2HIS) 

CALL  PLOT(XL+5.,0.,-3) 

GO  TO  aoo 
706  XA=0. 

YA  =  0. 

L  =1 

00  221  I =1 jNBAYLI 
00  211  K  =  1  ,NEH 

CALL  NUMBER (XA+. 5, YA- < .2*0NEW) , . 12, PLOT AX <L> ,90.,4hF5.2) 

L=L+ 1 


74 


318 

319 

320 

321 

322 

323 

324 

325 

326 
32/ 

328 

329 
331' 

331 

332 

333 

334 

335 

336 

337 

338 

339 
24C 

341 

342 

343 

344 

345 

346 

347 

348 

349 

350 

351 

352 

353 

354 

355 

356 

357 

358 

359 

360 

361 

362 

363 

364 

365 

366 
3o7 
366 

369 

370 

371 

372 

373 

374 

375 

376 

377 

378 

379 

380 
3di 


211  YA»YA*CNEH 
YA=0. 

X A=X A+CNEL 
221  CONTINUE 

XP0S=(XL/2.I-1. 

CAUL  PL0T(XP0S,-.35,3> 

CALL  PL0T<XP0S,-'.35,2> 

CAUL  SYM80L<XP0S,-.35, . 12* 21HAXIAL  FORCES  (  KIPS  ),0.,21) 

CAUL  SYMBOL(XPCS+.5,-.5fc,.10,7HPkOBL£M,Q.,7> 

CALL NUMBER (XPQS+1.25»-.50» . 1C » NPROB ,0 . >  2HI5) 

CAUL  PLOT<XL*5.,0.,-3> 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
C  PLOT  NUMBER  4 

C  INSIDE  STRESSES  AT  J-ENU  AND  K-END 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXA 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXX 
800  XA=3  *0 

CALL  PLOT (0««Q.»J) 

CALL  PLOT(O..Q.,2» 

DO  14  1=1* NBAYLI 
CALL  PLOT(XAtO.,YH,l* 

CALL  PLCT<XA*ON£L,YW,il 
CALL  PLOT<XA*ON£L»U. ,1> 

CALL  PLOT ( XA« 0  a « 1) 

CAUL  PLOT(XA*ONELtli.  ,3» 

CALL  PLOT ( XA+ONEL •  0.  •  21 
X A»XA+CNEL 

14  CONTINUE 

YA=0. 

LI=N8AYWI/2. 

DC  15  jsj,,LI 
CALL  PLOT  <0*tYA«3) 

CALL  PLOT ( 0.  « Y  A+ONfcH  » 3 ) 

CALL  PLOT (0.  «  YA40NEH f 2) 

CALL  PLOT(XL,YA+ONtH,l> 

CALL  PLOT  (XL  *  YA,>2*0NEH»3) 

CALL  PL0T(XL,YA42*0NEH,2) 

CALL  PLOT(U..YA+2*ONEH,l) 

YA=YA*2*ONEW 

15  CONTINUE 

C  CHECK  FOR  FULL  OR  HALF  ARCH. 

IF(MANGLE.EU«2)  60  TO  900 
XAsXL  1YA*0.  SL  =  1  • 

C  PLOT  NUMBER  FOUR. 

C  DEFINING  INSIDE  STRESSES  AT  J-END  ANO  K-ENO  CF  CLEMENT. 

00  310  1=1, N8L 
00  200  KaltNQAYHI 

C  STRSIJsINSIDE  STRESS  AT  J-END.  ' 

C  STRSIKsiNSIDE  STRESS  AT  K-END. 

CALL  NLMBER<XA-.28,YAM.i*0NEW)  ,  .12  .STRSIJ  (L)  .0.  , 4HF7.2) 

CAUL  NUMBER (XA-. 28 »Y A+ (. 8*0NEH > ,«12,STRSIK(L) • 0. *4HF7.2) 

L=L+ 1 

300  YAsYAfCNEH 
YA=0. 

XAaXA-CNEL 
310  CCNTINLE 

XPOSs (XL72 . )-2 , 

CALL  PLOT ( XPOS »-. 35, 3) 

CALL  PLOT ( XPOS  »-.35,2) 

CALL  SYMBOL  ( XPOS i -.35 » . 12  *43HC0HPARAT 1VE  STRESSES  CN  INSIDE  FIBERS 
1  <KSI) .0.,43)- 

CALL  SYMBOL(XFCS+1.5,-.50,.10,7HPROELEM,0.,7> 

CAUL  NUMBER (XPOS+2 .25 50  ,  •  10 »NPROB ,0 . , 2HI5 ) 


75 


342 

CALL  PLOT (XL+5.,0.  »”3) 

343 

GO  TO  1000 

344 

900 

XA  =  0. 

345 

YA=0. 

346 

L  =  1 

3a? 

00  321  I=l,rJBAYLI 

348 

oc'  311  K=1,NBX 

349 

CALL  NUMBER (XA+(.15*0NEL) ,YA-(.2*ONEW> ,.12,STRS1JU),90.,4HF7.2> 

390 

CALL  NLMBER(XA+(.85»0NEL) , YA- < .2*0NEH> , . 12 .STRS1K <L ) ,90. ,4HF 7 .2) 

391 

L=L»1 

392 

311 

YA=YA*CNEW 

393 

YA=0. 

394 

X A=X A+CNEL 

395 

321 

CONTINUE 

396 

XPGS=OL/2.>-2. 

397 

CALL  PL0T<XP0S,-.35,3) 

398 

CALL  PLOT (XPOS»-.35,2) 

399 

CALL  SYMOOL(XPCS,-.35,.12,43HCOMPARATlVt  STRESSES  CK  INSIDE  PieERS 

400 

1  (KSI),0.,43) 

401 

CALL  SYMBOL ( XPCS+ 1 .5 ,-.50».10» 7HPR0BLEM»C • « 7) 

402 

CALL  NLMBER(XPOS*2.25,-.5Q,.10,NFROE,0. »2HI5) 

403 

CALL  PL0TCXL*5.,Q.,-3> 

4i)4 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

405 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

406 

c 

PLOT  NLNBER  5 

407 

c 

OUTSIDE  STRESSES  AT  J-END  ANO  K-ENO 

408 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

409 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

410 

1000 

XA=0.0 

411 

CALL  PLOT(0.,0.,3) 

412 

CALL  PLOnO.,0.,2) 

413 

00  16  1  =  1 » NOA  YL I 

414 

CALL  PLOT(XA*Q. ,YH,1) 

415 

CALL  PLOT(XA*ONEL,YW,l) 

416 

CALL  PLOT(XA+ONEL,0. ,1) 

417 

CALL  PLOT(XA,0 •»!) 

418 

CALL  PLOT(XA+ONEL,0.,3» 

419 

CALL  F LO T ( XA+ONEL ,  0  • ,2) 

420 

X A=X A+CNEL 

421 

16 

CONTINUE 

422 

YA=  0 . 

423 

LI=N0AYHI/2. 

424 

OC  17  „=1*LI 

425 

CALL  PLOT (0. ,YA,3>  • 

426 

CALL  PLOT<0.  ,YA4-ONEW,3) 

427 

CALL  PLOT (0. ,YA+ONEW,2» 

428 

CALL  PLOT(XL,YA*ON£W»l) 

429 

CALL  PLOT (XL, Y A+2*0NEH ,3 > 

430 

CALL  PLOT(XL,YA  +  2*ONEW,;.') 

431 

CALL  PLOT ( 0, , YA+2*ONErt,l) 

432 

Y  A=Y  A  +  2*  ONEH 

433 

17 

CONTINUE 

434 

c 

CHECK  FOR  FULL  OR  HALF  ARCH. 

435 

IF(HANGlE.EQ.2)  GO  TO  1100 

436 

X A  =  XL  IYA=0.  1L  =1 

437 

c 

PLOT  NLHOER  FIVE. 

438 

c 

DEFINING  OUTSICE  STRESSES  AT  J-END  ANO  K-ENC  OF  ELEMENT. 

439 

00  41C  1=1, N8L 

440 

00  400  K=1 ,NBAYHI 

441 

c 

STRSOj=OUTSIOE  STRESS  AT  J-END. 

442 

c 

STRSOK=OUTSIOE  STRESS  AT  K-ENO. 

443 

CALL  NUMBER (XA-. 28, YA+(. l*ONEH) , . 12 , ST RSO J (L > , 0 . .4HF7.2) 

444 

CALL  NUMBER (XA-. 28, YA+(-8*0NEW> , . 12 ,STRSOK (L ) ,0. ,4hF7.2) 

445 

L  =L  +  1 
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<*46 

447 

446 

449 

450 

451 

452 

453 

454 
466 

456 

457 

458 

459 
46b 

461 

462 

463 

464 
4o5 

466 

467 

468 

469 
47b 

471 

472 

473 

474 

475 

476 

477 

478 

479 
48b 
481 

462 

463 
484 

465 

486 

487 

488 

489 

490 

491 

492 

493 

494 

495 

496 

497 

498 
459 

500 

501 

502 

503 

504 

505 

506 

507 

508 

509 


40b  Y A=YA»CNEH  > 

YA=0. 

XA=X A-CNEL 
41b  CONTINUE 

XP0S=(XL/2.l-2. 

CALL  PL0T(XP0S,-.35,3) 

CALL  Pl0T(XP0S,-.35,2) 

CALL  SYMBOL (XPCS,-. 35,. 12, 44HC0MPARATIVE  STRESSES  ON  OUTSIDE  FIBER 
IS  (KSI) ,0.  ,44) 

call  symbol (xpcs+i.5,-.50,.iq , 7h problem, o • ,7) 

CALL  NLMBER<  XP OS +  2. 25, -.50, .10 , NPROE  ,b • , 2HI5 ) 
call  plot <xl+  5 . , o. ,-3> 

GO  TO  1200 
110b  XA  =  0. 

YA=0. 

L=1 

OC  421  1=1 ,NSAYLI 
DC  411  K=l,N8h 

CALL  NUM6ER(XA*(.i5*ONEL>,VA-(.2*ONEN>,.i2,STRSOJIL>,90.,4HF7.2> 

CALL  NUMBER ( XA-M  . 85*0N£L> , YA- ( . 2*0NEN> ,.12,STRS0K(L>,9b. ,4HF7,2> 

L=L*1 

411  YA=YA+CNEH 
YA=Q. 

XA=XA+CNEL 
421  CONTINUE 

XPOS=(>L/2.)-2. 

CALL  PLO  T<  XPOS  *-.35,3) 

CALL  PLO} <XP0S,-.35,2J 

CALL  SYMI)OL(XPCS»-.35,.12 *44 HCOMPARA I I VE  STRESSES  CN  OUTSIDE  FIBER 
IS  (KSI), 0.  ,44) 

CALL.  SYMEOL( XPCS*1 .5 ,50 , . 10 , 7HPROULLK, 0 • , 7 ) 

CALL  XLMEER(XPCS«-2.25,-.50,  .10*NFROe,0.  ,2HIE> 

CALL  PLOT(XL*5.,0. ,-3) 

1200  CONTINUE 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
C  PLOT  NLMOER  6 

C  LOADING  VALUES 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
REA02  ,L 
2  FORMAT (110) 

X A=0  •  0 

CALL  PLOT (0. ,0 •  ,3) 

CALL  PLOT(0. ,0.,2) 

CO  18  1=1, NBA YLI 
CALL  PLOT(XA  +  0.,YH,’1> 

CALL  PlOT(XA*OKEL,YW,1) 

CALL  PbOT(XA*OKEL,u. ,1) 

CALL  PLOT(XA,0.,1) 

CALL  PLOT ( XA  +  OKEL ,0.  ,3 ) 

CALL  FLOT(XA  +  ONEL,0.  ,2) 

X A=X A+CNEL 
18  CONTINUE 
YA*Q  , 

LI=N8AYWI/2. 

'  OC  19  o=i, LI 

CALL  PL0T(0. ,YA,3) 

CALL  PLOT(0., YA+ONEH ,3) 

CALL  PLOT (0. * Y A+ONEW ,  2  ) 

CALL  PL0T(XL,  YA4-ONEW,l) 

CALL  PL0T(XL,YA+2*0NEW,3) 

CALL  PL0T(XL, YA*2*ONEW,2> 

CALL  PLOT(0. ,YA*2*ONEW,l) 
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510 

YA=YA+2*ONEH 

511 

19 

CONTINLE 

512 

C 

CHECK  FOR  FULL  OR  HALF  ARCH. 

513 

IF (MANGLE. EQ. 2 )  GO  TO  1400 

514 

NBAYRII  =  NBAYhIU 

515 

X  A=XL  JYA=0. 

516 

DO  521  J=1,N8L 

517 

OO  511  K=1,N8AYWII 

516 

CALL  NLMBER(XA-.26,YA* <.1*0NEH) ,. 12,0  CL) ,0.,4HF6.2> 

519 

L~L*5 

520 

511 

Y  A=Y  A+QNEW 

521 

YA  =  0. 

522 

XAsXA-CNEL 

523 

521 

CONTINLE 

524 

XPOS= (XL  72. >“.75 

525 

CALL  PLOT  (XPOS  .-.35. 3) 

526 

CALL  PLOT ( XPOS  35. 2) 

527 

CALL  SYMBOL  ( XPOS , 3 5 , . 12 . 1 6HL0A01NG  (  KIPS  >,u.,16) 

526 

CALL  SYMBOL (XPOS+.25, -.50 ». 10 , 7HPROBLEH.O . .7) 

529 

CALL  NLMBER<XPOS>1.0  *“.50* . 10.NPROB.Q. ,2HI5) 

530 

CALL  PLOT(XL*5.,0.,-3> 

531 

GO  TO  1600 

532 

1400 

CONTINLE 

533 

XA=0. 

534 

Y  A=Q  . 

535 

00  621  1=1 .NBAYLI 

536 

00  611  K=1,NBH 

537 

CALL  NLMBER(XA+(.15*0NEL> , YA“ < .2*ONEVO ,.12»0(L> ,90.,4HF6.2> 

538 

L=Lt5 

539 

6il 

XA=XAv'CNEL 

540 

X  A=0  . 

541 

621 

YA-YA+CNEH 

542 

XPOS=(XL/2.1-.75 

543 

CALL  PL0T(XP0S,-.35,3) 

544 

CALL  PLOT (XPOS .“.35. 2) 

545 

CALL  SYMB0LUP0S.-.35,  .12. 16HL0A0ING  (  KIPS  )»0.,16> 

546 

CALL  SYMBOL ( XPOS+ . 25 , -  .  5 1 1 . 10 , 7HPR0BLt M , 0 . .7) 

547 

CALL  KLHBER(XPCSn.O ,-.50, .10.NPROB.0. .2HI5) 

548 

CALL  PLOT(XL+5. ,0. ,-3) 

549 

1600 

CONTINLE 

560 

return 

551 

EHO 

552 

SUOROUTINE  SLALR(ONEL,CNEW,XL,YW.NBL,NEH) 

553 

C 

ONEL  =  ONE  INCRtMENT  IN  X  OIRECTION 

554 

C 

ONEW  =  ONE  INCREMENT  IN  Y  OIRECTION 

555 

C 

XL  =  TOTAL  LENGTH  X  OIRECTION 

556 

C 

Yfc  =  TOTAL  LENGTH  Y  OIRECTION 

557 

C 

NBAYH  =  NUHBER  OF  BAYS  HIDE 

558 

C 

N BAY L  =  NUMBER  OF  BAYS  LONG 

559 

C 

NBl,  =  NUMBER  CF  X  INCREMENTS  PLUS  1 

560 

C 

NBH  =  NUMBER  OF  Y  INCREMENTS  PLUS  1 

561 

0 

NSAYLI  =  NUMBER  OF  X  INCREMENTS 

562 

C 

NBAYWI  =  NUMBER  OF  Y  INCREMENTS 

563 

COMMON/3/CA.JERK.NBAYL , NBA YH, NUMPR, NPROB , MANGLE. JPRC8 

564 

2 

FORHAT (2110) 

565 

IF(NBAYW.GT.NBAYL)  GO  TO  1 

566 

c 

CHECK  TO  SEE  WHICH  UIMENS'ION  IS  LARGER 

567 

c 

SET  UP  SCALES  TO  A  10  INCH  VERTICAL  SCALE 

568 

Y  SCALE  =  1 0 . 

569 

X SCALE3 ( 10 ./NBAYH) *NBAYL 

570 

ONEL  =  XSCALE/NBAYL  $XL  =  0NEL*N8AYL 

571 

ONEW=YSCAlE/NBAYH  $YW=ONEH*NBAYW 

572 

N8AYLI=N0AYL 

573 

NUAYW1=NUAYW 
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574 

GO  TO  5000 

575 

1 

YSCAL£=10. 

576 

XSCALE=<10./NBAVL)*NBAVW 

57.7 

ONEL=XSC ALE/NBA YH  tXL=ONEL*NBAYN 

57B 

ONEW=YSCALE/NBAYL  $  YW=0NEW*N8AYL 

579 

N6AYLI-NBAYW 

580 

N8AYW1=NBAYL 

5Bi 

cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

562 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

5B3 

c 

PLOT  NUMBER  1 

584 

c 

GRIO  SYSTEM 

565 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX*XXXXXXXXXXXXXXXXXXXX 

566 

UXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

587 

5000 

X=0. 

568 

c 

ACTUAL  ORAWING  OF  STRUCTURE  WITH  ELEMENTS 

589 

CALL  PLOT ( 0 • • 0 • • 3) 

596 

CALL  PLCT<0. ,0.,2) 

591 

00  10  I=ltNBA YL1 

592 

CALL  PLOT < X+0 • tYH *  1) 

593 

CALL  PLOHX+ONEL,YH,i) 

599 

CALL  PLOT<X+ONEL,0.,1) 

595 

CALL  PLOT<X+0.,0.«1) 

596 

CALL  PLOT  <  X+O NEL  1 0  .  i  3) 

597 

CALL  PLOT ( X+ONEL 1 0 . » 2) 

598 

X=X+ON6L 

599 

10 

CONTINUE 

6t»0 

Y  =  0» 

601 

L=NBAYhI/2. 

602 

UO  11  v*l,L 

6u3 

CALL  PLOT (0 • • Y+C • |3) 

604 

CALL  PLCT(0.,Y«ONEN>3) 

605 

CALL  PLOT (0*  , Y+ONEWj2)  \ 

606 

CALL  PLOT<XL,Y+ONEW,ll 

607 

CALL  PLOT ( XL • Y+2*0NEH«3) 

608 

CALL  PLOT ( XL  *  Y+2*0NEH*  2) 

609 

CALL  PLOT<0.,Y+2*ONEH,1) 

610 

YsY+2+ONEH 

611 

11 

CONTINUE 

612 

T=NBAYXI/2.  $  LT=T  12=LT 

613 

IF(T.GT.2)GO  TC  16 

619 

102 

X=3.  S  V- 0  ■ 

615 

CALL  FLCT<Q.,Q.,3) 

616 

CALL  PLOT (0* |0»»2) 

617 

00  12  1= It  NBA YLI 

618 

IK«NBA YNI/2 

619 

00  20  K=lt IK 

620 

CALL  PLOT<X  +  ONELtY«-ONtH,l) 

621 

CALL  PLOT<X+0, ,Y+2*ONEH,l> 

622 

Y=Y+2+CNEN 

623 

20 

CONTINUE 

624 

X=X+0.  iY=YN 

625 

CALL  PLOT ( X+ONEL  >  YWt  3) 

626 

CALL  FLOT(X+ON£L,YH|2) 

62/ 

M=NBAYm/2. 

628 

00  15  «=1,M  . 

629 

CALL  PLCT(X+0.,Y-ONEH,1) 

630 

CALL  PL0t<X+0NELiY”2*0NEW,l> 

631 

Y=Y-2*CNt-H 

6o2 

15 

CONTINUE 

633 

X=X+ONEL  $V=0. 

634 

CALL  PLOT (XtY «  3) 

635 

CALL  PLOT (XtY  t2) 

636 

12 

CONTINUE 

637 

GO  TO  22 
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638 

16 

X*0 • IY*0  • 

639 

CALL  PLOT ( 0. , 0 . *  3) 

640 

CALL  PLOT ( 0  « , 0  •  ,2) 

641 

00  17  l=i,NBAYLI 

642 

IK=NBAYWI/2 

643 

00  19  K= 1 , IK 

6*4 

CALL  PL0T(Xf0NEL,Y4-0Nt«,l» 

645 

CALL  PLOT (X*Y+2*ON£H,l) 

646 

Y=Y+2.«ONEW 

647 

19 

CONTINLE 

648 

CALL  PLOT (XJ-OXELfY+ONtHfl) 

649 

X=X*0.  $Y=YW 

650 

CALL  PLOT <X,YH,3) 

651 

CALL  PL0T(X,YH,2> 

652 

M=NBAYhI/2 

653 

00  21  J=1,U 

654 

CALL  PLOT(X*ON£L»Y-ONtNfl) 

655 

CALL  PLOT(X,Y-2.*ONEW,i» 

656  ' 

Y=Y-2.*ONEW 

657 

21 

CONTINLE 

658 

CALL  PLOT(X+ONEL,Q.,l» 

659 

X  =  X+  0  NEL  $  Y  =  0. 

6b0 

17 

CONTINLE 

6bl 

22 

CONTINLE 

662 

C 

NUMBERING  OF  ELEMENTS 

663 

500 

X  =  XL  SY=0.  SSCALE=1. 

6b4 

NQL=N8AYLI+1 

6b5 

N8h*N8  A Y  hi *1 

666 

00  550  1=1, NBL 

667 

00  540  K=1,N8M 

668 

CALL  NLM0ER(X-.O7,Y  +  O.,.1O,SCAL'E,O.  ,4HF4.0) 

6o9 

SCALE  =SCALE  +1. 

670 

Y=Y*ONEH 

671 

U1 

4? 

C* 

CONTINLE 

672 

Y  =  0.  ' 

673 

X=X-0NEL 

674 

550 

CCNTINLE 

675 

X  =XL  $NUM8=i$Y  =  0« 

676 

00  2003  K= 1* N8 AYLI 

677 

OG  200C  I=1,NBAYWI 

678 

CALL  NUM8ER(X, YV,4*ONCH|.07,NUMB,9J. ,2HI3> 

679 

Y  =Y  +  ON£W 

680 

NLMB  =  NLM8-*-l 

681 

2000 

CONTINLE  • 

682 

2001 

Y=3. 

683 

OG  2002  J=1,NBAYHI 

684 

CALL  NLMBER(X-.5*ONEL, Y- . 05*ONEW , . 07 ,NUHB ,90 

685 

NUN8  =  NLM»>1 

666 

CALL  NLMBER(X-.25*ONEl,Y».25*ONEW,.07,NUHB,1 

687 

NUMB=NUMB+1 

688 

CALL  NUMBER (X-. 375*ONtL, Y* . 625*0NEH , .0 7, NUME 

689 

Y=Y+ONEH 

690 

NUM0=NUM8+1 

691 

2902 

CCNTINLE 

692 

CALL  NLMBERt X-.5*0N£L, Y-.Q5*ON£W,.Q7,NUM6,9Q 

693 

X=X-ON£L 

694 

Y=O.JNLMB=NUMBU 

695 

2003 

CONTINLE 

696 

Y  =0 .  J X  =  0  • 

697 

00  2004  1= 1,NB AYWI 

698 

CALL  NL’ri8ER(X  ,  Y+ •  4*0N£H ,  •  0  7.NUMB ,  90  •  , 2HI3) 

699 

Y=Y+ONEW 

700 

NUM8=NLM8+1 

701 

2J04 

CCNTINLE 

80 


m 

703 

704 

705 

706 

707 

708 

709 

710 

711 

712 

713 

714 

715 

716 

717 

718 

719 

720 

721 

722 

723 

724 

725 

726 

727 

728 

729 

730 

731 

732 

733 

734 

735 

736 

737 

738 

739 

740 

741 

742 

743 

744 

745 

746 

747 

748 

749 

750 

751 

752 

753 

754 

755 

756 

757 

758 

759 

760 

761 

762 

763 

764 

765 


GO  TO  aoo 
800  CONTINUE 

XP0S=(XSCAL£/2.)-.5 
CAUL  FL0T(XP0S,-.35,3) 

CALL  PLQT(XPQS,-.35,2) 

CALL  SYN80L(XPCS»-.35,.12» 11HGRI0  SYSTEM, 0. ,11) 

CALL  SYM8OL(XPCS,".50» .10, 7HPR0BLEM, 0 » , 7) 

CALL  NUMBER (XPCS+. 65 50 , .10 ,NPROB, 0. ,2HI5> 

CALL  PLOT (XSCALE+4.  ,  0. »”3) 

6000  RETURN 

ENO 

SUBROUTINE  ARCHPLT (RK,X , Z, D»M ANGLE, PLOTN , IB) 

C0NM0N/3/CA. JERK, NBA YL, NBA YW,NUMPR,NPROB 
CALL  MAP  (  *1  •  ,  20  .  , “1  •  ,  2C  «  ,  0  «  ,1  •  ,  0  •  ,  1  • ) 

C  RK  =  1.0  /  PCENT  OF  LOAD 

C  X=X  FROM  MAIN  PROGRAM. 

C  Z=Z  FRCM  MAIN  PROGRAM. 

C  0=DISPLACEMENT  • 

C  MANGLE=FULL  OR  HALF  ARCH. 

DIMENSION  X ( 1 > ,2(1) ,0(1) 

0 IME  NS  ION  XXQ(21) , ZZ 0 (21) , XOR < 21 ) ,Z0R(21) ,ZCCR(21) 

C  NUNX  =  NUMBER  CP  NOOES  TO  PLOT  ARCH 

C  NUMZ  =  NUMBER  OF  NOOES  TO  PLOT  CROWN  DEFLECTIONS 

READ  1,NUMX,NUKZ 
1  FORMAT (2110) 

DO  10  !=1,NUMX 

C  NOO  *,  NODE  NUMBERS  TO  BE  USEO  TO  PLOT  THE  ARCH 
READ  1 »NCO 
MOO=5*NOO-4 
MOUl=5*NOO-2 

C  '  CALCULATION  OF  COORDINATES  OF  CENTERLINE  ARCH.  AKO  DEFLECTED  CENTERLINE 
C  ARCH. 

C  XXU=OEFLECTEQ  COORDINATES. 

C  ZZO=UEFLECTEO  COORDINATES. 

C  XOR=ORIGINAL  CCOROINATES. 

C  ZOR=ORIGINAL  COORDINATES. 

XXD(I)=(0(M00) *RK) •10.+X (NOO) 

ZZO  ( I)  =  (0  (MODI)  *RK)  *10  . +-Z(NOO) 

XOR (I)=X (NOO) 

ZOR ( I ) =  Z (NOO)  • 

10  CONTINUE 

C  MANGLE  =  1  OR  2  DEPENDING  ON -WHETHER  IT  IS  A  FULL  ARCH  OR  HALF  AACP 

IF (MANGLE. EQ. 2 ) GO  TO  10G 
DO  11  1= 1, NUNX 

C  SCALING  OF  COORDINATES. 

XXO(I)=XXO(I>/ (-50.)  . 

ZZO (I) =ZZD (I ) /50 . 

XOR(I)sXOR(I)/(-50.) 

ZOR (I) =ZOR  (I>/50. 

11  CONTINUE 

C  PLOTTING  AXES. 

YY=10. 

CALL  PLOT (0. ,YY,3) 

CALL  PLOT (0. ,YY,2) 

00  12  IT=1,10 
CALL  PLOT ( 0. , YY, 1) 

CALL  PLOT(0. ,YY-1, ,1) 

CALL  PLOT (*.1,YY”1,,1) 

CALL  PLOT (0. , Y Y-l • , 1 ) 

12  YY=YY-1. 

xx-o. 

DO  13  IX-1,16 
CALL  PLOT( XX , 0 , , 1) 

CALL  PLOT (XX+1 .,0. ,1) 
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766 

767 
766 

769 

770 

771 

772 

773 

774 

775 

776 

777 

778 

779 
7d0 
7ai 

782 

783 
764 

785 

786 

787 
7d8 

789 

790 

791 

792 

793 

794 

795 
79b 
797 

796 

799 

800 
801 
802 

803 

804 

805 

806 
807 
608 

809 

810 
811 
812 

813 

814 

815 

816 

817 

818 
819 
8*0 
821 
822 

823 

824 

825 

826 

827 

828 


CALL  PLOT(XX*l.,-.l,l) 

CALL  PIOT(XX+1.,0.,1) 

13  XX=XX+i. 

C  NUMBERING  OF  AXES. 

SCAL1NC=50.  IYY=10.  JYSCALE=10.*SCALlNC 
OC  14  IY-1  ,11 

CALL  NUMBER <-l.tYYi.lU»YSCALE,0. ,4HF8.2) 

YY=YY-1. 

14  YSCALE= YSCALE-SCAL1NC 
XSCAl1*50<  IXX=0.  SXSCA=0. 

00  15  1X=1,17 

CALL  NLMBER<XX-.25,-.25,.10,XSCA,0. ,4HF8.2) 

15  XSCA=XSCA+XSCAL1 
CALL  PLOT ( 0 . *  0  •  »3) 

CALL  PLOT (0  . ,  0. ,  2) 

00  16  1=1, NUMX 

C  PLOTTING  DEFLECTEO  ARCH. 

CALL  SYMBOL <XXC<I) , ZZO <11 , . 07, 1,0. ,-l) 

16  CONTINUE 

CC  17  I=i,NUMX 

C  PLOTTING  ORIGINAL  ARCH. 

CALL  SYMBOL  ( XOM I )  , ZOR  (I ) , . 07 , 2, C . , -1 ) 

17  CONTINUE 
IF(PLOTN)200,200,3G3 

200  CALL  SYMBOL  (6  .  ,  -  .51)  ,  .  12, 49H0EFLECTI0NS-L INEAR  SOLUTION  (MAGN1FICA 
1TICN  OF  10 ) , 0  . «48) 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXkXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
C  PLOT  NUMBER  7 

C  CENTER  ARCH  ANC  DEFLECTED  CENTER  ARCH (LINEAR  SOLUTION) 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
GO  TO  400 

300  call  SYMBOL (6  <  •‘*’•50  i  « 12* 52HCEFLECTI0NS-NCNL INEAR  SOLUTION  (NAGN1F1 
1CATION  OF  10) ,0.  ,51) 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX»XXXXX 
C  PLOT  NUMBER  9 

C  CENTER  ARCH  ANC  OEFLECTEO  CENTER  ARCH (NONLINEAR  SOLUTION) 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXAXXXXXXXXXXXXXXXXXXXXXXXXX 
400  CONTINUE 

CALL  SYMBOL (7 • 5 65 « • 10 , 7HPROBLEM , 0 . ,7) 

CALL  NUMBER (8. 5,-<65,.lG,NPR0B,Q.,2HI5) 

CALL  PLOT (20 • , 0. , “3) 

IF  ( IB. EU .1 ) 22 , 25 
100  00  26  1=1, NUMX 

XXO ( I ) =XXO (1) /50  . 

ZZO ( I) =ZZO (I) /50. 

XOR(I) =XOR(I) /50.  ( 

ZOR ( I ) =Z0R ( I ) 750  • 

26  CONTINUE 

C  PLOTTING  AXES. 

YY  =  10.  IX X=1 6  • 

CALL  PLOT  <  XX , YY, 3) 

CALL  PL0T(XX,YY,2) 

00  27  IT=1 , 10 
CALL  PL0T(XX»YY,1) 

CALL  F10T(XX, YY-1.,1) 

CALL  PL0T(XX-.l, YY-1.,1) 

CALL  PLOTtXX, YY-1.,1) 

27  YY=YY-1. 


* 


I 


830 

831 

832 

833 

834 

835 

836 
83  7 

838 

839 

840 

841 

842 

843 

844 

845 

846 

847 

848 

849 

850 

851 

852 

853 

854 

855 

856 
657 

858 

859 

860 
861 
862 

863 

864 

865 

866 

867 

868 
8t>9 
d/0 

871 

872 

873 

874 

875 

876 

877 

878 

879 

880 
881 
882 

883 

884 

885 

886 

887 

888 
8o9 

890 

891 

892 

893 


UO  28  IX*1,32 
CAU.  PLOT(XX,0.,1) 

CALL  PLOT(XX>1.*0. »1» 

CALL  PLOT (XX41 • *“• 1,  1) 

CALL  PLCT(XX41.,0.»1) 

28  XX*XX+i. 

C  NUMBERING  OF  AXES. 

YS=50.  m=10.  SYSCA=10.*YS 
00  29  IY=l.li 

CALL  NLHBER<7. ,YY,.1Q,YSCA,0.,4HF8.2) 

YY=YY-1. 

29  YSCAsYSCA-YS 

CALL  SYM80L(-1.5,4.,.iO,iHZ,9Q.»l> 

XSCAIN=50«  SXX*0.  $XSCA=-400. 

00  30  IX=i *17 

CALL  NUrtBER<XX-.25,-.25,.iO,XSCA,0. ,4HF8.2> 

XX=XX41. 

30  XSCA=XSC A4XSCAIN 

CALL  SVN8OL<8.,-l.,.i0,lHX,Q.,i) 

XX=8.SYY*0. 

CALL  PLOT <XX»Y7»3) 

CALL  PLOT <XX,YY»2) 

00  31  I=1,NUMX 

C  PLOTTING  OEFLECTEO  ARCH. 

CALL  SVH0OL<XXG<I)4XX,ZZO<I)',.O7#1,O.,-1> 

31  CONTINUE 

DO  32  1*1 »NUMX 

C  PLOTTING  ORIGINAL  ARCH. 

CALL  SYH8GI.  ( XOK ( I ) +XX ,  ZOR(  I  > »  •  07  .2.0. .-I) 

32  CONTINUE 
IF(PLOTN)201, 201,301 

201  'CALL  SYH80L  ( 1 2. , - .50 , . 12 , 49H0EFLECTI0NS-L INfcAR  SOLUTION  (HAGNIFIC 
1ATION  OF  10) , 0 . ,48) 

GO  TO  401 

301  CALL  SYMBOL  (12., -*56, .12, 52H0EFLECT IONS 'NONLINEAR  SOLUTION  (MAGNI 
1FICATICN  OF  10 ) ,0 . ,51) 

401  CONTINUE 

CALL  SYMBOL (15. ,-.65, . 10 , 7HPROBLEM, 0 . ,7) 

CALL  NIMBER(16.,-.65,.10,NPROB,0.,2HI5> 

CALL  SYM80L(9. ,5. ,.12, 3HRK= , 0. , 3) 

CALL  NLMBER(9.5,5.,.1’2,RK,0.  ,4HF5 .2 ) 

CALL  PLOT (40 • , 0. * — 3 ) 

IF(ia.EQ.l)GO  TO  22 
25  00  18  I=1,NUMZ 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
C  PLOT  NUMBER  8 

C  CROHN  CEFLECTIONS 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
Cr.XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
C  NOO=NOOE  NUMBERS  TO  BE  USEO  TO  CALCULATE  CROHN  DEFLECTIONS. 

REAOl.NOO 

M0Cl=5«N00-2 

C  CALCULATION  OF  CROHN  DEFLECTIONS. 

C  ZC0R=CR0HN  DEFLECTIONS. 

ZC0R(I)*(D(M0Q1)*RK) 

16  CONTINUE 
C  ORAHING  AXES. 

XX=0.  SYY=5. 

CALL  PLOT ( XX , YY, 3) 

CALL  PLCT(XX»YY»2) 

TZ=8./NUMZ 
OC  191=1, NUMZ 
CALL  PLOT (XX»YY»1) 
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894  CALL  PL0T(XX+TZ,YY,1) 

895  CALL  PLOT  (XX+TZ*YY-«1,1) 

896  CALL  FLOT(XX+TZ,YY,l) 

897  19  XX=XX+1Z 

898  XX=0.  JYY=5. 

899  CALL  PLOT(XX,YY,3)  > 

9u0  00  21  Isl.NUKZ 

901  C  PLOTTING  OF  CROHN  OEFlECTIONS. 

902  CALL  SYHQOL(XX,ZCOK(I)+YV,.Q7,2,0.,-1> 

903  21  XX=XX+1Z 

904  CALL  SYMB0L(2.S,-.35,.12,28HCR0WN  DISPLACEMENTS  (INCHES) ,0., 28) 

905  CALL  SYMB0L(2.,-.75,.12,32HVERTICAL  SCALE  -  1  INCH  =  1  INCH, 0. ,32) 

906  CALL  SYMBOL (3«,~,5,,1Q » 7HPROBLEM , 0  #  ,7) 

907  CALL  NLMBER(4.,-.5,.1U,NPROB,0.»2HI5) 

908  CALL  PLOT(12.,0.,-3) 

909  22  CGNTINLE 

910  RETURN 

911  ENO 

912  SUBROUTINE  STRESSP 

913  CALL  MAP{-i.,ll.,-l.,il.,6.,l.,0.,l.) 

914  COMMON/l /A (1000) ,8(1000) ,C(1000) ,0(1000) ,E (10) , 

915  1  F(1000) ,GC1000) , 

916  2  H(220) ,0(220) ,P (220) , S (231) ,U(23l) ,V(ll55)  ,H(1155) 

917  OIMENSION  X  ( 1  00)  ,  Y  (100 )  ,  T  ( 100)  ,R  (10  0  )  ,  R1  (IOC  )  ,  Xj.  (1000)  ,Y1(10JC)  ,X2 

918  1(1030)  ,Y2(1000) 

919  C  A-STRSCJ 

920  -C  6-STRS1J 

921  C  C-STRSCK 

922  C  O-STRSIK 

923  C  t (1)=CA,£(2) = JERK ,£(3)=NBAYL,E (4)=NBAYH,E (5)=NUMPR,£(6) =NPR08 

924  I=(E(4)/2.)*E(3)+1. 

925  NBAYL=E ( 3)  2  NBAYH=E(4)  $NUMPR=E(5)  $NPRCB*E(b)  SCA=E(1) 

926  NUH=I +NOAYL 

927  NUMT=NEAYL 

92o  J=NDAYL+ 1 

929  PRINT1 ,1 ,NUM,NLMT,  J 

930  1  FORHA1 (4 (5X,I10) ) 

9ol  F  (1) =A  (I )  JF(J)=C(NUH) 

932  G  (1 ) =B ( I )  SO ( J) =0 (NUM) 

933  00  10  K=2,NBAYL 

934  F  (K) - ( C 1 1 ) *A ( I +1 ) )/2« 

935  G  (K)  =  (C (I) +8(1+1) )/2. 

936  PRINT300  0,F(K),G(K>,C(I>,A(IU),0(I>,B(I*1) 

937  300o  F GREAT (5X,5HF ( K)  =  ,F1 0 . 4 , 5X , 5HG ( K) =, F 1 0 , 4 ,5X , 5HC ( I ) *  ,Fi 0 .4,5X ,7P A  ( I 

938  l+l)=,F10.4,SX,5HO(l)=,FlQ.4,5X,7HB(I+i>=,F10.4) 

939  1=1+1 

940  10  CGNTINLE 

9+1  PRINT  3001, CA 

942  3001  FORMAT (5 X , 3HCA= «F1G« 4) 

943  Y (1) =0  . 

9+4  X (1) =8  • 

945  00  12  K  =  2 , 81 

946  Y (K) =Y (K-1J+.1 

947  X(K)=SCRTF<64-Y(K)+Y(K>) 

948  12  CCNTINLE 

949  PI=3. 14159265 

950  8ETA=(Fl-CA)/2, 

951  L=0 

952  READ  2 ,  SF 

953  2  FORMAT (F10 .4) 

954  00  13  K=l, J 

955  T <K) = ( (Pl-CA) /2.)+L* (CA/20, )  ' 

956  R(K) =8.+ (SF*F(K) ) 

957  R1(K)=8.+(SF*G(K)) 


84 


Rft . 

. 

~ .  .  • . 

fj 

1 

958 

PRINT  2002, TOO, RU>  ,*1(K> 

1 

959 

3002  F CRM AT  <5X»5HT(IO  =  »F1Q.4*5X,5HROO  =  *F10.4»5X>  6HR1  (KT=,F10.4) 

1 

960 

L=Ltl 

£  v 

9ol 

13  COM  INLE 

962 

DO  14  K=1,J 

•« 

963 

TAK=T(K) 

§ 

964 

Xl(K)=fi(K)*  COS(TAK) 

965 

Y1 (K ) =R ( K) *  SIMTAK) 

966 

X2(K)=R1<K)*  CCS(TAK) 

* 

r 

967 

Y2<K>=fii <K)*  SIN (TAK> 

9 

968 

14  CCNTINLE 

l 

969 

2EROX=fi.*  COS  CT(1>  > 

& 

970 

2EROY=8.*  SIN(T(1>) 

1 

971 

LP=Q 

? 

fc 

972 

.1000  YY  =  8.  |LP=LP+1 

| 

973 

CALL  PLOT(0.,YY,3) 

1 

974 

CALL  PLOT (d. , YY,2) 

I 

975 

OO  15  Ml, 8 

k 

976 

CALL  PLOT(0.,YY-1.,1) 

* 

977 

CALC  PLOT(-.l,YY-l.,l> 

l 

978 

CALL  PLOT (0. , YY“1.»1> 

\ 

979 

YY=YY-1. 

> 

96(i 

15  CCNTINLE 

1’ 

981 

xx=o. 

t 

982 

OO  16  Ml, 9 

s. 

983 

CALL  PLOT<XX»1.,0.,1> 

1 

9a4 

CALL  PLOT«XXtl.,-«l»l) 

s 

;  \ 

965 

CALL  PLOT  <  XX*1  • , 0 • «1> 

1 

9o6 

•XX=XX»1. 

t# 

967 

16  CONTINUE 

1 

986 

YY=8.  1  XX=-.5 

l 

939 

NLH  =  8 

*  * 

99u 

OO  5C0 0  1-1,9 

\ 

991 

CALL  NLM8ER(XX,YY,.1Q,NUH,0.,2HI3> 

t 

592 

NUM=NUF-1 

\ 

993 

5000  YY=YY-1. 

X 

i~ 

994 

NUMX=0  %  YY=-.3  $  XX=J 

Sr 

995 

OO  5001  1=1,10 

996 

CALC  KLMBER(XX ,YY , .10 , NUMX , 0 • ,2H13) 

£' 

s 

997 

NUMX=NLMX+1 

998 

■  5001  XX=XX4-1 

S»r 

999 

CALL  PLOT<3.,0.,3) 

r 

1000 

CALL  PLOT ( 0.  ,  0  • ,  2) 

i, 

1001 

OO  17  K= 1 , 81 

L 

1C02 

call  plot <xoo  ,'y(ki,i> 

i 

1003 

17  CONTINLt 

r. 

1004 

CALL  PLOT (  0.  ,  0  .,  3) 

1005 

CALL  PLOT<0. ,0.,2) 

s 

1006 

CALL  PLOT <ZEROX, ZERO Y, 1) 

1007 

CALL  PLOT ( 0. ,  0  • ,  3) 

! 

1008 

CALL  PLOT(0.,0.,2) 

1009 

UNIT=1,/SF 

* 

!> 

1010 

IF (LP.hQ  *2) GO  TO  2000 

? 

fc. 

1U11 

GXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

f 

1012 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

> 

1013 

C  PLOT  NLM8ER  10 

5 

1014 

C  AVERAGE  OF  J-ENO  ANO  K-tNO  OUTSIDE  STRESSES  ALONG  CENTER  ARCH 

1 

1015 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

1016 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

[? 

1017 

OO  18  K*1,J 

- 

g  f 

1016 

CALL  S YM80L< XI <K) ,Yl<Kf ,.07,2,0. ,-l) 

l  If 

l  I  • 

1019 

18  CONTINLE 

1020 

CALL  SYM80L<3.,-.50, .12,16H0UTS1DE  STRESSES, 0. ,i6) 

i  J 

1021 

CALL  SYMBOL (3. 2, -.70 , , 10 , 7HPROBLEM , 0 . , 7 > 

s  f 

1 

85 

- - - - - - - -  ~~  ”  ~ 

1C22 

1023 

1024 

1025 

1026 
1027 
1026 
1025 

1030 

1031 

1032 

1033 

1034 

1035 
10o6 
1037 
1030 

1039 

1040 

1041 

1042 
10*3 

1044 

1045 


CALL  NLN8£R<4„0,-.70  ,  . 10 , NPR08.0 . ,2HI5) 

CALL  SYM80L(J.,-.9,.lu,7H  UNIT  =,0.,fl) 

CALL  NUMBER <3. 7, -.9, . 10 , UNIT ,0 . , 4HF5 .2 > 

CALL  SYMBOL <4.0,-. 9*. 10, 3HKSI,0. ,3) 

CALC  PLOT ( 12. , G. ,-3) 

GC  TO  1000 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
C  PLOT  NUMBER  11 

C  AVERAGE  OF  J-ENO  ANO  K-END  INSIOE  STRESSES  ALONG  CENTER  ARCH 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXU 
200c  OO  19  K=l»  J 

CALL  SYM80L(X2(K),Y2(i(), .07,2,0,,-!) 

19  CONTINLE 

CALL  SYM80L(3.,-.5U,.12,15HINSI0E  STRESSES, 0., 15)  • 

CALL  SYMBOL <3. 2 ,-.7, .10 , 7H PROBLEM, 0 . ,7) 

CALL  NLM8£R(4.0,-«7,.10, NPROB, 0.,2HI5) 

CALL  SYMBOL<3.,-«9,.10,7H  UNIT  =,0.,6) 

CALL  NLM6ER(3.7,-.9;.10,UNIT,0.,4HF5.2> 

CALL  SYMBOL (4.0, -.9, .iO,3HKSI,Q. ,3) 

CALL  PLOT ( 12. • G. , -3) 

RETURN 

END 
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APPENDIX  II 


INPUT  INSTRUCTIONS  FOR  COMPUTER  PROGRAMS 


1.  INTRODUCTION 

The  basic  program,  SANOS,  can  be  executed  independently  but  will  yield 
only  printed  results;  to  obtain  plots  of  the  results,  the  SAN1PLT  program 
must  be  executed. 


2.  DATA  FOR  SANOS 

Format  _ Information _ 

I  10  (1)  NUMPR  -  number  of  problems  in  this  computer  run 

8110  (2)  NPROB,  NCODE,  NRAYL,  NBAYW,  NLOAD,  NSUPN,  NRK,  MANGLE 

where:  NPROB  =  problem  identification  number 

NCODE  =  0,  nodes  numbered  in  direction  of 
shell  axis 

=  1,  nodes  numbered  along  arch 
NBAYL  =  number  of  bays  along  shell  axis 
NBAYW  =  number  of  bays  along  arch 
NLOAD  =  number  of  loads  on  shell 
NSUPN  =  number  of  nodes  which  have  supports 
NRK  *  number  of  Rungc-Kutta  (loading) 
intervals 

MANGLE  =  1,  one-half  arch  (symmetrical  load¬ 
ing) 

-  2,  full  arch  (unsymmetrical  loading) 

3F10.0  (3)  SPAN,  HEIGHT,  SENGTH 

where:  SPAN  =  span  of  shelter  (total)  (inches) 

HEIGHT  =  height  of  shelter  (inches) 

SENGTH  =  length  of  shelter  (inches) 

7F10.0  (4)  *THICKM,  THICKW,  ULTM,  THICKA,  ULTSTR,  E,  R 

where:  THICKM  =  equivalent  thickness  of  shell¬ 
bending-strong  axis  (inches) 

THICKW  =  equivalent  thickness  of  shell¬ 
bending-weak  axis  (inches) 

UL1M  =  ultimate  moment  (inches -kip/inch) 
THICKA  =  equivalent  thickness  of  shell  - 
axial  (inches) 

ULTSTR  =  ultimate  axial  stress  (ksi/inch) 

E  =  modulus  of  elasticity  (ksi) 

R  =  curve-fitting  parameter 


This  information  is  obtained  from  laboratory  tests. 
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Format 


Information 


5. 


3F10.0 


8F10.0 


(5)  COUT,  CIN,  RATNDD 

where:  COUT  =  distance  from  neutral  axis  to  outside 
fibers  (inches) 

CIN  =  distance  from  neutral  axis  to  inside 
fibers  (inches) 

RATMOD  =  E  strong/E  weak  (ratio  of  moduli) 

(6)  PCENT  (I),  I  =  1,  NRK 

percent  of  load  (in  decimal)  on  each  increment  of 
loading 


6110 


110,  F10.0 
313,  Ill 


(7)  NODE,  NP(J),  J  =  1,5 

support  information  0  indicates  no  support 

1  indicates  supported  against  displacement  in  x 
direction 

2  indicates  supported  against  displacement  in  y 
direction 

3  indicates  supported  against  displacement  in  z 
direction 

4  indicates  supported  against  rotation  about  x  axis 

5  indicates  supported  against  rotation  about  y  axis 

(8)  I,  D(I)  NLOAD  times 

loading  conditions  (I  indicates  coordinate  number) 

(9)  MN,  JM,  KM,  MTYPE 

where:  NN  =  member  number 
JM  =  j  end  of  member 
KM  =  k  end  of  member 
MTYPE  =  1,  interior  member  along  arch 

2,  exterior  member  along  arch 

3,  interior  horizontal  member 

4,  exterior  horizontal  member 

5,  diagonal  member 


F10.0  (10)  CONTROL 

where:  CONTROL  =  0  --  no  plots  required 
=  value  --  plots  required 


DATA  FOR  SAN1PLT 

Format 

Information 

F10.0 

(1) 

PLOTCON 

If  PLOTCON  =  0,  signals  last  set  of  plots 
If  PLOTCON  =  value,  more  plots  to  come 

,110 

(2) 

L  (Determines  type  of  loads) 

L  =  1  --  loads  in  x  direction 

L  =  2  --  loads  in  y  direction 

L  =  3  --  loads  in  z  direction 

L  =  4  --  moment  about  x  axis 

L  =  5  --  moment  about  y  axis 
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Format 


Information 


2110 

(3) 

NUMX,  NUMZ 

NUMX  =  number  of  nodes  around  arch  desired  for 
plotting 

NUMZ  =  number  of  nodes  along  crown  for  plotting 

no 

(4) 

NOD  (node  numbers  for  arch  deflection  plot  of  linear 
solution,  NUMX  of  them) 

no 

(5) 

NOD  (node  numbers  for  crown  deflection  plot  of  linear 
solution,  NUMZ  of  them) 

no 

(6) 

NOD  (node  numbers  for  arch  deflection  plot  of  non¬ 
linear  solution,  NUMX  of  them) 

F10.4 

(7) 

SF  (scale  factor  desired  for  plots  of  inside  and  out¬ 
side  stresses) 

APPENDIX  III 


f 

? 

s 


r 


I- 


I 


SAMPLE  INPUT  DATA  FOR  SANOS 


Test  problem  5  is 

presented  since 

it  is  an 

example  of  unsymmetric  loading 

5 

0 

10 

20 

lo 

23  6  2 

590. 

290. 

360. 

2.4$ 

0.15 

15.0 

0.021 

81.5 

300 O0 .  7.0 

2.00 

17.0 

7.5 

1  V 

.166? 

.1667 

.1667 

.1667 

.1667 

■.  1667 

1 

1 

2 

3 

4 

2 

1 

2 

3 

4 

3. 

1 

2 

3 

4 

4 

1 

2 

3 

4 

4 

1 

7 

3 

4 

6 

1 

2 

3 

4 

7 

1 

2 

3 

4 

8 

1 

7 

3 

4 

0 

1 

2 

3 

4 

1* 

J. 

2 

3 

4 

n 

1 

2 

3 

4 

??i 

1 

2 

3 

4 

277 

1 

2 

3 

4 

273 

1 

2 

3 

4 

‘ 

224 

1 

2 

3 

4 

2?5 

1 

2 

3 

4 

276 

1 

2 

3 

4 

277 

1 

0 

1. 

3 

4 

??fl 

1 

2 

3 

4 

270 

1 

7 

3 

4 

730 

1 

2 

3 

4 

231 

1 

2 

3 

4 

13  6 

0 

2 

0 

0 

&.S  P 

-2.00 

463 

-2.00 

473 

-2.00 

403 

-2.00 

503 

-1.52 

51« 

-1.57 

570 

-1.5? 

5  0.0 

-1.52 

53  3 

-2  .OQ 

550 

-4.17 

563 

-4,17 

568 

-4.17 

573 

-4.17 

57P 

-4,17 

603 

-4  •  1  7 

500 

-4.17 

503 

-4.17 

5o  p 

-4.17 

6O0 

-2. Oo 

90 


1  301002 

2  302  303 . 

3  303  304 
4304305 
5305306 

6  306  307 

7  307  308 
8308309 
9309310 

10310311 
1100lbl2 
12301(013 
13  302  012. 
14302013 
15002014 
16003013 
17003014 
18303015 
19004314 
20004315 
21004016 
22005015 
23005016 
24005017 
25006016 
26006)017 
27006018 
28  00701.7 
29007018 
30007019 
31008018 
32008019 
33008020 
34009019 
35009020 
36009021 
37010020 
38010021 
39010022 
40011021 
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APPENDIX  IV 


SAMPLE  OUTPUT  DATA  FOR  SANOS 


Test  problem  5  is  presented  since  it  is  an  example  of  unsymmetrical  loud- 


PRORLFM  NUMREO  5 

NONLINEAR  ANALYSIS  OF  ORTHOTROPIC  SHELLS 
HATERIAL  NONLINEARITY  »•  RICHARO  FORMULATION 
PLASTIC  ANALYSTS  *•  SMITH  PLASTIC  MODEL 
0E0MFTOIC  NONLINEAP.ITY  ••  INCREMENTAL  TECHNIQUE 


METHOD  OF  NUMBERING  NODES  »*•••••»•*•• 
NUHAFP  OF  SAYS  LONG  »•*••***•••»•••»»• 
NUMBER  OF  PAYS  NIOE  »♦*»»♦*»♦*♦***»»»* 
NUMBER  OF  LOADS  *•*•••*••*•*••*•**•**• 
NUMBER  OF  NODES  WITH  SUPPORT'S  *••••♦»• 
NUMBER  OF  RUNGE-KUTTA  INTERVALS  •*•*»• 
AMOUNT  OF  STRUCTURE  MnCFtpo 


SPAN  OF  STRUCTURE  »• 
hfight  of  rtructuke 
length  of  SMELL  •••• 


THICKNESS  OF  SHELL  IRFnoTNG) 

COUBLF  THICKNESS  of  METAL (TMICKH  CAL)* 
ULTIMATE  MOMENT  (IN-KtPS/INI  ••♦••*•♦• 
THIOKNFSS  OF  SHELL  < AXI AL )•***•**»*•* * 
ULT,  AXIAL  STRESS  (KIPS/SO.  IN/INCH)  •• 
MODULUS  OF  PLASTICITY  OF  SHrLL  ******* 
CURVE  FITTING  PAPAMETFP  •**♦.*.**»*.*• 


DISTANCE  FROM  N.A,  TO  OUTSIDF  FIBERS  * 
OTSTANCE  FROH  N.A.  TO  TNSIOr  PIPERS  »» 
modular  »atio(Uncobr,  to  rcRR.)  ****** 


NCOIE 

s 

0 

NRAYL 

= 

10 

NRAYW 

s 

20 

NLOAO 

s 

1R 

NSUPN 

c 

21 

NIK 

s 
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H  ANGLE 

s 

2 

SPAN 

r 

500.00 

HEIGHT 
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205.00 

SENGTIt 
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150,00 

THir.<M 

s 

2.5H 

THICK 

s 

.15 

ULIH 

s 

15.  OU 

THIC<A 

s 

.02 

ULTSIP 

S 

At.  50 

t 

5 

30000.00 

R 

= 

2.00 

COUT 

_ 

2.00 

Cl  N 

= 

12.00 

"ATHOO 

= 

7.50 

NUMBER  OF  MEMBERS  •*«*•«*•••**•••»*♦»• 
NUMBER  OF  NOOES  »•**♦•**#«*»•»•*»»*»** 
NUMRFP  OF  FOUATIONS 

BAND  MOTH  ••»***••••*•**»•••*•»•♦»»»* 
HALF  RANO  HtnTH  PLUS  OIACOHAL  **»•••»• 
NUMREP  OF  SUPPORTS  ••*»»»••»*»♦•*••»•» 
SMALLER  NOOE  NliNPFR  ON  ENO  OF  MFMDFR  * 
LARGER  NOOE  NUMBER  ON  fro  of  MEHRFR  *» 
THICKNESS  OF  SHELL  (RFNOINC-HEAK  AXIS) 
CENTRAL  angle- 'DECREES* *•**•»** •»»•»*• 

length  oe  longitudinal  MrHREPS  •***••• 

RADIUS  OE  SHF l L  **•••»•••*•••••*•»•••• 


■NM 

= 

830 
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231 

NO)5 
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11*55 

NTA'JO 
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1?9 

NHS  /M 
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65 

NSUPS 

= 

89 
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13 
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46.291 

59.642 

53.642 
46.291 

58.642 

53.642 
46.291 

58.642 
58.642 
46.291 
58.642 
58.642 
46.291 
58.642 

58.642 
46.291 

53.642 

58.642 
46.291 
58.642 

58.642 
46.291 

53.642 
53.642 
46.291 
36. 000 
36.09U 
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700  140  141 

.472 

70,227 

*.1*0 

694.363 

36. 000 

3 

701  141  142 

.472 

79.227 

0.1*0 

644.363 

36.000 

3 

702  142  143 

.972 

74.227 

*.1*0 

694.363 

36.900 

3 

703  143  144 

.972 

79.227 

*.1*0 

644.363 

36.000 

3 

704  144  145 

.972 

79.227 

8.1*0 

694.363 

36.000 

3 

705  145  146 

.472 

79.227 

8. 180 

694,363 

36.  000 

3 

706  146  147 

.972 

79.227 

*,180 

694.363 

36.000 

3 

707  197  190 

.472 

79.227 

8.1*0 

644.363 

36.000 

i 

700  18*  149 

.37* 

30. *07 

22.179 

270.000 

46. 2°1 

? 

709  18*  200 

.010 

10.000 

.010 

10,000 

51.642 

710  189  199 

.010 

10,000 

.010 

10.000 

58.642 

5 

711  1*9  200 

.756 

61.614 

45.759 

540,000 

46.241 

t 

712  1*9  201 

.010 

10.000 

.010 

10.001 

58.642 

5 

713  140  ?00 

.010 

10.000 

.010 

10.000 

51.642 

b 

714  140  201 

.  756 

61.614 

45.759 

540.000 

46.291 

1 

715  190  ’0’ 

.010 

10.000 

.010 

in, 000 

58.642 

5 

716  191  201 

.010 

10.000 

.010 

10.000 

51.642 

* j 

717  191  202 

.756 

61.614 

45.759 

540.000 

46. 291 

1 

VI*  141  203 

.010 

10.000 

.010 

10.000 

51.642 

5 

719  192  »9? 

.010 

10.000 

.110 

10.000 

51.642 

5 

720  192  203 

.756 

61.614 

45.759 

540.000 

46.291 

l 

721  192  204 

.'010 

10.000 

.010 

10.000 

51.642 

5 

722  143  203 

.010 

10.000 

.110 

10.  000 

51.642 

V 

723  143  ’04 

.756 

61.614 

45.759 

541.000 

46.291 

1 

724  193  205 

.010 

10.000 

.110 

10.000 

58.642 

5 

725  194  ’04 

.010 

IO.OCO 

.010 

10.000 

58.642 

5 

726  194  205 

.756 

61  .614 

45.759 

540.000 

46.291 

t 

727  194  206 

.010 

10.000 

.010 

10.000 

51.642 

5 

728  195  205 

.010 

10.000 

.010 

10.000 

51.642 

724  195  206 

.756 

61.614 

45.759 

540.000 

46.291 

1 

730  195  207 

.010 

lo.oao 

.110 

10.000 

51.642 

5 

731  196  706 

.010 

10.000 

.010 

10.000 

58.642 

5 

7?2  196  707 

.756 

61.614 

45.75a 

540.000 

46.291 

1 

737  196  200 

.010 

10.000 

.010 

in.000 

51.642 

5» 

734  147  207 

.010 

10  .000 

.010 

10.000 

51.642 

5 

735  107  210 

.756 

61.614 

45.754 

540.000 

46. ’91 

1 

736  197  209 

.010 

10.000 

.010 

10.000 

51.642 

5 

737  19*  200  ■ 

.010 

10.000 

.110 

10.00U 

51.642 

5 

738  19*  704 

.378 

30. *07 

22.179 

270.000 

46.291 

? 

7’9  199  200 

.97? 

79.227 

*.1*0 

694.363 

36.000 

3 

740  200  201 

.97? 

74.227 

*.t*0 

694.363 

36.000 

3 

741  201  20’ 

.«72 

74.227 

8.1*0 

694.363 

36.010 

3 

742  202  203 

.472 

79.227 

*.l*0 

694.363 

36.000 

3 

743  203  ’04 

.472 

79.227 

1.1*0 

694.363 

36.000 

744  204  205 

.972 

’9,’?7 

*.180 

694.363 

36.000 

3 

745  205  206 

.97? 

79. >27 

8. HO 

694.363 

36.000 

3 

746  206  ’07 

.972 

79.227 

*.l*0 

694.363 

36.000 

3 

747  207  200 

.972 

74.227 

0.1*0 

694.363 

36.000 

3 

74*  208  209 

.972 

79. ’27 

1.100 

644.363 

36.010 

3 

744  199  211) 

.37* 

30  .507 

’2.179 

270.000 

46.241 

2 

750  199  211 

.010 

•  10.000 

.110 

10,000 

51.642 

b 

751  700  210 

.010 

10.000 

.  no 

10.000 

5*»o42 

b 

752  200  211 

.756 

61 .614 

45.759 

540.000 

46.291 

1 

753  200  212 

.010 

10.000 

.010 

10.000 

51.642 

b 

754  201  211 

.010 

10 .000 

.010 

10.000 

58.642 

0 

755  201  212 

.756 

61 .614 

45.759 

540.000 

46.291 

1 

756  201  213 

.010 

10.000 

.010 

i'n.uoo 

51.642 

5 

757  202  212 

.010 

10.000 

.010 

10,000 

58.642 

b 

.75*  202  213 

.756 

61.614 

45.759 

540. OCO 

46.201 

t 

759  202  714 

.010 

10.000 

.010 

IO.OCO 

51.642 

5 

760  203  213 

.010 

10.000 

.no 

10.000 

51.642 

5 

761  20.3  714 

.  756 

61.614 

45.759 

540,000 

46.291 

1 

762  203  215 

.010 

10.000 

.010 

10.000 

51 ,642 

5 

763  204  ’14 

.010 

10.000 

.010 

10.000 

51.642 

5 

i 
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THESE  Apr  THE  ELASTIC  N39AL  DISPLACEMENTS 


HOOF 

DELTA  X 

DELTA  Y 

1 

0. 

0. 

2 

0. 

0. 

3 

0. 

8. 

4 

8. 

0. 

5 

0. 

0. 

6 

0. 

0. 

7 

8. 

8. 

ft 

0. 

0. 

9 

0. 

0. 

10 

0. 

8. 

11 

0.  . 

0. 

12 

-3.4?230E-01 

-4.27863E-04 

13 

-3.62242F-01 

-1. 81 915E-0  4 

14 

-3*  84183E-81 

-i. 14901E-04 

19 

-3. 993652-01 

-2.27447E-84 

16 

-4.03493F-01 

-1  .262382-04 

17 

-4.108292-01 

-1.8143 3E-05 

in 

-4.01498E-01 

8  ,933092-05 

19 

-3.99369F-01 

1.910  71E-0  4 

20 

-3. 84150E-01 

2.784852-04 

21 

-1.6230  IE-01 

3  .454932-04 

?4» 

-3.42264E-01 

3.9147  52-04 

23 

-5.86048E-01 

-6.197592-04 

24 

-6.218HE-01 

-5.642562-04 

25 

-6.59906E-01 

-4.71277F-04 

26 

-6.859J7E-Q1 

-1.44548E-04 

27 

-7 .01291 F-gl 

-1.956252-04 

2  * 

-7.05196E-01 

-3.610552-05 

29 

-7.0 1296E-0 1 

1.23413r-04 

30 

-6.859252- 0 1 

2. 7®31 12-34 

31 

“6. 59854F-0 1 

3.939662-04 

3? 

-6.?t9J6r-01 

4.919272-04 

33 

-5.86115F-01 

5.474652-04 

34 

-6. 72590E-O1 

-7.010432-04 

35 

-7.  1605PE-01 

-6.435102-04 

36 

-7.61896E-01 

-S.39103F-04 

37 

-7.906512-81 

-3.967632-04 

3  n 

-8.07782P-01 

-2.382916-04 

39 

-3.120012-01 

-5.260032-05 

40 

-8.077832-01 

1.25  0792-04 

41 

-7. 906312-81 

9158  92 -04 

42 

-7.60839F-01 

4.33721E-04 

43 

-7.162252-01 

5.38079^-04 

44 

-6.727472-01 

5  .97643E-Q4 

45 

-5.936992-01 

-6.543832-04 

46 

-6.34495F-01 

-5.979J9F-04 

47 

-6.757862-01 

-5.003362-04 

4  n 

-7.015732-01 

-3.710492-04 

49 

-7.161912-01 

-2.227732-04 

50 

-7. 19416F-01 

-6.60  8752-05 

51 

-7.161852-01 

9.057352-05 

52 

-7 .015452*01 

2.337682-04 

53 

-6.757312-01 

3.678862— 04 

54 

-6.347152-01 

<•.654052-04 

.55 

-5.939152-01 

5.213712-04 

5b 

-3,863562-01 

-4.644782-04 

57 

-4. 144122-01 

-4.192502-04 

53 

-4.434622-81 

-3.439332-04 

CELTS  7 

0. 

o. 

0. 

0. 

a. 

o. 

o. 

o. 

o. 

o. 

o. 

2. 093702-0? 
2.262432-0?, 

2.402532-0? 

2.501632-02 

2.552JOE-02 

2.56725E-0? 

2*  557252-0? 

2.50159E-C? 

2.402922-02 

2.26275E-0? 

2.09O34E-0? 

7.309592-0? 

7.552672-02 

5*34411 2-02 

5.650942-02 

3.55735E-0? 

5.906612-02 

5.55729C-0? 

5.630642-02 

5. 34403E-0? 

7.554162-02 

7.311812-0? 

1.02H6E-01 

1.I0972P-01 

1. 173642-01 

1.  22516E-01 

1.247542-01 

1.25340F-01 

1.  247512-01 

1.22510E-01 

1.  173642-01 

1.106052-01 

1.  021702-01 

4*  666442-0? 

5.29560F-02 

5. 756092-02 

5.94513E-0? 

5.93847  2-0? 

5.98302S-0? 

9.  95763  E-02 

5.947092.-02 

5.75667E-0? 

5*  30 149E-0? 

4. 67692E-0? 
-1.376582-01 
-1. 43421E-01 
-1.499482-01 


THETA  X 

0. 

0. 

a. 

o. 

o. 

o. 

o. 

o. 

0. 

0. 

0. 

4.64779E-05 
4.33631E-55 
3.39546F-05 
2.  05631E-09 
S.55063E-06 
-1. 11059E-10 
-5.66729F-05 
2.05609E-05 
-3.34916F-05 
-4.32743E-05 
-4.638892-05 
1.51931E-04 
1.432122-04 
1.15165E-04 
7.01262E-05 
2.95436E-05 
-1.637652-09 
-2.99575F-05 
-7.019342-05 
-1.1757UF-04 
-1.47534E-04 
-1.91824E-04 
2.33100C-04 
2.26344E-04 
1 • 6944  3T  -04 

9.24349F-05 
3.580812-05 
-3.362552-08 
-3.59153E-05 
-9. ?439?E-05 
-1.63803E-04 
“ 2*  255042-04 
-2.38355C-04 
1.821232-04 
1.59360E-04 
8.99955F-05 
2.59065E-05 
5*8 1539E-07 
-1.194822-07 
-1. 09941E-06 
-2.532832-05 
-8.59050E-05 
-1.57562E-04 
-1.809612-04 
-1.  56582C-04 
“ 1  *  67553C-04 
-  1. 959552-04 


THETA  y 
-7.76342E-03 
■8.21096,;. 03 
-H.70662E-03 
-9.05075E-03 
-9.29902E-03 
-9.31233E-03 
-9.259l4f-l)3 
■9. 050652*03 
*8.705552-03 
212402-03 
-7 , 7641 1C-03 
-6.690572-03 
-7.09681E-J3 
-7.530152-03 
-7.827772-03 
-8.00359t.-03 
-fl.04863c-03 
-8.003662-03 
-7.  82 764E-03 
-7.529S6E-03 
-7.09816t-r.3 
-b. 691492-03 
-3.844762-03 
-4, 123572-03 
-4.39020E-03 
-4.5blt7E-03 
-4.653302-03 
-4.675622-0? 
-4.65326t-03 
-4.561002-83 
-4.38913^-03 
-4.124772-03 
-3 . 846 l®2-03 
1.52529c-05 
-4. 912282-05 
-8.340092-05 
-7.644202-05 
-6. 155932-05 
-5 .24292E-05 
-6.140212-05 
-7 .b2314E-85 
-8. J4396E-05 
-5.026142-05 
1.337272-05 
4.119502-03 
4.340552-03 
4.56348E-03 
4.76852E-02 
4.88400E-03 
4.925312-03 
4. 884242-03 
4. 768762-03 
4.568382-03 
4  •  33938E-0? 
4.117362-0? 

7 .71976t-03 
8.2S563E-03 
8.737 35E-U3 


Trtl  TA 


121 


59 

63 

61 

63 

63 

66 

65 

66 
67 

65 

69 

70 

71 

72 

73 
76 

75 

76 

77 
7* 
79 
60 
61 
62 
63 
66 
69 

66 
67 
66 
69 

90 

91 

92 

93 
96 

95 

96 

97 
96 
99 

.100 

101 

102 

103 

106 

105 

106 
107 
1 06 

109 

110 
111 
112 
113 
116 

115 

116 
117 
116 

119 

120 
121 
122 


-6.69Jb6F-01 
-6*  61 7692-01 
-6.701075-01 
-6.63773E-01 
-6.  593289-01 
-6.6J6975-91 
-6.166765-01 
-3. 666669-01 
-1.177559-01 
-l.26663f.oi 
-1.33676E-01 
-1.62526E-01 
-l.66161£-'Jl 
-1.65625E-01 
-1.66117E-01 
-1.626779-01 
-1.366169-91 
-1.269669-01 
-1.161269-01 
1 .366229 -0 1 
1.657039-01 

1.51296E-01 
l .556639-0 1 
1.576569-01 
1.5  6560E-0 1 
1.576909-01 

1.557209-01 
1.513679-01 
1.65373F-01 
1.360109-01 
3.15H«9£-01 
•3.361609-01 
3.5036  69-01 
3.S7666E-01 
3. 6361 6E-0 1 
3.6366  09-01 
3.636519-01 
3.577299-01 

3.506399-01 
3.337919-01 
3.156619-01 
3.959019-01 
6.126609-01 
6.765259-01 
6.366679-01 
6.631969-01 
6. 636019-0 1 
6.632369-01 
6. 369519-01 
6«?6637r-oi 
6,121169-01 

3.956766-01 
6.026639-01 
6.150379-01 
6.267039-01 
6,36926  F-0 1 
6.625569-01 
6.61636F-01 
6.625929-01 
6.3696  6F-01 
6.2066  09-01 
6.1<,6569-01 
6  .0 ’0669-, 01 
6. 19  1929-01 


-’.636069-06 
-1.706589-06 
-7  .67675^-05 
2.09263E-05 
t .137559-06 

1  .990632-01, 
2 .69260E-06 
3.16506E-06 
- t . 52369E-06 
-1 .’59239-06 
-1.017379-06 
-6.661035-05 
-7.956019-05 
-7.676969-05 
-7.60  7169-05 
-6.717169-05 
-5. ’31679-05 
-2.628199-05 
-1.603509-06 
2.200699-06 
’.’35059-06 
1.912279-06 
1.231209-06 
3.161969-05 
-7.069919-05 
-1  .7306  19-06 
-2.667379-06 
-3.331639-06 
-3.656629-06 
-3.671709-06 
5.375659-06 
5.216555-06 
6.607919-06 
3.077639-06 
1.312739-06 
-5.569759-05 
-’.62337F-06 
-6.139959-06 
-5.516969-06 
-6.336619-06 
-6.693669-06 
6.333909-06 
6.116819-06 
5 • 192679-06 
3.663269-06 
1.761719-06 
-3.105589-05 
-2.363669-06 
-6  .266639-06 
-5.619609-06 
-6.763629-06 
-6.962169-06 
6.639569-06 
6.716569-06 
6.107669-06 
3.011039-06 
1. 565599-06 
0. 

■1.567209-06 
■2.979919-06 
•6.070955-06 
■6.678595-06 
■6.  79  7559-06 

.’.66  3249-04 


-1.566149-01 

-1. 611615-01 

-1.62671E-01 

-1.611989-01 

-1.568355-01 

-1.  6993bE-01 

-1.633295-01 

-1.376665-01 

-6.59322E-01 

-6.89735E-01 

-5.173929-01 

-5.365675-01 

-5.500659-01 

-5.  560109-01 

-5.500769-01 

-5. 386’IE-Ol 

-5. 173765-01 

-6.  89608E-01 

-6.590769-01 

-8.810655-01 

-9.65150E- 01 

-1. 002535600 

-1.037995*00 

-1.058609600 

-1.063265*00 

-1. 056665*00 

-1.  036O6<r*OO 

-1. 002525*00 

-9.669929-01 

-6. 6  07307-01 

-1.321209*00 

-1.61330E*  00 

-1.699155*00 

-1.  561259*00 

-1.57317^*00 

-1.576109*00 

-1.573229*00 

-1.561329*00 

-1.699155*90 

-1.613125*00 

-1.320665*00 

-1.66263EVOO 

-1.  75  766E*00 

-1.665539*00 

-1.692009*00 

- 1. 9  2526E*  00 

-1.928065*00 

-1 . Q2532C*00 

-1.897069*00 

-1. 865555*00 

-l. 757675*00 

-1.66  2269*00 

-1.  763609*00 

-1. 835529*00 

-1.900669*00 

-1.965205*00 

-1.971725*00 

-1.979005*00 

-1.  971775*00 

-1.965785*00 

-1.900679*00 

-1.  835379*00 

-1.763299+00 

-1.528819*00 


-l.b35565-06 
-3.513339-05 
-2.617925-07 
6.676259-05 
1.638265-06 
1.976365-06 
1. 6992bE-06 
1. 587199-06 
-8.570035-06 
-8.203235-06 
-6.990215-06 
-6.  536609-06 
-2.077739-06 
-6.666629-07 
2.071659-06 
6.533619-06 
7.016065-06 
8.239605-06 
'  8,602275-06 
-1.792609-03 
-1. 756665-03 
-1.306259-03 
-7.620779-06 
-3.  357825-06 
-6.570379-07 
3.36912E-06 
7.626169-06 

1. 30735E-Q3 

1. 7611 3E-03 
1. 79675E-03 
-2.517515-03 
-2.633395-03 
-1.757569-03 
-9.936699-06 
-6.357965-06 
-8.35J35E-07 
6.36627F-06 
9.961665-06 
1. 761255-03 
2.663755-03 
2.522515-03 
-2.6 1 115E-03 
-2.693595-03 
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-1.448?4E*0? 
-1.827256-0? 
-1.041396-0? 
-1.441346*0’ 
-l*  651106-u? 
-2.078596-0? 
-1*  471646*92 
-1.4763bE-0? 
-2.  170746-02 
-t*  392046*0? 
*1  •  27914002 
-7.  25068F-0’ 
-1.360?4F*0? 
-1  •  30972F.-02 
-’. 10/27F-02 
-6*  636226*91 
-2.556?4F-0l 
•6* 128726-31 
-6.433096-01 
-5.312306-01 
-4.22078F-91 
-5*  40127E-01 
-6*  63601E-U1 
-6*  773926-01 
-2.5’30S6-0l 
1.77771F-02 
-7.615746*01 
-7.19646E-0? 
-1.  41642E-0? 
-1. 54155E*0? 
-2.33760E-02 
-1.339006-0? 
-1. 56339F  *02 
-?,?6616c-02 
-1.49146F-02 
-t.5"5496*02 
-2. 16755F-02 
-1.660076-0’ 
-l • 60191F+02 
-2.00S69F-07 
-1 . 056486-02 
-i.6n6oi,:*o’ 
-1.  0604 86 -02 
-7.00440E-0? 
-1 . 60 1 95T  *07 
-1.66405F-0’ 
-2.164296-0? 
-1.  56551E  *02 
-1. 4941 0E-02 
-2. 26507E- 02 


-2.569776-04 
1.44429F-05 
-8.72 647E-01 
-2.06030E-04 
7.06173F-05 
-8.50590F-01 
-2  *  50  7  9  7E-04 
2  .62622E-05 
-0.41293E-01 
-2.098246-04 
-2.64005r-05 
-0,41 192E-Q1 
-1.50053E-04 
-9.54256E-05 
-0. 428232-01 
-9.67101E-05 
-1.49  584F-04 
-0.41161E-01 
-2.76052F-05 
-2.085806-04 
-0.41 179E-01 
2.77791E-05 
-’.503106-04 
-8.5045or-oi 
7.  32215E-05 
-?.8O70Hr-04 
-0.72569r-0t 
3,3303or-05 
-2,558590  04 
-4,5126  9  0  01 
-7.59327F.oi 
-9,61 810E-U3 
-5.28  773F-04 
3. 35024E-0  3 
3.15735F-03 
-3.27770E-03 
-3.55237F-03 
4. 94335E-04 
t .037722-02 
7 .19567003 
-2.11022E-01 
-1.801106-04 
4 , 31 365F-05 
-3.91031E-01 
-2.263736-04 
9.157296-05 
-3.67905F-01 
-1.97633E-04 
6.19249F-05 
-3 i 53926F-01 
-1  .59406E-04 
1 ,73?24r-05 
-3.460796-01 
-1  .013656-04 
-4.665256-05 
-3  ,4499  IF-  o  1 
-4 .79  259T -05 
-1.00056E-04 
-3 ,460 15F-01 
1,61 41QF-05 
-I.50220E-04 
-3.53071F-01 
6. 1328  5F- 05 
-1.97007F-04 


-1.20672F-U2 
8.0Q432E-Q2 
2.950676*00 
-4.866006-02 
9.747Q0E-02 
3.  030056*00 
-3.  36315E-02 
0.304046-02 
3. 04220E*00 
- 1. 43900E-C2 
6.40976E-02 
3. 076236*00 
1.20989E-U2 
3. 697046-02 
3.069616*00 
3.640366-02 
1.347266-02 
3. 076146*00 
6.356006-42 
-1.  306626-02 
3.047O0F*00 
8.279386-02 
-3.  134226-02 
3,033976*00 
9.854706-02 
-4.975  326-07 
7.  94492F*00 
7.94045r-07 
-3.232996-07 
1.406606*00 
-4.1U70F-U2 
-6.502076-02 
-0.61710E-02 
-1.008306-01 
-1.  082116-01 
-1.082776-01 
-1.010>3E-01 
-8.643606-02 
-6. 509306-02 
-4.104966-02 
1.52515E+00 
-1.294586-02 
6. 30916r- 02 
3 .07222E*00 
-2. 828926-02 
7. 876696-02 
3.112016*00 
-1.847196-02 
6. 91164E-02 
3.112366*00 
-4.  689306-03 
5.537436-02 
3.131146*00 
1 , 60299E- 07 
3.469126-02 
3.125016*00 
3.473196-02 
1.649516-07 
3,  1 30  93F  *00 
5.412866-02 
-4.24442F-03 
3.117636*00 
6. 88730F-07 
-1 .820690  02 
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650 

-1.393006*02 

-1.563206*0? 

-3.678426-01 

3. 11459E*UU 

651 

-1.662545-02 

-1.331426-02 

9.340906-05 

7.954436-02 

652 

-4.98410I--03 

-2.335726-0? 

-2.287306-04 

-2.906156-02 

653 

- 1.  360315*02 

-1.541266*0? 

-3.90897E-01 

3.067406*00 

654 

-1.66510E-02 

-1.414156-02 

4.194246-05 

6.262446-02 

655 

-1 • 046016- 02 

-2.142046-0? 

-1 *8641 IF-04 

-1.24970E-0? 

656 

-6.636446*01 

-7.613146*01 

-2.110046-01 

1.526166*00 

657 

2.076106-02 

-2.417426-02 

-1.36778F-03 

-4.106276-02 

656 

-3. 895286- 0? 

-3.  814306-01 

-9.51747E-03 

-6. 277086-02 

655 

-3.420126-01 

-5.264476-01 

-3.  73570F-03 

-8.344866-02 

660 

-5.345356-01 

-5*654326-01 

-7.610466-04 

-9.950886-02 

661 

-5.662246-01 

-5.141566-01 

1.30760E-03 

-1.076406-01 

662 

-5.142106-01 

-5.64625E-01 

-1.400516-03 

-1.07705E-01 

663 

-5.554406-01 

-5.44S46E-01 

5.789826-04 

-9.  969676-02 

664 

-5 .166466-01 

-4.046466-01 

3.665296-03 

-8.  421296-0? 

665 

-3.440716-01 

-2.601856-02 

1.022376-02 

-6.2.84766-02 

666 

-1.621646-07 

2. 08  2566-32 

1.024076-03 

-4.102606-02 

667 

-7.615446*01 

-7.436416*01 

3.079226-02 

1.  584246*00 

666 

-1.143406-02 

-1.67409E-02 

-1.160066-04 

3.902686-03 

665 

-1.466516-02 

-1.27124F-0? 

3.370666-05 

4.753306-02 

670 

-1.541636*02 

-1.442216*02 

1.067746-01 

3.104446*00 

671 

-1.064346-02 

-2.007056-02 

-1.57348E-Q4 

-8.  887326-03 

672 

-1.654446-07 

-1.157676-02 

8.460646-05 

5. 980116-02 

67  3 

-1.567476*07 

-1.505436*0? 

1.253866-01 

3.123326*00 

674 

-1.156556-02 

-1.453146-02 

-1.355006-04 

-3.  123096-03 

675 

-1.656 1SE-02 

-1.264266-02 

6.662706-05 

5.392136-0? 

676 

-1.505566*02 

-1.520666*02 

1.402066-01 

3, 11394E*  00 

677 

-1.346286-02 

-1.864116-02 

-1.04999E-04 

5. 416886-03 

67  6 

-1.611176-02 

-1.346546-0? 

3.660066-05 

4.525086-02 

674 

•1.602006*02 

-1.532716*02 

1.49&95E-01 

3.121706*00 

660 

-1.781156-02 

-1.714706-02 

-5.667056-05 

1.91071E-02 

661 

-1.444116-07 

-1.57073E-02 

-1.306686-05 

3.150256-02 

662 

-1 .608 1 nE*  0? 

-1.537046*02 

1.53416E-01 

3.114196*00 

68  3 

-1.441446-07 

-1.57436E-U2 

-t. 420926-05 

3.114796-02 

684 

-1.3  84596-0? 

-1.71108E-0? 

-5.567456-05 

1.946276-0? 

6  65 

-1.402046*112 

-1.532716*12 

1  .497706-01 

3. 121006*00 

686 

-1.60870002 

-1.400146-0? 

3.544756-05 

4.489406-02 

687 

-1.252366-02 

-1. 860666-0? 

-1.03726004 

5. 7/086E-0J 

668 

-t. 545606*07 

-t,  520646*1)? 

t  ,402226-01 

7.114U7F*00 

664 

-1.65443C-0? 

-1.266566-02 

6.51407F-05 

5.768976-02 

640 

-1.161  8QF-0? 

-1.451726-02 

-1.347036-04 

-2.877336-03 

641 

-1.563376*07 

-1.505766*02 

1.257136-01 

3.125146*00 

6«2 

-1.654466-0? 

-1.151046-0? 

6.67847F-05 

6.079046-02 

643 

-1.074476-0? 

-2.013146-02 

-1.59283C-U4 

-9.391426-03 

644 

-1.641 74E*oo 

-t. 451836*0? 

1.06945F-01 

.7.  105906*00 

645 

-1.464046-07 

-1.273646-02 

3.24342F-05 

4.714036-02 

656 

-1.147136-02 

-t. 165407-0? 

-t.  14726E-04 

4.277246-03 

647 

-7.617446*01 

-7.474196*01 

3.072126-02 

1.58488E*0p 

658 

1.440456-0? 

1.003666-01 

2.76281F-03 

-.7.862126-02 

644 

4.304066-02 

-2.  010776-01 

-8.169926-03 

-5.  57629E-02 

700 

-2.002526-01 

-7.376676-01 

-7.5670*E-03 

-7.39404E-02 

701 

-3. 440056-01 

-4,148266-01 

-2.106146-03 

-3.  775986-02 

702 

-4. 237446-01 

-4 . 074126-01 

5.508806-04 

-9.513526-02 

703 

-4.074666-01 

-4. 24995E-01 

-5.980256-04 

-9.519066-02 

704 

-4.716846-01 

-3. 50605E-01 

1.97442E-03 

-8.792356-02 

705 

-7.443546-01 

-7.  706266-01 

3  .436966-03 

-7.417696-02 

706 

-2.123836-01 

1.050606-01 

8.517846-03 

-5.582476-02 

70  7 

1.123366-01 

1.895626-02 

-2.593886-03 

-3.  876416-02 

706 

-7.4.76716*01 

-6.088712*01 

2.912026-01 

1.585356*00 

70  5 

-5.704076-03 

-1. 27  335E-02 

-5.63  9666-05 

1.  573576-02 

710 

-4.577466*03 

-9.135016-03 

1.266076-05 

3.574236-02 

711 

-1.442246*0? 

-1.212806*02 

6.037686-01 

3.061066*00 

712 

-*.356866-03 

-l  .38384E-0? 

-9.279266-05 

6.970626-03 

713 

-1.121776-02 

*7,  984236-07 

5 .507066-05 

4.711046-02 
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714 

-1«50553£*02 

.220296*02 

6.161926-01 

3.062856*00 

715 

-8  .918376-0? 

-1.342636-02 

-7.587256-05 

9.848376-03 

716 

-1. 132305-02 

-3.57424E-03 

4.  63  7456-05 

3.998126-02 

717 

-1.520776*0? 

-1 . 229056*02 

6.30195E-01 

3.045256*00 

714 

-9.40805F-03 

-1.27752E-0’ 

-5.741346-05 

J .40727F-02 

714 

-1.12404E-0? 

-9. 38945E-Q3 

3,156336-05 

3.552356-02 

720 

-1.53?80£*02 

-1.236736*0’ 

6.39  59  56-01 

3.045806*00 

721 

-1.114406-0? 

-1. 165786-02 

-2.53133E-05 

2.142946-02 

722 

-1.06783E-02 

-1. 06121F-02 

1.129526-06 

2.80479F-02 

723 

-1.537186*0? 

-1. 23915E  *02 

6.438256-01 

3.  036856*90 

724 

-1.U6420E-0? 

-1.064256-07 

1.616656-07 

2.778236-0? 

725 

-1.01726E-0? 

-1.162676-07 

-2.479736-05 

2.  169146-02 

726 

-1.5328  tE*  02 

-l.23671E*0? 

6.396546-01 

3.045426*00 

727 

- t* 121 38F-0? 

-9, 47026E-03 

3.053536-05 

3.6*24946-02 

728 

-9.44366E-03 

-1 . 77437r-02 

-5.627567-05 

1.43425F-02 

729 

-1.520796*02 

-1.2290',E*0? 

6.302156-01 

3. 045246*00 

730 

-1.130?3E-0? 

-8.596696-13 

4.613996-05 

3.976726-02 

731 

-8.95515E-03 

-1.34O67E-07 

-7.591116-05 

1.007396-02 

73? 

-1.505466*0? 

-1.22030E*02 

6 . 160136-01 

3.064466*00 

733 

-1.12388E-02 

-7.93778E-03 

5.62914F-05 

4.333116-02 

734 

-8.36886E-0? 

-1.38804E-0? 

-9.399656-05 

6.735376-03 

735 

-1.491  926  *  02 

-1. ’12396*02 

6.038686-01 

3.057746*00 

736 

-9*828056-03 

-9. 163496-03 

1  ,t3324E-05 

3.502266-02 

737 

-9. 34255E-03 

-1 . 7696  0E-02 

-5.718476-05 

1.603546-02 

738 

-7.43449E*01 

-6. 0871 7E  *01 

2.910566-01 

1.585766*00 

739 

1.23457F-0? 

9.42774E-02 

2.275886-03 

-3. 517146-02 

740 

4.00575F-0? 

-8. 33408E-02 

-4.816626-03 

-4.612867-02 

741 

-8. 83540E-02 

-1.61261E-01 

-2.025196-03 

-5.919566-02 

742 

-1.651 61E-01 

-7.211056-01 

-1.554007-03 

-O.93135C-02 

743 

-’.’3223F-01 

-2. 148006-01 

2.33  947F-04 

-7. 478196-0? 

744 

-2.14850E-01 

-7.238506-01 

-2.000097-04 

-?.4820tF-0? 

745 

-?.’18?3E-01 

-1. 68.34  0F-01 

1.98594F-03 

-6.983146-0? 

746 

-1.O4532F-01 

-9.525*66-0’ 

t  .9’427E-U3 

-4.937»5r-02 

747 

-9.02637F-0? 

9. 76755E-0’ 

5.21915r-03 

-4.616976-0? 

748 

1.01«13E-01 

1.23980E-02 

-2.433767-03 

-3.611277-0? 

749 

-6 .08894E  *  0 1 

*3, 59174E  *0 1 

5.394706-01 

1. 5*2676  *00 

750 

-4.33272E-03 

-6.163266-03 

-3.  l?’437-05 

2. t 9735r-0? 

751 

-3. 884  79F-0? 

-4.678106-03 

-1.352806-05 

2.7’0836-0? 

752 

-1. >12976*02 

-7. 110196*01 

1 .084136*00 

2.932506*00 

75? 

-3.928016-0’ 

-6*  77666F-03 

-4.857736-05 

1.760326-02 

754 

-4.50915E-03 

-3.  921726-03 

1.001746-05 

3. 03992E-02 

755 

- 1. 220  36E*0? 

-7.144576*01 

1,09288F*00 

2.929416*00 

756 

-4.206467-03 

-6.513’4F-03 

-3.9.33696-05 

1. 858416-02 

757 

-4.64578E-0? 

-4.110006-03 

9.136526-06 

2.912746-02 

758 

-1.229126*02 

-7.179356*01 

1.104.306*00 

2.906346*00 

759 

-4.35539E-03 

-6.163966-03 

-3.083936-05 

1.986596-02 

760 

-4.750?8F-Q3 

-4.454546-03 

5.043216-06 

2.757226-02 

761 

-l.’36«16*02 

-7.216116*01 

1,11 2957*  oo 

2.903356*00 

762 

-4.56552E-0? 

-5.576776-03 

-1.724466-05 

2.246926-02 

763 

-4.65756E-03 

-5.045346-03 

-6.612776-06 

2.48260F-02 

764 

-1. ’39236*02 

-7. 274436*01 

1.116386*00 

2.892896+00 

765 

-4.034826-0.3 

-5.070286-03 

-7.425786-06 

2.462616-02 

766 

-4.599336-03 

-5.  551266-03 

-1.640516-05 

2.266546-02 

767 

-1. ’36706*02 

-7.215876*01 

1.112966*00 

2.902936*00 

768 

-4 , 7?672E- 03 

-4.480416-03 

4.200387-06 

2.736316-02 

769 

-4. 38366F-03 

-6. 137006-03 

-2.989926-05 

2.00700E-02 

770 

-1.229136*0? 

-7.179356*01 

1. 104326*00 

2.90625E*00 

771 

-4.624186-03 

-4.130126-03 

8, 425056-06 

2.893736-02 

772 

-4.24033E-03 

-6.439976-03 

-3.83  6256-05 

1.87854E-02 

773 

-1. ’20377*02 

-7.145156*01 

1.09  2787*  00 

2.93102E*00 

774 

-4.50S51E-03 

-3.906166-03 

1.022046-05 

3,039226-02 

775 

-3.929966-03 

-6.790126-03 

-4.877527-05 

1.759626-02 

778 

-1. ’1246660? 

-7.  106997  *01 

1 .083926*00 

2.929??F*00 

777 

-3.9449tr-03 

-4 . 70  70  37  -03 

-t.4703tr-05 

2,69608r-0? 
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77 4 

-4.36009E-03 

-6. 13446E-03 

-3.02590E-05 

2.219747-02 

779 

-6.347447*01 

-3.590977*01 

5.39  3  0  2  7  -  01 

1.572447*00 

710 

4.03457E-O3 

9.40400E-07 

1.491517-04 

-3.06277E-02 

741 

7.585097-03 

-1 • 0797 1E-Q9 

-5. 102277-04 

-3.1.90997-0? 

742 

-t. 303437-02 

-4.64905E-0? 

-9.40339F-04 

-4.127116-02 

7*  7 

-4.461517-0? 

-5.767417-02 

-2.51 639E-04 

-4.62630E-02 

m 

-5.462347-02 

-o. 35994E-07 

-1.779157-04 

-4.494437-0? 

78  5 

-6.363377-0? 

-5.97645E-02 

1.35257E-04 

-4. 900547-02 

786 

-S.79070F-02 

-4.999337-02 

2.225967-04 

-4.632307-02 

787 

-4.92b  19E-0? 

-1.362667-02 

9.620997-04 

-4.137767-0? 

r«* 

- 1. 14241F-02 

9.41232L-03 

5.510117-04 

-3, 49231E-02 

749 

1. 117217-0? 

4. 096767-03 

-1.69 036F-34 

-3.05774E-02 

790 

-3.591447*01 

t. 641597-04 

7.759407-01 

1.432797*00 

791 

4.49290F-04 

-’.03H4E-07 

-4.22904E-05 

2.27523E-0? 

792 

7.53097E-04 

-1 .635407-03 

-4.073047-05 

2.32614E-02 

793 

-7.110497*01 

3.  519977-04 

1.53605E+00 

2.731407*00 

794 

4.23441E-04 

-2. 09999E-07 

-4.294367-05 

2. 176507-02 

795 

6 ,691967- 04 

-1.47450E-03 
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9. 106596*  J.) 

60n 

-5. 64511,2-01 

-2. 26955E4  OQ 

7.76004E*00 

9,  l?*.24f  *00 

603 

-5.288582-01 

-2.  230742  *00 

7.74723C*00 

9#  1U  H/3F  * jO 

6  06 

-4.868J3F-01 

-2,  18372E  *  00 

7.63420F*00 

h*  vtt/srtun 

600 

-3.54258F-01 

-2.07495E400 

/.507lbE*00 

n*34$7ir*oo 

612 

-2.44427F-01 

-’•049852*  00 

7.239252*00 

*H  *65 159 C*  fj 0 

615 

-2.537922-01 

-2.  07985F  *00 

6.90114E4Q0 

4*  36199E  *uO 

626 

-l.768582»00 

-2.622452*00 

8.67748E*00 

9*  36J6QF*  04 

620 

-l.88?75E*00 

-2.67390E400 

8.821U1E*00 

9.453937*00 

63? 

-t.97236E*00 

-2.7 1672E*  00 

8.  98743E  *00 

9.5429 2F *  00 

635 

-».0«678  Eton 

-2. 812862*00 

9 . 087002*  00 

9*  66j66 F* J O' 

634 

-2.157982*00 

-2.859312*00 

9.18l47E*00 

9.7429HE*  30 

641 

-2.19657£*00 

-2.89457E400 

9.19776E400 

9. 75ol6C *  JO 

644 

-2.15«56E*08 

-2.  36025F.4Q0 

9.18145F.400 

9.74241E*'J0 

647 

-2.09664£*00 

-2.A1261F400 

9.08635E400 

9.6611?E*00 

650 

-l.96861E*00 

-2. 71283E*  00 

6,99058E*00 

9*  56597E *00 

653 

-1.08746E*00 

-2.67334E400 

«.81449E*00 

9.44/2  0  F*  0  0 

656 

-1.76374F.400 

-2.61757E400 

8.67844E*O0 

9.36150E*  00 

667 

-2 .46637  2*0  0 

-2. 30940F*  00 

9.51709E*00 

9*  39151E  *00 

670 

-2.625032*00 

‘  -2. 4090 1 E*00 

9 • 503465*00 

9.33J64E*  00 

673 

-2. 70214E*00 

-2.  448  452*00 

9.59820E*00 

9. 39525E*  00 

676 

-2.8 11 16E+00 

-2.52749E400 

9.66305E400 

9*4  3614E+10 

670 

-2.873212*00 

-2.57034E400 

9.73023E*00 

9«44/94F*gO 

642 

-2.9OO29E400 

-2.59883E400 

9.  74216E*00 

9*49344E*00 

645 

8737q£*0Q 

-2, 57077E*  00 

9.729982*00 

9«44/56E*30 

644 

-2.811112*00 

-2.  527412  *00 

9.663332*00 

9.43637F*  00 

601 

-2.69O77E400 

-2.44574E400 

9.600242*00 

9. 39/41E* JO 

604 

-2.62*442*00 

-2.412062*00 

9.49776E*00 

9  •  3?4bbE  *00 

60? 

-’.462432*00 

-2.305752*  00 

9,51702E»00 

9 .391676  «■  00 

70  4 

-2.30673F*00 

-1. 128352  *00 

9.394652*00 

4.4b 197 F*  00 

140 


711 

-2. *73396*00 

-1.  251916*00 

9.266956*00 

9.299706*00 

71* 

-7. 529956*90 

-1.2921*E*00 

9.315S9E*00 

8.71927E+0Q 

717 

-2.61977E+00 

-1.3*3736*00 

9.3*5606*00 

9.-7  2  5  5  76  *  00 

72(1 

-2.670*7E*00 

-1.  77657E*00 

9.399*16*00 

9.357166+00 

723 

-2.701616*00 

-1,399Q9E*00 

9.39199F*00 

9.3*9796*00 

726 

-2.671176*00 

-1.77699E*00 

9.30793E*00 

0.35259F+00 

729 

-2.519976*00 

-1.3*3776*00 

9.3*5666+00 

8.723596*00 

732 

-7.526*76+03 

-1.290076*00 

9.317*66+00 

9.370.796*00 

735 

-2,*7615E*90 

-1.  *5*31E*00 

9.261266+03 

9.293956*00 

7'9 

-2.30769E*00 

-1.1259*F*00 

9.39*196*00 

9. *51996+00 

7*9 

-1.267996*00 

9.14993E-01 

9. .712926+00 

6«56o*56+  03 

757 

-1.*2217E*00 

7. 71301E-01 

0.119996+00 

6.365116+00 

755 

-1. *59975*00 

7.521916-01 

9.1*1976*00 

6.3730*F*)0 

759 

-1.57779F+00 

7.06471F-01 

9.1*210E+00 

6. .75*696+00 

761 

-1.56573E*90 

6.96**7E-01 

9.165006+00 

6.363596*00 

76* 

-l.59977E*00 

6.69956F-01 

9.159616+00 

6.3526*6*00 

767 

-1.56591E+00 

6.95997E-01 

5. 16*396+00 

6.3629*6+00 

779 

-1.5  2  7  96E*00 

7.  063*76-01 

9.1*2016+00 

b.  35-.56C  +oO 

773 

-l.*5691E*0  0 

7.54052E-01 

9.1**136  +  011 

6.37537E*  >0 

776  - 

-l.*2*69E*09 

7. 60756F-01 

0.11*096+00 

6.359656+00 

779 

-1.266976*00 

9.  162076-01 

9.312296*00 

6.566*b6+Q0 

790 

'  6. 506766-01 

3.  790*16*00 

6.302326+00 

3. 79j*5F  *00 

79' 

5.05157E-01 

3.612976*00 

6.099206+00 

.3,612956*  10 

796 

*.799?6E-01 

.7.601716*00 

6.1000*6+00 

3.60l75r*00 

799 

*.325096-01 

3. 570566*00 

6.080996+00 

.7.570536  *o0 

902 

*.095766-01 

3.562706*00 

6.095976+00 

.7.56267F  +  00 

905 

7.91'97E-01 

7. 5*915E*00 

6.0753.76+00 

7.5*912F*00 

909 

*.0919*F-91 

7.56212F*90 

6.095316  +  90 

3  .56209F  *00 

911 

*«37**0£-0 1 

3.57050E*00 

6.080916+00 

3.570*76*00 

01* 

*.9  0597E-0 l 

3.603616*00 

6.102156+00 

3. 60 1666*00 

917 

5.025156-01 

3.609976*00 

6.09*036*00 

•7.60190r*JO 

929 

6.51127E-01 

3.79077E*90 

6  .30160  6  +  0  0 

3. 790266*00 

POOALEN  NUM97R 
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RUNGF-KUTTA  INTERVAL  *  6 


TOTAL  PERCENT  OF 
PERCENT  OF  TOTAL 


LOAO  INCLUDE*)  AT  THIS  TIME  s 
LOAO  TAKFN  IN  THIS  INTERVAL  = 


100 

1? 


NOOE 

1 

Z 

3 


A 

A 

10 

11 

1Z 

13 
1A 
IS 
1C 

17 

18 

14 

zo 

zi 

zz 

Z3 

ZA 

Z5 

76 

Z7 

ZA 

Z4 

-*0 

31 

3Z 

33 

34 

35 

36 

37 
3« 
35 
AO 
A 1 
AZ 
A3 
AA 
AS 
A6 
A7 
AA 


DELTA  X 
0. 

0. 

0. 

a. 

o. 

o. 

o. 

o. 

0. 

0. 

A. 

-3.7S71AE-01 

•A*  059726-0 1 

-A.37Z31F-01 

-A.60AA5r-01 

-A.77034F-01 

-A.A4Z0TE-01 

-4.8T479F-01 

-A.7A34Zr-u; 

-4.571946-01 

-4.3335ZE-0I 

-4.1Z4376-01 

-6.39Z79E-01 

-4.4636  3F-01' 

-7.5T544E-01 

-7.96S97E-01 

-A.Z60A4E-01 

-A. 347466-01 

-A.403T5E-01 

-A.ZA53AE-01 

-7.946716-01 

-7.5ZZA6F-01 

-7.1T666F-01 

-7.17A77C-01 
-7.90500E-01 
-8.617997-01 
-4. 1506AE-01 
-9.51507E-01 
•4.6OZ13F-01 
-4.71A46E-01 
-5.55A54F-01 
-4.21796E'-01 
-8.714847-01 
-A.ZJ36TF-01 
-6.043137-qi 

-6.7349’F-Ol 

-7.AZ73AE-01 

-7.434947-01 


DISPLACEMENTS 
DELTA  Y 

0. 

0. 

0. 

0. 

I). 

0. 

A. 

0. 

0. 

0. 

0. 

-5.05345E-03 
-5.01A63E-0T 
-A.955A0E-03 
-4.A6730E-03 
-4.756906-0 3 
-4'.6-*4?of-03 
-4.510356-0’ 
-4.T96T36-01 
-4 • TO  2796-0  3 
-4.736307-0’ 
-4.143T1E-03 
-9.1T171C-03 
-4.11 134E-03 
-9.04530E-03 
-S.916A5E-03 
-8.7(.47ZE-03 
-8.54A03E-03 
-8  .TAA70E-03 
-8.168066-03 
-8.0Z6A6E-03 
-7.4414AE-03 
-7.90047E-03 
-1.16 114E-02 
-1.16136E-0? 
-1.15535E-0? 

-1 . 14 10AE-0? 
-1.1Z1Z5E-0Z 
-1.09A37E-0Z 
-1.07498E-0? 
-1.05’73E-02 
-1.03721E-02 
-I.OZ7A3E-OZ 
-1 .02416E-0? 

-t.  18347E-0Z 
-1.1839ZE-07 
-1  .17837E-n? 
-t.i6';?3E-n? 


OELTA  7 

0. 

0. 

0. 

0. 

0. 

0. 

0.' 

0. 

0. 

0. 

0. 

Z.1Z  I74E-0? 

2. 41 1 75E-02 

Z.6306AE-0? 

Zi  7A1A1E-0Z 

2. 87S19E- 0? 

Z.9Z15TE-0Z 

Zi  42166C-07 

Z.A7056E-0Z 

Z.76Z61E-0Z 

Z.6176ZE-0Z 

Z.A6Z43E-0? 

7.  607AAE-0Z 

8.6Z361E-0Z 

9.AAZ7ZE-0? 

1.0QX10E-01 

1.  0A1T4E-01 

1.06  080E-01 

1.06ZA7F-01 

li04501E-01 

li 0  059AE-01 

9. 5070ZE-0Z 

Ai  927696-07 

1.00717E-01 

1.  17901E-01 

1..T1759E-01 

li  A166ZE-01 

1.411ZAE-01 

1.51644E-01 

1, 577636-01 

li  50OT0E-01 

1.A5619F-01 

l«3710ZE-0t 

1.Z7407E-01 

Z.Z'ZOOE  07 

3.813337-0,. 

5. 041 357-0? 
S.A5574E-07 


THETA  X 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

8.603A1E-O5 
7 .490706-05 
5.3A837F-05 
3.54677C-05 
Z.0Z39ZE-05 
6.66625E-06 
-7.5Tb97E-Q6 
-Z.A7AA7F-05 
-3.93Z1 16-05 
-4.5SZ90E-05 
-A.5AA65E-05 
Z. 949566-04 
'Z.63694E-0A 
1.977466*04 
1.361A7E-04 
8. 0A9A1E-05 
3.03670E-05 
-Z.13691E-05 

-8, 129036-05 

-1.39Z6ZE-0A 
-1.65090E-04 
-  li  6Z65  7E-U4 
4, 943907-04 

Ai40566E-04 
3.ZA580E-04 
Z. 247126-04 
1  • T6795E- 04 
6.428386-05 
-7.7Z070F-06 
-9.7SZ9ZE-05 
-li9A617E-04 
-Z.6Z1Z6E-0A 
-Z.7Z43ZE-04 
4.36Z0ZE-0A 
3, H96617 -04 
2.785676-04 
I. 750446-94 


THETA  Y 
-8.56328E-03 
-9.171ZAE-U3 
-9.A6Z74E-03 
-li03A73E-0Z 
-1.074796-OZ 
-1.09065C-0? 
-1.0899AE-02 
-1.068426-07 
-1. 03006E-0Z 
-9.77Q40E-03 
-9.309Z5c.-03 
-7.3Z070E-03 
-7.979l9t-03 
.623846-03 
-9.1 00596-03 
-9.A397AE-03 
-9,590/66-0? 
-9.592526-03 
-9.A069AL-03 
-4.0(,ZJ4L-03 
-A. 578406-03 
-A.14Z4ZE-03 
-3.99A7TC-03 
-A  .53Z6ZE-03 
-5. 01096C-U3 
-t>.3f>ZA7E-03 
-5.6OZA5E-0.3 
-5,730446-03 
-5.  765526-03 
-S  .6  854  06-03 

-5,485266-03 
-5.166626-03 
-A .833966-03 
5.369236-04 
4.250686-04 
Z. 718156-04 
1.734436-04 
1.04094C-04 
A. A«773t-il5 
-1.545026-05 
-7 ,841236-05 
-1.1 72506-04 
-8.059296-05 
1.2/854c-05 
5.326146-03 
5.714246-03 
5.48666C-03 
A.Z7Z7JC-03 


49 

-4.21134E-01 

-1.14713E-02 

6.36550E- 02 

1.04769E-04 

6.37506E-03 

90 

-4.45344E-01 

-1.12644E-02 

*  6. 64733E-02 

7.39247E-05 

6.441742-03 

51  / 

-4.49464E-01 

-1.10610E-02 

6.9.U93E-02 

4. 61640E-05 

6.349662-03 

52 

-8.37S34F-01 

-1  .54689E-02 

7.00278E-02 

-2.51960E-07 

6.22534E-03 

53 

-4.09627E-01 

-t.071412-0? 

6. 64650E-02 

-9.13160E-05 

5.94543E-03 

54 

-7.640*12-01 

-1.062122-0? 

6.32356E-02 

-1. 432572-04 

5.651452-03 

55 

-7.19975E-01 

-1.0S835E-02 

5.  57627E-02 

-2, 13649E-04 

5.39130E-03 

56 

-3.45922E-01 

-9.66246E-03 

-2.07497E-01 

-3.448912-05 

9,429532-03 

5? 

-3.96589E-01 

-9.64487F-03 

-2.  06404E-01 

-4. 78299E-05 

1. 014  76E-02 

54 

-4.49721E-01 

-9.597.15E-03 

-2.11027E-01 

-9, 73205E-05 

1 .0  44  74E-02 

59 

-4.86787E-01 

-9.50  705E-03 

-2.  15038E-01 

-1.121192-04 

1. 13468E-02 

60 

-5.12401E-01 

-9.39430E-03 

-2. 18092E-01 

-6. 31 650  F.- 05 

1.17250E-02 

61 

-5.24445E-01 

-9.27076E-Q3 

-2.14007E-01 

2. 264072-05 

1.190242-02 

62 

-5.30070E-01 

-9.1459  7E-03 

-2.16192E-01 

1.26950E-04 

1.1490  7E-02. 

63 

-5.23484F-01 

-9.074752-03 

-2.10006F-01 

2.22796E-04 

1.170102-02 

64 

-5.07372F-01 

-8.92779E-0I 

-2. 01094F-01 

2.51907E-04 

1.12758E-0? 

65 

-4.76431F-01 

-4 .8576  IE-03 

-1.93043E- 01 

2.01664E-04 

1.07078E-02 

66 

-4.46201F-01 

-8.81945E-03. 

-1.86592E-01 

1.76750F-04 

1.00453E-02 

67 

-2.39070E-02 

-5 .56  67  7E-0? 

-5.96128E-01 

-9.24741F-04 

1  .204 12E-02 

60 

-4.9S849E-02 

-5.5629  IE-03 

•6.  29061E-01 

-9.030022-04 

1.307952-02 

69 

-7.93984F-02 

-5.55777E-03 

-6.603562-01 

-8. 31 836E-U4 

1.403512-82 

70 

-9.71244E-0? 

-5.55059E-03 

-6.  471.UE-01 

-6.75  2  05  2-  04 

1.473322-02 

71 

-1.11744E-01 

-S.54  547E-03 

-7 • 04540F- 0 1 

-3. 60549E-04 

1.517452-02 

72 

-t.18059F.01 

-5.54200F-03 

-7. 12916F-01 

-8.74310F-05 

1.538842-02 

73 

-1.22837E-01 

-5.53976E-03 

-7.10594E-01 

2. 03495E-Q4 

1.54220t-07 

74 

-1.714S1E-01 

-5. 53629F-03 

-6.9  7  4  93E-01 

5.239452-04 

1.52357C-0? 

75 

-1.17961F-01 

-5  .53  867E-03 

-6.  73325E-01 

4.09904E-04 

1  ,4791  6!'-02 

76 

-1.06145F-01 

-5.52176E-03 

-6.41650E-0! 

9. 370452-04 

1.403466-07 

77 

-9.71805F-0? 

-5.51119F-03 

-6.07441E-01 

9.604552-04 

1.316192-0  2 

76 

2.74576E-01 

-5.75856E-04 

-l.  09666  0  00 

-2.153421  -03 

1 ,2/2712-07 

79 

2.72714F-01 

-6.S4192E-04 

-1.174S1E*00 

-2.  174U2E-II3 

1.349782-0? 

60 

2.64054E-Q1 

-7.65411E-04 

-1.24844F*  00 

-1.41247F-03 

1.507772-07 

61 

7.6S374F-01 

-4.46525E-04 

-1.3Q272E+00 

-1.25244E-U3 

l. 569692-07 

62 

7.67690F-81 

-1.01904E-03 

- 1. 339322*00 

-7.157472-04 

1.673102-0? 

43 

2.60554E-01 

-1.15420E-03 

-1.354422*00 

-2.266502-04 

1.63483E-02 

44 

2.57549E-01 

-1  .76462E-03 

-1. 35553E»00 

2.69440F-04 

1.64496E-02 

45 

2. 54172F-01 

-1.41600E-03 

-1. 33453E*00 

4.365032-04 

1.67420E-02 

46 

2.44795E-01 

-1.579632-03 

-1.29439E400 

1.47383E-03 

1.595092-02 

47 

2.42005E-01 

-1 .63130E-03 

-1.23016E*00 

1.95736E-03 

1.50917E-02 

64 

2.31415E-01 

-1.70156E-03 

-1. 15464E700 

1.99770E-03 

1.41923E-U2 

49 

4.82437F-01 

3. 47298E-03 

-1.61357E400 

-3.52503E-03 

1.09197E-U2 

90 

4.98457E-01 

.  3.62381E-03 

-1, 74291E*  00 

-3.71444E-03 

1.197612-02 

91 

5.13132E-01 

3.30664E-03 

-1.86859E*00 

-2.83441E-03 

1.29429E-02 

92 

5.201S6E-01 

3. 04 291 E -03 

-1.94649E*00 

- 1. 89384E-03 

1.355522-07 

93 

5.26345F-01 

2.76  244E-0? 

-2.00655E700 

-1.07826E-03 

1.401012-02 

94 

5.25937F-01 

2.50216E-03 

-2.0255ZE400 

-3.75459E-04 

1.41524E-02 

95 

5. 2562 OE- 01 

2.24443E-03 

-2.03341E700 

3.081442-04 

1.427462-02- 

96 

5.10412E-01 

1.99697E-03 

-2, 00217E*  00 

1. 04393E-03 

1.41374E-02 

97 

5.09515F-01 

1.77775E-0? 

-1. 95371E*  00 

2.02124E-03  • 

1.39865E-02 

96 

4.90569E-01 

1.53549E-03 

-l.45596E*00 

2.99831E-03 

1.344766-02- 

99 

4.69383F-01 

1.75653E-07 

-1.75067EfOO 

2.881572-03 

1.29599E-J2 

too 

5.7484 7E-01 

5.61100E-03 

-2.  00915E*00 

-4.59259E-03 

6 . 16001E-03 

101 

S.95780E-01 

5.2(1 322E-03 

-2.17506E700 

-4.63350E-03 

6.456392-03 

102 

6.16367E-01 

4.74362E-03 

-2.32969E100 

-3.575082-03 

fc. 630732-03 

103 

6.20424E-01 

4.37548E-03 

-2. 43129E*  00 

-2.41043E-03 

7.09773E-03 

104 

6.34011E-01 

4.01629E-03 

-2. 50  525  E*0  0 

-1.40752E-03 

7.212672-03 

105 

6.40049F-01 

3. 69483E-03 

-2.53433F700 

-5.  43055F-04 

7.542402-03 

106 

6. 40143E-0 1 

3.34709E-03 

-2. 54497E*  00 

2.70049E-04 

7.449142-03 

t07 

6.33322E-01 

3.09044E-03 

-2.514722*80 

1.135482-03 

7.711252-03 

108 

6.23503E-0 1 

2. 441972-03 

-  7. 461612*00 

2.12724E-03 

7.69159E-03 

109 

6.04576F-01 

2.57175E-03 

-2. 36t 14E*  00 

3.  100972-03 

4.169302-03 

110 

5.84792F-01 

2.36425E-03 

-2.25031E+00 

3.062452-03 

4.5  0  7  7  92-  03 

111 

5.4’175E-01 

2.14126E-03 

-2 . 1  09 19  E*  0  0 

-4.50633E-03 

-2.21234E-03 

112 

9.99742F-01 

l .79606E-03 

-2.  26685F*00 

-4.122842-03 

-4.02994E-03 

143 


113 

6.1755AF-01 

1.375436-03 

-2.349766'.  00 

-3.071556-03 

-5.042606-03 

114 

6.29746^-01 

4.424336-04 

-2.49838E*00 

-2.291616-03 

-5.596126-03 

115 

6. 3*21 *6-0 1 

4. *41956-04 

-2.565247*00 

-1.40710c-03 

-5.8*2*66-0? 

116 

6.408456-01 

0  . 

-2.600446*00 

-5.876456-04 

-5 .9 179*l-03 

117 

6.40597E-01 

-3.550056-05 

-2.608096*00 

1.873736-04 

-5.901346-03 

ll« 

6.3470  AF-0 1 

-9.5**966-05 

-2.5*7236*00 

9.409596-04 

-5.6414*6-03 

111 

6.25431E-01 

-1.745976-04 

-2.541146*00 

1.631*76-03 

-5. 1430*4-03 

120 

6.09541E-01 

-2.642546-04 

-2.  472967*00 

2.061506-03 

-4.7*321t-03 

121 

5.4, *467-01 

-7.4307*6-04 

-2.397206*00 

2.  124676-03 

-3.241416-03 

122 

6. 0  3964E-0 1 

2.461076-07 

-l.«0*’«F*00 

-2,?65t  1 6-03 

-9.447646-03 

123 

6 .73707E-0 1 

2.3680  76-0  3 

-1. *90076*00 

-2. 749626-03 

-1 .1565*6-07 

124 

6.43141E-01 

2.246996-03 

-1.96*456*00 

-2.06004E-03 

-1.267*96-07 

125 

6 .5455*6-01 

2.110927-03 

-2.035067*00 

-1.610096-03 

-1 .*41406-0? 

126 

6. 691776-01 

1.975356-03 

-2.  08  2416*00 

-1. 059046-03 

-1.. 3*4576-02 

127 

6.7242*6-01 

1 .*49206-03 

-2.  1107.37*00 

-4.637290-04 

-1.4050*6-02 

126 

6. 721646-0 1 

1.73*526-03 

-2.116936*00 

8.378546-05 

-1.406046-07 

124 

6.65A27E-01 

1 .649376-03 

-2.104*96*00 

5.718926-04 

-1,. 307556-02 

130 

6.55327C-01 

1 .5*0*46-03 

-2.  077*06*00 

8.944936-04 

-1.3455*6-07 

131 

6.37671^-01 

1.4673*6-03 

-2.0437*7*00 

.9.91219004 

-1.2*2396-07 

132 

6.22276E-01 

1.447776-03 

t2.  007*36*00 

9.890356-04 

-1.209336-07 

133 

7.29399F-01 

-1.00**56-03 

-1.759446*00 

-7.1*5746-04 

-1.373226-0? 

134 

7.59171E-01 

-1.006356-03 

-1.2*7706*00 

-8.476276-04 

-1.454516-02 

135 

7,  90  91  tF-0 1 

-1.010456-03 

-1.  320736*00 

-9. *71596-04 

-1.54  3  1  66  -  02 

136 

*.12424r-01 

-1.023196-03 

-1.  355206*00 

-8.982916-06 

-1.609636-02 

137 

4.263506-01 

-1.0429*6-03 

-1.3*3246*00 

-6.480296-04 

-1.6S2046-0? 

130 

4.31*606-01 

-1 .065536-03 

-1.400*56*00 

-3.242026-04 

-1.67339C-0? 

134 

4.31722E-01 

-1.0*7616-03 

-1.406636*00 

-5.598906-06 

-1.674196-0? 

140 

4 . 237*4E-0 1 

-1.10642p-03 

-1.402366*00 

2.276136-04 

-1.654526-02 

141 

4.093*46-01 

-1.1143*6-03 

-1.392246*00 

3.0990  7E-04 

-1.624766-02 

142 

7.067516-01 

-1.195076-03 

-1.  3*1906*00 

2.635010-04 

-1 .5*0446-0? 

143 

7.6704 *E-01 

-1.174946-03 

-1.373546*00 

2.227446-04 

-1.5*3046-0? 

144 

4.73252F-01 

-S.60747E-03 

-6.506  336-01 

3.741716-05 

-1.409356-02 

145 

1.013"3E*(I" 

-5.5*6306-03 

-6.  44  7946-01 

-2.403390-05 

-1.451616-02 

146 

l. 05S*2r»00 

-5.54  3646-03 

-6.530066-01 

-1.792750-04 

-1.5l530t-02 

147 

1.047416*00 

-5.4770*6-03 

-6.61*757-01 

-2.94*1 3F-04 

-1.566*54-02 

144 

1. 107" 76*00 

-5.74*397-03 

-6.  727446-01 

-2.  931*16-04 

-1 .603166-0? 

144 

t.  11644»00 

-5.1146  06-03 

-6. *16176-01 

-1.968506-04 

- 1.  b?2t  36-0? 

150 

1.117057*00 

-5.231157-03 

-b. *674*6-01 

-*.*69556-05 

-1 .624126-0? 

151 

1.1066*F*00 

-5.15235F-03 

-6.0*9326-01 

-4.25166F-05 

-1.6114  9l-0? 

152 

1.0*744E*on 

-5.U8360E-03 

-6.4117*6-01 

-8.951*16-05 

-1.5«431t-U? 

153 

1.054116*00 

-5.033576-03 

-6.46?*6'-01 

-1.74*7  JF- 04 

-1.5nln?4-02 

154 

1.034166*00 

*5 . 00  24  7E-0  3 

-7.034*06-01 

-2.  2214  /0  04 

-1  ,5.55596-0? 

155 

1.242146*00 

-1.043307-07 

-1.291496-01 

4.453100-04 

-1.16*016-02 

156 

1.331236*00 

-1.041426-02 

-1.  114376-01 

4.479036-04 

-1.200606-0? 

157 

1.3*417F*00 

-1.076146-07 

-4.773276-02  • 

3. 054046-04 

-1.2.31124-0? 

15* 

1 . 4255 IF  *0  0 

-t.0’6047-02 

-9.007376-02 

1.146466-04 

-1.262526-0? 

154 

1.453347*00 

-l  .012*36-07 

-*.*450*6-00 

-1.419276-05 

-1.7*1614-0? 

160 

1.466546*00 

-4 .4*  03SF-G3 

-9.03*306-02 

-4.320326-05 

-1,304446-0? 

161 

1.467376*00 

-4.332456-07 

-9.  414476-02 

-1.263716-04 

-1.308304-02 

167 

1.454406*00 

-9.649726-03 

-9.9*9976-02 

-1.937670-04 

-1.303436-02 

163 

1.431876*00 

-4.541376-03 

-1.  086966-01 

-2.49914F-04 

-1.294246-02 

164 

1.349546*00 

-4.520516-03 

-1.709916-01 

-3.729*96-04 

-1.2*5544-02 

165 

1.371337*00 

-9.4*2176-03 

-1.352956-01 

-4.0097*6-04 

-1  .279316-0? 

166 

1.562136*00 

-1.422026-07 

2.114106-01 

7.  285430-04 

-7.105 126-03 

167 

1.61*427*00 

-1.4217*6-02 

2.369586-01 

6.576096-04 

-7.193036-0? 

16* 

1.676406*00 

-1.416356-02 

2.579116-01 

4.906126-04 

-7.25655E-03 

169 

1.723466*00 

-1.405176-02 

2. 723826-01 

3.0*08.36-04 

-7.341056-03 

170 

1.757256*00 

-1.3*9766-02 

2.805456-01 

1.452926-04 

-7.453*36-03 

171 

t. 774746*00 

-1.372076-02 

2.  *33206-01 

1.168546-05 

-7 • 557576-U3 

172 

1  . 776*5E*  0  0 

-1.354256-02 

2.8161*6-01 

-1.0529  46-04 

-7.623786-03 

173 

1.764346*00 

-1  .338756-0? 

2.7572*6-01 

-2.177366-04 

-7  .662234-03 

174 

1.740437*00 

-1.325656-07 

2.659586-01 

-3.  156036-04 

-7  .698326-03 

175 

176 

1. 701076*00 
1.6*0516*00 

-1.317516-02 

-1.317067-02 

2. 530136-01 
2.396246-01 

-3.693866-04 

-3.7491*6-04 

-7.766436-03 

-7.861*66-03 

144 


m 

178 

179 

190 

191 
197 
193 

m 

199 

196 

197 
199 
199 

190 

191 
197 

193 

194 
199 
196 
1«7 
199 
199 

?no 

201 

70? 

203 

204 

205 
706 
207 
209 
209 
?10 
211 
712 
213 
714 
215 
?16 
717 
719 

219 

220 
271 
222 
273 
224 
725 
226 
227 

229 
279 

230 
731 


1.73496E*0Q 

1.760477*00  ' 

l.M6607*0G 

1.962*07*00 

1.99706E*00 

t.O161SF*00 

1.92073F+00 

1.910497*00 

1.099327*00 

1.66007e>no 

1.839317*00 

1.624047*00 

1.66603^*00 

1.71106F*00 

1.749797*00 

1.779867*00 

1.797699*00 

1. 803347*0  0 

1.796377*00 

1.79031E+00 

1.759617*0  0 

1.742967*00 

1 . 29  337E*  0  0 

1.304fl7F*o0 

1.333307*00 

1.360367*00 

1. 397707*00 

1.399627*00 

1.410567*00 

1.3945?E*00 

1.3961 3F*n  o 

1.371797*00 

1. 367997*00 

7.10106F-01 

7 .158467-01 

7.7990  97*01 

7.435117-01 

7 .549407-01 

7.61991F-01 
7.6*6657-01 
7 .627997-01 
7.5774 9r-01 
7.516397-01 
7.46607F-01 
0. 

0. 

0. 

0. 

0. 

0. 

9. 

0. 

0. 

0. 

0. 


-1  .S9652E-02 

-1.595077-02 

-1.579497-02 

-1.569697-07 

-1.553927-0? 

-1.536747-0? 

-1 .51925E-0? 

-1 .503497-02 

-1 .491117-02 

-1.492757-02 

-1 .47795E-02 

-1.476147-02 

-1  .473127-02 

-1.467677-02 

-1 .4S966E-02 

-1.44663E-02 

-1.432607-02 

-1.41917E-Q2 

-1.404967-02 

-1.39477E-07 

-1.396937-02 

-1.3a  1947-02 

-1.11939E-02 

- 1 . 11439E-Q  2 

-t .109947-0? 

-1.10336F-02 

-1 .09495E-Q  ? 

-1.095177-0? 
-1.0  75007-02 
-1 .06546f-02 
-1.05747C-02 
-1.05141F-02 
-1 .046957-02 
-6. 119197-03 
-5.991557-03 
-5.947247-03 
-5.905597-03 
-5.955007-03 
-5. 797666-03 
-5.737007-03 
-5.679797-03 
-5.677337-0.3 
-5.593037-03 
-5.544917-03 
0. 

0. 

0. 

0. 

0. 

0  . 

0. 

0. 

n. 

0. 

0. 


3.430947-01 

3.671957-01 

3.960947-01 

3,947307-01 

4.0976.37-01 

4.132447-01 

4.133767-01 

4.  094107-01 

4.  020967-01 

3.929717-01 

3.935307-01 

3.022427-01 

3.171957-01 

3.29137E-01 

3.390407-01 

3.445197-01 

3.492237-01 

3.490517-01 

3.471137-01 

3.430007-01 

3.379427-01 

3.330177-01 

1.705197-01 

1.760257-01 

1.911247-01 

1.  95’977-ni 

1.99434r-oi 

1.90347E-01 

1.909327-01 

1.902197-01 

1. 99496E-01 

1.963437-01 

1.944957-01 

4.40027E-02 

4.513507-0? 

4.639937-02 

4.74030F-0? 

4. 914207-02 
4.960197-07 
4, 973697-0? 
4.956707-02 
4.915137-02 
4.766737-02 
4.725497-0? 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


6. 993597-04 
6.091427-04 
4.46923C-04 
3.139127-04 
1. 974417-04 
6. 310277-05 
-5.535047-05 
-1. 614927-04 
-2.302317-04 
-2.629917-04 
-2.561557-04 
4.349967-04 
3.016097-04 
2.990037-04 
2.191247-04 
1.442997-04 
6.359027-05 
-1.752537-05 
-9.950067-05 
-1.. 390957-04 
-1.439397-04 
-1.306107-04 
1. 560597-04 
1.534147-04 
1.  3435S7-04 
1.05391F-J4 
7. 325797-05 

3.572637-05 
-3.199697-06 
-3.779137-05 
-6.06105E-05 
-5.96540^-05 
-4.962727-05 
3.047967-05 
3. 593927-J5 

3.57953E-05 

2.624197-05 

1.520127-05 

9.738077-06 

-1.^25037*06 

-1.014197-05 

-1.552627-05 

-1.44640F-05 

-1.176407-05 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


-i.i«9iit'-u3 
-9.299237-04 
-7.4575/7-04 
-b.33459t-04 
-5.961972-07. 
-6.162307-04 
-6.724967-04 
-7.619177-04 
-8.040542-04 
-1.041977-0.3 
-1.21S14CJ3 
5.046737-03 
5.559977-03 
5.90545i-03 
6.144607-93 
6,296607-03 
6.360942-03 
6. 346937-03 
6.259607-03 
6.114297-03 
5.934677-03 
5,  76470&-0.3 
1.056707-02 
1.104027-0? 
1.140917-0? 
1. lo99?c-92 
l  ,!9l43t-0? 
l.20369t-0? 
1.706717-07 
1.20  04  lt-0? 
1.11/397-0? 
1.1/0517-02 
1.157432-OP 
1 . 4439*7-02 
1.465447-02 
1.497217-0? 
1.526997-02 
1.551167-02 
1.566097-02 
1.571667-02 
1.567292-02 
1.555947-02 
1.540297-02 
1.530622-02 
1.596067-02 
1.5913U-02 
1.62099E-U2 
1.650617-02 
1.675277-02 
1.690647-02 
1 .696767-02 
1.692767-02 
1.6H193F-02 
1.66660c-0? 
1.658407-02 
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MOMENTS 

NO, 

.  MON.  AT  J  TNO 

MO«.  »T  K  fno 

SHEAR 

1 

0. 

0. 

0. 

2 

0. 

0. 

0. 

3 

0. 

0. 

•  0. 

4 

0. 

0. 

0. 

5 

0. 

0. 

0. 

6 

0. 

0. 

0. 

7 

0. 

9. 

0. 

4 

0. 

0. 

0  . 

1 

0. 

0. 

0. 

10 

0. 

0. 

0. 

11 

•4,  361  MF+01 

-1.507886*01 

6.165016-01 

12 

2 .971146-03 

-t.  S7313E-J7 

-7 .750636-05 

13 

-2.25256E-03 

7, 38  4*56 -04 

5. 10072F-05 

14 

1 . 99139F ♦ 02 

-1.64  23  96*92 

-7 . 855626*00 

15 

J.22279F-03 

-2. 05707E -03 

-9.00609C-QS 

16 

-2.tT419E-03 

7.40  41  IE-04 

4. 9157  7E-05 

17 

2.  12207E»02 

-t.68028E*J2 

-7.99901F.00 

n 

3,08645£-0T 

-2.  261166-31 

-9, I?0b2£-G5 

19 

-1.926»f>r-0* 

9.161O0E-04 

4.  04  729F-0S 

?0 

1.«1210E*02 

-l.65270E*9» 

-7./0105E.00 

21 

2.A9?7T£-03 

-2.37243E-0T 

-0.18  130F-OS 

22 

-1.74568F-0? 

1.06O50E-07 

4. 79925F-05 

23 

1.79526E*02 

-1.616316*02 

-7.170S’**00 

24 

?.o3690E-03 

-2.426846-31 

-8. 636156-05 

25 

-1.513AI.F-07 

1. 165196-01 

4.561916-05 

26 

1.42750E»0? 

-t.54517F»07 

-6.05556E.00 

27 

2.43273E-03 

-2.50972E-03 

-O.42927.E-05 

24 

-1.3156AF-0’ 

1. 247916-03 

4.37191E-05 

2« 

1.46976E»02 

-1.468746*0'’ 

-6.74076F.OO 

30 

2. 18421F-03 

-2.6U02C-03 

-0.17016E-05 

31 

- 1 .  06494E-  03 

1.32056E-03 

4.06690E-05 

32 

1.  24X65E*02 

-1 .34804**3’ 

-S.59944E*00 

33 

1.950A9E-OA 

-2.68h97E-0T 

-7.90634E-05 

34 

-0.940  «4F-04 

1.47’8  IE-03 

4.07006F-05 

35 

1. 02994£. 02 

-1.2’642E*0’ 

-4.07492E.00 

36 

1.69177F-0X 

-2 ,6’S/6£-03 

-7.76256F-05 

37 

-7.70737F-04 

1.64  742E-03 

4.1374t*-05 

3  4 

0.23620Ft0l 

-1.10225F+J2 

-4. 16O90E.00 

39 

1.69117E-0T 

-2.57994F-01 

-7.20425F-05 

40 

-9.’44  97E-04 

1.012O6C-0T 

4.66OUE-05 

41 

4.’4409F*01 

-5.501086*01 

-2.10547E.no 

4? 

2.0O1T7F-07 

-1.70U5bE-01 

-5.30 194F- Q  T 

4  3 

-1.601786-01 

-1.35579E-JI 

b  .  03  296  C- 04 

44 

-t.  U971F-01 

-1 .15714F-01 

4,4O945r.0^ 

46 

-1.107016-01 
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